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INTRODUCTION 


The  radar  cross  section  of  jet  engine  configurations  on  modern 
aircraft  is  a  most  perverse  problem  analytically.  An  ability  to  predict 
the  radar  cross  section  of  loaded  cavity  structures,  even  when  such 
structures  are  only  very  crude  models  of  the  actual  engine,  has  been 
seriously  lacking.  At  the  same  time  it  is  well  known  that  above  a 
certain  spectral  range  the  engine  configuration  largely  controls  the 
radar  cross  section  of  the  aircraft  over  a  large  span  of  nose  and  often 
stern  aspects.  Recently,  the  radar  signal  modulation  produced  by  the 
engines  has  become  of  important  analytical  concern.  The  research 


reported  in  this  document  on  the  radar  cross  section  of  a  loaded  cavity 
structure  is  felt  therefore  to  answer  significant  needs  of  the  Air 
Force, 

This  techical  report  is  separated  into  two  parts.  Part  I,  "‘Ray 
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Analysis  of  EM  8ackscatter  from  a  Cavity  Configuration""  contains  the 
major  thrust  of  our  research  and  basically  details  an  asymptotic 
frequency  domain  analysis  of  various  loaded  cavity  structures.  Part  II, 


o 

""Canonical  Response  Waveforms  of  Finite  and  Open  Circular  Waveguides"  is 
essentially  a  time  domain  analysis  of  the  same  types  of  problems.'  Both 
Part  I  and  Part  II  are  self  contained.  ^The  parts  are  interrelated 
however  in  the  sense  that  results  from  one  are  used  in  the  other.  For 
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example.  Part  II  reproduces  a  portion  of  the  exact  Wiener-Hopf 
computations  for  an  open  circular  waveguide.  These  results  were  used  to 
verify  a  portion  of  the  theory  in  Part  I.  Conversely,  frequency 
computations  from  the  analysis  in  Part  I  are  used  in  Part  II  as  part  of 
a  Fourier  synthesis  procedure  to  obtain  canonical  response  waveforms. 

It  is  suggested  that  the  reader  peruse  Part  I  and  Part  II  in  that 

i 

order.  Only  when  the  real  complexity  of  the  problem  is  understood  dan 

the  approximations  used  in  Part  II  be  appreciated.  Both  Part  I  and  Part 

i 

II  have  separate  preliminary  discussions  and  separate  summaries.  A 
final  section  of  this  report  summarizes  our  main  accomplishments  and' 
makes  recommendations  for  future  research. 


RAY  ANALYSIS  OF  EM  BACKSCATTER 
FROM  A  CAVITY  CONFIGURATION 

Ching-Chao  Huang 
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SECTION  1 


PRELIMINARY  DISCUSSION  FOR  PART  I 

The  jet  inlet  structure  is  a  significant  scatterer  that  must  be 
taken  into  account  when  computing  the  radar  cross  section  (RCS)  of 
modern  aircraft.  While  this  is  in  general  a  difficult  problem,  it  is 
strongly  felt  that  this  difficulty  can  he  substantially  reduced  if  one 
proceeds  from  relatively  simple  to  more  complex  inlet  geometries,  and 
develops  a  thorough  understanding  of  the  basic  scattering  mechanisms 
involved  in  such  geometries.  These  analyses  can  then  be  modified  or 
extended  so  that  a  realistic  analytical  model  of  the  actual  inlet  (or 
intake)  should  be  obtainable.  With  the  above  view  in  mind,  the  low 
frequency  RCS  of  a  simplified  jet  inlet  configuration  is  analyzed  in 
this  report.  The  method  of  analysis  employed  here  is  based  on  the 
Uniform  Geometrical  Theory  of  Diffraction  (UTD)  [1]  ray  technique  and 
its  modifications  which  are  required  within  caustic  regions,  together 
with  the  self-consistent  multiple  scattering  method  (MSM).  The  modifi¬ 
cations  of  the  UTD  at  caustics  which  are  incorporated  in  this  work 
involve  the  use  of  equivalent  currents  and  aperture  integrals;  it  is 
noted  that  the  aperture  integrals  referred  to  here  are  more  closely 
associated  with  the  Physical  Theory  of  Diffraction  (PTD)  [24],  However, 
the  PTD  in  general  can  be  shown  to  be  related  to  the  Geometrical  Theory 
of  Diffraction  (GTI) )  [21]  and  its  uniform  version  (UTD)  [1]  if  one 
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evaluates  the  integrals  of  PTD  asymptotically  in  the  high  frequency 
limit.  Due  to  the  latter  relationship  between  GTD  (UTD)  and  PTO,  it  is 
convenient  to  think  of  the  use  of  equivalent  currents  and  aperture 
integrals  simply  as  modifications  of  the  UTO  at  caustics,  rather  than  as 
being  a  part  of  the  PTO. 

The  simplified  jet  inlet  model  chosen  here  consists  of  a 
semi-infinite,  perfectly-conducting  hollow  circular  pipe  (or  waveguide) 
with  a  planar  periodic  blade  structure  placed  within  it  at  a  distance, 
l  ,  from  the  open  end  of  the  pipe  as  illustrated  in  Figure  1.  The  axis 
of  the  blade  structure  coincides  with  the  axis  of  the  circular  pipe. 

The  inlet  model  is  illuminated  by  an  external  electromagnetic  plane  wave 
as  shown  in  Figure  1,  and  the  medium  surrounding  the  inlet  is  assumed  to 
be  free  space.  In  the  present  case,  the  blade  geometry  in  Figure  1  only 
crudely  models  either  the  actual  fan  In  the  case  of  a  turbo-fan  type  jet 
engine,  or  the  first  compressor  ring  of  blades  in  a  conventional 
(non-turbo-fan)  jet  engine.  A  second  configuration  which  is  also 
considered  in  this  work  is  one  which  contains  a  combination  of  two 
planar,  non-identical  sets  of  blades  extremely  close  together  to 
simulate  a  stator-rotor  combination  in  which  the  planar  stator  vanes  are 
fixed  and  the  planar  rotor  blades  are  allowed  to  turn  relative  to  the 
fixed  stator.  The  stator-rotor  jet  Inlet  model  is  illustrated  in  Figure 
2.  An  analysis  of  the  RCS  for  the  problem  in  Figure  2  yields  an 
estimate  of  the  RCS  modulation  resulting  from  the  motion  of  the  rotor 
blades  relative  to  the  stator  vanes.  It  is  noted  that  the  stator-rotor 
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(a)  SIDE  VIEW 


(b)  FRONT  VIEW 


W,  :  WIDTH  OF 

STATOR  VANE 

W2  ;  WIDTH  OF 

ROTOR  BLADE 


(c)  ENLARGED  FRONT  VIEW 


Figure  2.  Geometry  of  a  simple  jet  inlet-model  with  a  planar 
stator-rotor  combination  within.  Stator  and  rotor 
lie  in  the  same  plane  for  convenience  of  analysis. 
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blade  configuration  of  Figure  2  serves  to  improve  the  basic  model  of 
Figure  1  because  a  set  of  stator  vanes  is  generally  present  immediately 
behind  the  fan  in  the  case  of  a  turbo-fan  type  jet  engine,  and  they  are 
also  generally  present  immediately  in  front  of  the  first  compressor  ring 
of  blades  in  a  conventional  (or  non-turbo  fan)  jet  engine.  Additional 
problems  which  are  also  considered  in  this  work  deal  with  the  RCS 
analysis  of  finite  length  hollow  circular  cylinders  open  at  the  front 
end,  but  either  closed  (shorted),  or  open  at  the  back  end  as  shown  in 
Figures  3  and  4,  respectively.  The  solutions  to  the  latter  problems  are 
of  interest  in  that  they  allow  a  useful  comparison  of  the  effect  on  the 
RCS  resulting  from  the  scattering  by  the  blade  structure  in  the 
configurations  of  Figures  1  and  2,  as  opposed  to  the  short  and  open 
terminations  of  Figures  3  and  4.  In  the  present  work,  the  RCS  is  of 
interest  for  the  range  of  angles  0  which  satisfy  6  <  n/2;  the  angle  0  is 
illustrated  in  Figures  1-4. 

Some  of  the  earlier  published  work  available  in  the  open  literature 
on  the  types  of  problems  considered  in  this  report  is  briefly  reviewed. 
The  electromagnetic  (EM)  radiation  from  open-ended,  semi-infinite 
circular  and  rectangular  waveguides  was  treated  by  Chu  [2]  via  a 
Kirchhoff  approximation  for  the  fields  in  the  aperture  at  the  open  end. 

A  formally  exact  solution  for  the  problems  of  radiation  and/or 
scattering  from  an  open-ended,  semi-infinite,  hollow  circular  pipe  was 
obtained  via  the  Wiener-Hopf  procedure  by  Levine  and  Schwinger  [3], 

Jones  [4],  and  Noble  [5]  for  the  acoustic  case,  and  by  Pearson  [61  for 
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the  perfectly-conducting  EM  case.  Later  Weinstein  [7]  treated  in  great 
detail  both,  the  acoustic,  and  EM  problems  of  radiation  and  scattering 
by  an  open-ended  semi -infinite  hollow  circular  pipe  via  the  Weiner-Hopf 
technique;  furthermore,  Weinstein  discussed  the  limitations  of  the 
Kirchhoff  method  for  solving  these  problems.  Subsequently  Bowman  [81, 
Lee  et  al  [9],  Mittra  et  al  [10],  Chuang  et  al  [11],  and  Johnson  and 
Moffatt  [12]  essentially  dealt  with  improving  the  numerical  efficiency 
and  accuracy  of  the  asymptotic  high  frequency  approximation  to  the 
Wiener-Hopf  factors  given  by  Weinstein  in  his  Wiener-Hopf  solution  for 
the  EM  radiation  and  scattering  from  an  open-ended  circular  waveguide 
conf iguration.  Witt  and  Price  [13]  analyzed  the  EM  scattering  by  finite 
length  perfectly-conducting  hollow  circular  and  rectangular  pipes  using 
a  Kirchhoff  type  approximati on.  However,  their  procedure  appears  to  be 
incorrect  because  the  dominant  term  in  their  solution  which  pertains  to 
the  scattering  by  just  the  open  front  end  exhibits  the  wrong  frequency 
dependence  for  the  circular  pipe  geometry  when  compared  with  an  accurate 
asymptotic  approximation  of  the  Wiener-Hopf  solution  for  the  same 
problem  [12],  Moll  and  Seecamp  [14]  analyzed  the  RCS  of  an  inlet 
geometry  modeled  by  an  open-ended,  semi-infinite  hollow 
perfectly-conducting  circular  pipe  with  a  planar  blade  structure  inside; 
their  blade  model  is  somewhat  similar,  but  not  identical  to  the  one 
chosen  in  the  present  work  (as  illustrated  in  Figures  1  and  2). 

However,  Moll  and  Seecamp  [14]  employ  the  procedure  of  Witt  and  Price 
[13]  which  appears  to  be  incorrect  as  mentioned  earlier;  furthermore, 
they  include  additional  approximations  (e.g.,  they  neglect  all  TM  modes 
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inside  the  duct)  which  appears  to  lead  to  additional  errors.  Kao  [IS] 
presents  a  numerical  type  solution  for  the  scattering  by  finite  length 
cylinders;  however,  his  work  is  restricted  to  broadside  incidence  on  the 
cylinder.  Finally,  it  may  be  remarked  that  Mittra  et  al  [10]  also 
consider  the  effect  on  the  RCS  resulting  from  a  uniform  planar 
termination  inside  the  semi -inf ini te  hollow  circular  pipe  in  which  the 
planar  termination  is  characteri zed  by  a  surface  impedance  type  boundary 
condition.  Their  analysis  of  this  problem  is  based  on  combining  their 
Wiener-Hopf  solution  with  the  Generalized  Scattering  Matrix  Technique 
(GSMT).  It  is  noted  that  the  MSM  based  procedure  employed  in  the 
present  work  is  essentially  the  same  as  the  GSMT  [16,17].  Both,  the 
GSMT  and  the  self-consistent  MSM  allow  one  to  account  for  the  multiple 
wave  interactions  between  the  open  end  of  the  circular  pipe  and  the 
termination  (or  discontinuity)  placed  within  the  circular  pipe.  A  ray 
analysis  of  the  EM  radiation  from  an  open  ended  circular  waveguide  has 
been  presented  by  Narasimhan  [18];  whereas,  Felsen  and  Yee  [19]  have 
performed  a  ray  analysis  of  the  acoustic  modal  reflection  coefficient 
associated  with  an  open  ended  rigid  circular  pipe  which  is  excited  from 
within.  The  ray  analyses  of  Narasimhan  [18],  and  Felsen  and  Yee  [19] 
will  be  briefly  reviewed  in  Section  3,  wherein  their  work  will  also  be 
compared  and  contrasted  with  the  present  UTD  ray  analysis  and  it's 
modifications  for  the  EM  radiation  and  reflection  from  an  open-ended 
circular  waveguide.  The  latter  analysis  is  necessary  in  the  present 
development  of  the  explicit  expressions  for  the  RCS  of  the  geometries  in 
Figures  1-4. 
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In  the  present  work,  a  ray  analysis  based  on  the  Uniform 
Geometrical  Theory  of  Diffraction  (UTD)  and  its  modifications  is 


combined  with  the  multiple  scattering  method  (MSM)  to  arrive  at  a  very 
efficient  solution  for  evaluating  the  RCS  of  the  inlet  configurations  in 
Figures  1-4.  This  method  of  analysis  will  simply  be  abbreviated  as 
UTD-MSM,  for  convenience.  Basically,  the  MSM  allows  one  to 
systematical ly  take  into  account  all  interactions  between  the  open  front 
end  of  the  pipe  and  the  termination  (or  discontinuity)  at  the  back  end 
located  a  distance  L  from  the  open  end  as  shown  in  Figures  1-4.  In  the 
MSM,  the  wave  interactions  are  described  by  "scattering  matrices" 
pertaining  to  the  canonical  scattering  events  which  occur  at  the  front 
(Z=0)  and  back  (Z=-l).  These  multiple  interactions  can  be  summed  in  a 
closed  form  via  a  self-consistent  procedure.  The  elements  of  the 
scattering  matrices  in  the  MSM  analysis  are  obtained  here  by  using  the 
UTO  and  its  appropriate  modifications  which  are  required  within  caustic 
regions,  and  within  the  confluence  of  caustic  and  ray  optical  shadow 
boundary  transition  regions.  It  is  important  to  note  that  the 
conventional  waveguide  modal  fields  in  the  pipe  region  need  to  be 
converted  into  a  set  of  equivalent  rays  near  the  edge  (or  rim)  of  the 
open  end  in  order  to  systematically  employ  the  UTD  ray  technique  for 
obtaining  the  elenents  of  the  scattering  matrices.  The  UTD  ray  analysis 
provides  simpie  and  accurate  expressions  for  the  elements  of  the 
scattering  matrices,  and  hence  for  the  scattered  and  diffracted  fields. 
The  UTD  analysis  also  provides  a  physical  description  for  the  scattering 
processes  in  terms  of  rays. 
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These  UTD  based  scattering  matrix  elements  for  the  EM  radiation  and 
scattering  by  an  open-ended  semi -inf inite  hollow  circular  pipe  are  found 
to  agree  extremely  well  with  those  obtained  from  the  exact  Wiener-Hopf 
solution  for  the  same  problem;  the  Wiener-Hopf  solution  is  available  in 
the  open  literature  [7,9,10,12]  as  mentioned  earlier.  This  agreement 
holds  up  even  for  the  dominant  mode  regime  in  the  circular  pipe.  A 
distinct  advantage  of  the  UTD  solution,  besides  being  physically 
appealing  and  accurate,  is  that  it  Is  far  less  complicated  to  use  than 
the  Wiener-Hopf  solution.  Furthermore,  the  present  UTD  analysis  can  be 
extended  to  jet  intake  or  inlet  shapes  other  than  those  which  can  be 
analyzed  by  the  Wiener-Hopf  method.  Nevertheless,  the  Wiener-Hopf 
solution  for  the  hollow,  semi-infinite  circular  pipe  is  very  useful  in 
that  it  is  a  formally  exact  solution,  and  it  therefore  provides  an 
important  check  on  the  approximate  UTD  ray  solution. 

A  general  description  of  the  self-consistent  MSM  formulation  of  the 
solutions  to  the  RCS  problems  shown  in  Figures  1-4  are  described  in 
Section  2.  The  elements  of  the  various  scattering  matrices  required  in 
the  MSM  are  obtained  via  the  UTD  and  its  appropriate  modifications  in 
Section  3.  Several  interesting  numerical  results  based  on  the  UTD-MSM 
solutions  for  the  problems  in  Figures  1-4  are  then  presented  and 
discussed  in  Section  4.  Here,  the  UTD-MSM  solution  therefore  implies 
that  the  backscattered  field  or  the  RCS  of  any  of  the  problems  in 
Figures  1-4  is  obtained  via  the  self-consistent  MSM  based  formulation 
which  contains  scattering  matrices  whose  elements  are  obtained  via  the 
UTD  ray  technique  and  its  appropriate  modifications.  In  addition,  some 
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topics  worthy  of  future  investigation  which  would  serve  to  improve  the 
jet  inlet  model  of  Figures  1  and  2,  and  to  extend  the  present  analysis 
to  inlet  shapes  other  than  the  ones  considered  in  Figures  1  and  2,  are 
also  briefly  discussed. 


SECTION  2 


SELF-CONSISTENT  MSM  FORMULATION 

As  mentioned  in  Section  1,  the  method  of  RCS  analysis  of  the 
configurations  in  Figures  1-4  is  based  primarily  on  the  Uniform 
Geometrical  Theory  of  Diffraction  (UTD)  ray  technique  and  it's 
modifications  which  are  used  in  conjunction  with  the  multiple  scattering 
method  (MSM).  In  the  MSM  based  analysis,  the  basic  scattering 
mechanisms  are  "isolated",  and  "identified"  as  being  associated  with  the 
scattering  junctions.  For  example,  these  junctions  are  marked  (I)  and 
(II)  in  Figures  3-5.  The  multiple  scattering  between  the  junctions  is 
calculated  via  a  self  consistent  procedure.  Such  a  procedure  requires  a 
knowledge  of  the  generalized  scattering  matrices  [Sji],  [S^],  [S21], 
[S22],  and  [Sb]  for  the  scattering  junctions  (I)  and  (II).  These 
generalized  scattering  matrices  [16,17]  are  directly  associated  with  the 
canonical  scattering  events  shown  in  Figures  6,  7,  and  8.  The 
scattering  matrix  in  microwave  circuit  analysis  [20]  arises  in  the 
description  of  interior  region  scattering  (e.g.,  in  waveguides)  and  it 
is  defined  in  terms  of  interior  propagating  modes;  whereas,  the 
polarization  scattering  matrix  [31,32,33]  is  defined  for  exterior  region 
scattering.  The  generalized  scattering  matrix  discussed  here  extends 


Figure  5 


Scattering  junctions  I  and  II  for  the  geometries  in 
Figures  1-4. 


gure  6.  Canonical  scattering  problem  for  determining  Sii  and  S?i. 
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Figure  7.  Canonical  radiation  problem  for  determining  2  and 


Figure  8.  Canonical  problem  for  determining  associated  with 
the  configuration  in  Figures  1-4. 
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the  concept  of  microwave  circuit  scattering  matrix  to  include  evanescent 
modes,  and  to  also  include  the  polarization  scattering  matrix  for 
exterior  scattering.  Thus,  the  generalized  scattering  matrix  is 
essential  for  solving  scattering  problems  which  involve  a  coupling 
between  interior  and  exterior  regioqs  as  in  the  problem  of  the 
backscatter  from  jet  inlet  configurations.  The  elements  of  these 
generalized  scattering  matrices  which  will  hence  forth  be  referred  to 
simply  as  the  scattering  matrices  are  defined  next  in  the  paragraphs  to 
fol low. 

_i 

Let  E  denote  the  electric  field  intensity  of  the  incident 
electromagnetic  plane  wave.  Using  the  coordinates  given  in  Figures  1-4, 
one  may  express  E"*  as 

E1  =  (0  El  ♦  l  E1)  ejk(*  Sin9  +  2  C0S  0>  -  (1) 

'  u  <J> 

i  i  .  « 

where  E0  and  ^ 4>  are  the  9  and  $  polarized  components  of  F1  , 
respectively  at  z=0,  and  k  denotes  the  free  space  wave  number.  The  unit 

A  A 

vectors  0  and  z  lie  in  the  plane  of  incidence  defined  by  the  cylinder 

axis  (z)  and  the  incident  ray  direction.  The  unit  vector  $  is 

orthogonal  to  the  plane  of  incidence  as  shown  in  Figures  1-4.  Likewise, 

iet  T$  denote  the  total  far-zone  backscattered  electric  field  and  ^ 

o 

denote  the  far-zone  electric  field  backscattered  from  only  the  open 
front  end  (at  z=0)  of  the  configurations  in  Figures  1-4.  One  can  also 

A  A 

express  Fs  and  in  terms  of  their  9  and  $  components  as 
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One  may  now  describe  the  scatteirng  matrix  [Sji].  The  scattering 

matrix  [S^]  relates  the  field  p  to  Ft  (at  z=0)  as  follows. 

o 


(4) 


where  R  denotes  the  distance  from  the  origin  to  the  far  zone  observation 
point  in  the  R  direction  as  shown  in  Figures  1-4.  Thus, 
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The  scattering  matrix  [S12]  converts  the  circular  waveguide 
(intake)  modal  fields  that  are  incident  on  the  opening  (at  z=0)  from  the 
region  z<0  into  the  fields  radiated  by  these  modes  from  the  open  end. 

The  modal  electric  field  Fg  within  the  circular  waveguide  region  (z<0), 
i.e.,  within  the  Intake  region,  may  be  represented  in  the  usual  manner 
by  [20]  , 


where 


r9-n 

m  n 

e'nm  and  ¥nm  denote  the  transverse  (to  z)  electric  vector  mode 
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functions  for  the  TEnm  and  TMnm  modes,  respectively.  Likewise,  0' 

nm 

and  8nm  denote  the  propagation  constants  of  the  TEnm  and  TMnm  modes, 

respectively.  As  usual,  n  and  m  denote  the  modal  indices  associated 

with  the  ci rcumferential  and  radial  variations  of  the  modes  (or 

eigenfunctions).  Also,  e^m  denotes  the  z  directed  TMnm  electric  mode 

± 

field.  The  corresponding  magnetic  field  TTg  within  the  circular  guide 


_ ±  t  +i8'z  ±  +i R  z 

=  1  l  ( ^nm  +  ^nmz)  nm  +  ^nm  (^nm)6  nm  ]  •  (7) 


The  superscripts  ±  in  Eqs.  (6)-{7)  refer  to  modes  propagating  in  the  ±z 

A  I 

direction.  To  be  specific,  let  behave  as  fcos  n * '  'i  in  the 

'-sin  n<(>'  > 

A  I 

waveguide  aperture  at  z=0  (see  Figure  1-4).  Then  <t>*  behaves  as 

A  A 

/-sin  n 4> *  ^ .  likewise  <f> ' ‘^nm  behaves  as  f“sin  n*'-,  an(j  4,'.^ 

COS  n •}>  J  *  rnc  > 


cos  n<{> 


behaves  as  rcos  In  Appendix  A,  a  list  of  the  circular  waveguide 

'sir  n4>  > 

modes  and  also  the  values  of  ka  at  mode  cutoffs  are  given  in 
chronological  order. 

Let  T@  and  T^,  denote  the  0  and  $  directed  components  of  the 
electric  field  radiated  by  the  waveguide  fields  (Tg,  TTg )  after  they 
impinge  upon  the  open  end  (z=0)  from  the  region  z<0.  Then 


[Sl?l 


e-jkR 
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2x«  jV  3 
nmJ 


where  a  typical  subscript  of  the  type  MxN  on  the  matrices  implies  M  rows 
and  N  columns,  and 


w  * 


The  scattering  matrix  [S21  ]  describes  the  transformation  or 
coupling  of  the  incident  plane  wave  field  into  the  waveguide  modes  as 
illustrated  in  Figure  6.  Symbolically,  [S21 D  is  defined  by  the  relation 
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lBnJ 

=  [$21 1 

ri 

“Xl 

«x2 

Ei> 

_  _ 

(ID 


2  x  1 


where  Anm  and  Rnm  denote  the  amplitudes  of  the  TEnm  and  TMnm  modes 

transmitted  into  the  waveguide  (i.e.,  in  the  -z  direction)  by  the 

i  i A  i A 

external  plane  wave  field  T  =  ^q0  +  which  is  incident  on  the  open 

end  of  the  semi-infinite  circular  pipe.  From  Eq.  (11),  it  follows  that 
[S21 ]  is  given  by 


[s21 


[ S 9nm  ] 

[S<t>nm] 

[ 5 0nm  ] 

[  S  .{inm  ] 

(12) 


«x2 


Clearly  the  problem  of  determining  [S21]  is  the  reciprocal  of  the 
problem  of  determining  [S12]  (see  Figures  6  and  7).  Thus,  a  knowldege 
if  [Si 2];  i.e.,  S gnm,  S$nm,  Sonm.  and  s4>nm.  together  with  the  use  of 

the  reciprocity  theorem  for  electromagnetic  fields  allows  one  to 
calculate  Senin,  S^,,,,  S0nm,  and  S^nm,  and  to  hence  obtain  [S21]  in 
terms  of  [S 1 2 ] •  ln  short,  the  elements  of  [S2 1 ]  are  simply  related  to 

n 
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the  elements  of  [Si2]^5  here  T  denotes  the  transpose  matrix  operator. 
The  precise  relation  between  [S21]  and  [S12]  will  be  discussed  in 
Section  3. 


The  scattering  matrix  [S22]  is  a  modal  reflection  coefficient 
matrix  which  is  associated  with  the  interaction  shown  in  Figure  7.  In 
particular,  the  elements  of  [$22!  describe  the  reflection  coefficients 
associated  with  TEnm  and  TMnm  modes  reflected  back  from  the  opening  (at 
z=0)  into  the  semi -inf inite  waveguide  region  (z<0)  when  either  a  TEnm, 
or  a  TMnm  mode  is  incident  on  the  open  end  from  within  the  waveguide 
region  .  The  matrix  [S22]  is  defined  symbolically  by  the  relation 


lAnm) 

=  [S22] 

[A+  ] 

L  pqJ 

JBn'J 

®X“ 

“Xl 

[b+  ] 

L  p^U 

(13) 
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where,  it  is  noted  that  [S22]  can  be  further  symbolically  expressed  as. 


[S22I  = 


fRhh  1  ! 
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[Rhe  ] 

1  nm;pq  1 

[Reh  ]  | 

1  nm;pq  J  1 

[Ree  ] 

1  nm;pq  J 

00  x  00 

• 

(14) 

The  meaning  of  R 


he 

nm;pq : 


for  example,  in  the  above  equation  is  the 


following.  A  TMpq  (or  e)  type  mode  with  amplitude  RDa  which  is  incident 


pq 


at  the  open  end  (z=0)  is  partly  transformed  (or  coupled)  into  a 

he  + 

reflected  TEnffl  (or  h)  type  mode  with  amplitude  Anm  =Rnm;pqRpq* 


The  scattering  matrix  [S(j],  like  [$22]»  is  also  a  reflection  type 
matrix  associated  with  the  discontinuity  at  junction  (2)  in  Figure  8. 
From  Eqs.  (13)-(14),  one  may  therefore  write 
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~he  he 

ir  which  Rnm;pq,  for  example,  has  the  same  meaning  as  Rnm;pq  of 

,  ,  ~he 

tq.  (14).  In  general,  however,  the  values  of  Rnm;pq  are  quite  different 

he 

from  those  of  Rnm;pq* 

At  a  given  operating  frequency,  the  inlet  duct  or  waveguide  region 
in  Figures  1-4  can  support  a  finite  number  of  propagating  modes,  and  an 
infinite  number  of  evanescent  (or  non-propagating)  modes.  Consequently, 
the  matrices  associated  with  [Si?],  [S21],  [S22L  and  [Sbl  are  of 
infinite  order  to  include  the  infinite  number  of  evanescent  modes.  It 
may  be  remarked  tnat  although  the  matrices  [S12L  C$2l]>  [$22].  and 
are  of  infinite  order  in  a  formal  sense,  one  needs  to  retain  only  a 


finite  number  of  the  elements  of  these  scattering  matrices  in  practice 
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because  the  distance  "L"  shown  in  Figures  1-4  is  generally  large  enough 
to  where  the  infinite  number  of  ''evanescent"  waveguide  modes  generated 
at  junction  (I)  are  not  observed  at  junction  (II),  and  vice  versa.  The 
finite  number  of  elements  retained  in  practice  thus  correspond  to  only 
the  finite  number  of  all  the  propagating  (non-evanescent )  modes  which 
can  exist  in  the  waveguide  or  the  cavity  region.  If  the  distance  L  in 
Figures  1-4  is  small  enough  so  that  lower  order  evanscent  modes  become 
important,  then  one  must  include  these  modes  but  still  ignore  all  the 
higher  order  evanescent  modes  in  practice  since  their  contribution  must 
become  vanishingly  small.  In  any  case,  one  retains  only  a  finite  number 
of  elements  in  the  scattering  matrices.  In  the  present  work,  "L"  is 

t 

chosen  such  that  the  contribution  from  all  the  evanescent  waveguide 
modes  can  be  ignored. 

Referring  to  Figure  9,  it  is  seen  that  the  incident  electric  field 
Ei 

[E1]=[  9]  is  scattered  by  junction  (I)  to  produce  a  scattered  field 

E1 

$ 

c 

[Eq]  in  the  region  A,  then 

CE*]  =  [Sn]  [E1]  •  07) 

Part  of  the  incident  field  which  is  scattered  by  (I)  into  region  R 
becomes  incident  at  junction  (II)  from  which  it  is  subsequently 
reflected;  this  reflected  field  is  incident  back  at  (I)  where  it 
undergoes  further  scattering  into  regions  A  and  B,  and  so  on,  thereby 
giving  rise  to  multiple  wave  interactions  between  junctions  (I)  and  (II). 
The  fields  resulting  from  these  multiple  interactions  may  be  expressed 
in  a  convergent  Neumann  series  as  done  by  Pace  and  Mittra  [16]; 


REGION  (a) 


REGION  (S) 


Figure  9.  Multiple  interactions  between  junctions  (I)  and  (II)  for 
calculating  the  total  backscattered  field. 


however,  an  alternate  procedure  based  on  a  self-consistent  method  leads 

directly  to  the  same  result.  The  latter,  i.e.,  the  self-consistent  MSM 

method  for  summing  up  the  multiple  interactions ,  will  be  used  in  this 

analysis.  Let  [Ej?]  represent  the  net  value  of  the  field  incident  at 

(I)  from  (II)  after  taking  all  the  multiple  interactions  into 

consideration.  Similarly,  let  [E21]  represent  the  net  value  of  the 

field  incident  at  (II)  from  (I)  after  taking  all  the  multiple 

interactions  into  account.  Then  the  total  scattered  field  in  region  A 

denoted  by  [Es]  consists  of  a  superposition  of  the  field  [Es],  and  [Es], 

0  m 

where  [Es]  is  the  field  scattered  into  A  when  [E12]  is  incident  on  (1). 
m 

Thus, 


[ES]  =  [ES]  +  [Es]  .  (18) 

0  m 

where  [Es]  may  be  expressed  as 
m 

[Es]  -  [Si2]  CEl2l  •  (19) 

m  ^ 

The  expressions  for  [E12]  and  [Eji]  are  given  in  terms  of  [S21]  and 
[S22]  by 

[E2ll  =  [S2i][Ei]  +  [S2Z][Ei2]  .  (20) 


[El2]  -  [P][Sb][P][E2l] 


(21) 


where  L"Sb]  's  the  generalized  matrix  corresponding  to  the  reflection 
coefficient  at  junction  (II)  as  indicated  earlier,  and  [P]  is  a  diagonal 
matrix  accounting  for  the  phase  delay  in  the  propagation  path  L. 
Eliminating  [El?]  between  Eqs.  (20)  and  (21)  yields 

([lMS22][P][Sb][P])  [E2i]  =  [S2i][Ei]  .  (22) 

where  [I]  is  an  infinite-order  identity  matrix.  Hence,  it  follows  that 

[E2j]  =  ([I]-[S22][P][Sb][P])-l  [S21][Ei]  .  (23) 

From  Eqs.  (21)  and  (23),  [Ei2]  becomes 

[Ei2]  =  [P][Sb][P]  ([I]-CS22][P]LSb][P])-l  [S21][Ei]  .  (24) 

Incorporating  Eq.  (24)  into  (19)  yields 

[esJ  =  [Si2][P][Sb][P]  (Cl]-[S22][P][Sb][P])-l  [S21][Ei]  e-jk«  . 

R 

(25) 

Finally,  combining  Eqs.  (17)  and  (25)  according  to  Eq.  (18)  yields  the 
sel f-cons istent  expression  for  the  total  scattered  field  [Es]  as 

[E^]  -  {rSii>[Si2j[P][Sb]CP]([lMS22][P]CSb][P])-l  [S2i]}  [Ei]  e-J^  . 

R 

(26) 

In  the  soecial  case  of  the  geometry  in  Figure  4  where  the  field 
backscattered  from  a  hollow  finite  cylinder  open  at  both  ends  is  of 
interest,  an  additional  interaction  becomes  important  and  it  must  be 
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included.  The  latter  interaction  corresponds  to  waves  coupled  from  the 
interior  waveguide  region  to  the  exterior  region  via  diffraction  at  the 
end  z=-L,  and  vice  versa.  The  fields  of  these  interactions  can  also  he 
found  via  the  UTD  in  a  manner  similar  to  that  done  for  [S^];  this 
field  contribution  to  the  backscattered  field  resulting  from  the 
internal -external  coupling  must  be  added  to  the  result  in  Eq.  (26). 

It  should  be  observed  that  the  latter  is  of  little  interest  in  terms  of 
jet  intake  configuration  of  Figures  1  and  2. 

It  now  remains  to  find  explicit  expressions  for  [Sn],  [Si?], 
[S21L  [$22]*  and  Csb3  t0  complete  the  calculation  in  Eq.  (26).  These 
scattering  matrices  are  formed  via  the  UTD  ray  technique  as  discussed 


r 


SECTION  3 

ANALYTICAL  DEVELOPMENT  OF  THE  ELEMENTS 
OF  THE  SCATTERING  MATRICES  IN  MSM 

The  scattering  matrices  [Sn],  [Sj^L  [S21L  [S22L  and  [St>] 
associated  with  the  canonical  scattering  events  shown  in  Figures  6-8 
were  defined  in  the  previous  section.  In  the  present  section,  the 
elements  of  these  scattering  matrices  are  found  via  the  uniform 
geometrical  theory  of  diffraction  (or  UTD)  ray  technique  [1]  and  its 
modi ficationb  at  caustics.  The  present  ray  analysis  is  restricted  so  as 
to  be  valid  only  in  the  range  0<ir/2,  where  6  is  shown  in  Figures  1-4. 

The  reason  for  this  restriction  is  that  in  the  present  work,  the  RCS  is 
of  interest  only  for  3<n/2.  This  analysis  can  of  course  be  extended  to 
0>ir/2;  however,  that  extension  is  not  reported  in  this  work.  The  UTD 
ray  analysis  is  appealing  because  it  leads  to  a  simple  localized 
physical  description  for  the  scattering  process  in  terms  of  rays,  and  it 
also  provides  relatively  simple  and  accurate  expressions  for  the 
scattered  fields,  and  hence  for  the  scattering  matrices.  In  order  to 
exploit  ray  methods  systematical ly,  the  modal  fields  within  the  circular 


30 


intake  (or  waveguide)  region  have  been  converted  into  a  set  of 
equivalent  rays,  and  vice  versa.  These  scattering  matrix  calculations 
based  on  the  UTD  and  its  modifications  at  caustics  are  found  to  agree 
extremely  well  with  those  obtained  from  the  exact  Wiener-Hopf  solutions 
available  in  the  open  literature  [7,9,10,11,1?]  for  the  problem  of 
diffraction  of  waves  by  a  hollow  semi -infinite  circular  pipe.  This 
agreement  holds  even  for  the  dominant  mode  regime  within  the  circular 
pipe.  It  is  noted  that  while  the  exact  Wiener-Hopf  solutions  for  the 
hollow,  semi-infinite  circular  pipe  provide  very  valuable  checks  on  thr 
UTD  solutions,  they  are  far  more  complicated  to  use  than  the  present  UTD 
solutions.  Furthermore,  the  present  UTD  analysis  can  be  extended  to 
intake  shapes  other  than  those  which  can  be  analyzed  via  the  Wiener-Hopf 
method.  The  reflection  from  the  blades  (i.e.,  blade  scatter)  in  Figure 
8  is  found  by  employing  a  geometric  optical  type  of  approximation  for 
the  blade  current  induced  by  the  modal  fields  which  are  excited  within 
the  intake  by  the  incident  plane  wave.  This  approximation  for  the  blade 
current  is  simple  to  apply,  and  it  is  expected  to  he  reasonably  accurate 
for  a  large  number  of  blades  closely  spaced  together  as  is  generally  the 
case  in  actual  jet  inlets.  It  would  be  worth  investigating  a  more 
rigorous  representation  for  the  blade  currents  in  a  future  study.  The 
UTD  ray  analysis  and  its  modifications  at  caustics  which  are  required  in 
the  development  of  the  elements  of  the  scattering  matrices  is  described 
next. 


The  scattering  matrix  [S]j]  describes  the  fields  scattered  from 
only  the  open  end  (z=0)  of  the  geometries  in  Figures  1-4  when  these 
configurations  are  illuminated  by  a  plane  wave.  In  particular,  [Sj j ] 
describes  the  fields  scattered  from  junction  (I)  corresponding  to  the 
canonical  scattering  event  as  depicted  in  Figure  6  pertaining  to  a 
hollow,  semi-infinite  circular  pipe  illuminated  by  an  EM  plane  wave. 

The  latter  problem  in  Figure  6  can  be  analyzed  efficiently  via  the  UTD 
ray  method  and  its  modifications  which  are  required  along  ray  caustics. 

Consider  the  plane  of  incidence  formed  by  the  incident  ray  and  the 
cylinder  axis  (z-axis)  and  let  the  x-coordinate  lie  in  this  plane.  This 
x-z  plane  of  incidence  intersects  the  rim  at  the  open  end  (z=0)  of  the 
cylinder  at  two  points  denoted  by  Qj  and  Q2  whose  coordinates  are  (x=a, 
y-0,  z=0)  and  (x--a,  y=0,  z=0),  respectively.  Only  Q<n/2  is 
considered  here  for  reasons  mentioned  earlier.  According  to  Keller's 
[21]  generalization  of  Fermat's  principle,  Qi  and  Q?  constitute  the  two 
points  of  first  order  or  single  diffraction  from  the  circular  edge  (rim) 
of  the  open-ended  cylinder.  Additional  higher-order  diffracted  rays 
which  emanate  from  the  curved  rim  after  undergoing  multiple  diffractions 
across  the  aperture  formed  by  that  rim  at  the  open  front  end  (z=0)  also 
exist,  and  their  contribution  must  be  included  in  calculating  the  total 
field  backscatterod  from  the  open  front  end.  In  general,  these  higher 
order  (or  multiply)  diffracted  ray  contributions  beyond  the  second  order 
(or  doubly)  diffracted  ray  interaction  may  be  ignored  for  backscattering 
aspects  away  from  the  z-axis  (or  0=0  direction).  However,  for  on-axis 


(0=0)  backscatteri ng  calculation  the  higher  order  diffracted  ray 
contributions  can  become  important  especially  as  the  radius  (a)  of  the 
cylinder  is  decreased.  It  is  noted  that  the  axial  (or  o=0)  direction 
is  a  caustic  of  the  rays  which  are  singly  and  multiply  diffracted  from 
the  circular  rim  (or  edge)  at  the  open  front,  end  (z=0)  in  the  case  of 
axial  backscatter. 

Initially,  the  UTD  analysis  of  the  problem  in  Figure  6  will  be 
described  for  the  first  order  or  singly  diffracted  fields  which 
contribute  to  the  backscatter.  This  analysis  will  then  be  followed  by 
an  analysis  of  the  multiply  diffracted  ray  fields  which  can  also 
contribute  to  the  backscatter. 


(i)  First  Order  (or  Singly)  Diffracted  Ray  Contributions  to  [Sji] 

The  UTD  ray  field  I115!  which  is  singly  diffracted  from  points  Q\ 
and  Q2  on  the  circular  rim  (edge)  at  the  open  front  end  (z=0)  is  given 
by  [1,21] 


_bsl  bsl  *  bsl  * 
E  =E<l>o  <t>  +  Ee0  0 


(27) 


where 


bsl  2  -j  -j2kasino 

E  4>o  =  ±  I  E 4.  DsU  » *,  60)  e 
0o  4=1  0  h  1  *■ 


1  1  a _  P<t(j  )  e 

2sin9 


-  j  kR 


(28) 


0  not  close  to  zero 


and 
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sec  F(2ktcos2(lil!i)l 


for  a  half  plane.  In  Eq.  (29),  the  “L“  is  an  appropriate  distance 

parameter  within  the  “F"  function  which  involves  a  Fresnel  integral; 

"L"  and  “F"  are  as  defined  in  [1].  When  2kLcos^(  *■  'l;t)>3, 

- 2 - 

F[2kLcos^( ’H  f  3al ;  and  in  the  following,  "F"  is  set  equal  to 
2 

unity  whenever  the  latter  condition  is  true. 

,  ,  .soft . 

In  Eq.  (28),  fls  denotes  {hard}  edge  diffraction  coefficient 
h 

[1,20]  at  (i.e.,  at  Qi  and  Q2  for  t=l  and  1=2,  respectively).  The 

I 

parameters  ,  and  30  at  for  the  backscatter  case  are  defined  as 

♦t  3  +’*  >  *0  =  V2  .  (30) 


where  or  ^2  35  in  figure  10.  The  angle  0^  is  defined  as 

-0,  if  t=l  (i.e.,  for  ray  diffracted  at  Qi) 

=1  I 

0,  if  t=2  (i.e.,  for  ray  diffracted  at  O2) 

The  factor  P ^ ( j )  is  defined  by 

.  1,  when  an  edge  diffracted  ray  from  Qa(ji=1,  or  2)  does 


not  cross  the  axial  caustic  (z-axis) 

P*(j)  =  (32) 

,j,  when  an  edge  diffracted  ray  from  Qt(t=l,  or  2)  crosses 
V  the  axial  caustic  (z-axis) 

The  factor  Py(j)  is  required  since  the  cylinder  axis  is  a  caustic  of  the 

rays  diffracted  from  the  circular  edge  (at  z-0).  Upon  crossing  such  a 

caustic,  the  phase  associated  with  that  diffracted  ray  undergoes  a  jump 

of  +  n/2  radians  which  is  properly  taken  into  account  via  the  factor  j 

appearing  in  the  RMS  of  Eq.  (32).  The  quantity  R  in  Eq.  (28)  denotes 

the  usual  distance  from  0  to  the  far  zone  observation  point  in  the  R 

direction.  From  Eq.  (28),  Eqs.  (4)-(5),  and  the  explicit  form  of  Ds  in 
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Figure  10.  Geometry  of  rays  diffracted  from  points  Qi  and  Q? 
the  edge. 


Figure  11.  Aperture  and  rim  coordinates  at  z=0. 


[1,21],  it  readily  follows  that,  for  a  unit  amplitude  plane  wave 
incident  on  junction  (I), 
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Near  the  caustic  direction  (O+O),  Eq.  (28),  and  hence  Eqs.  ( 33 ) - ( 34 ) 

break  down.  A  uniform  asymptotic  approximation  for  S90  and  which  is 

valid  near  and  at  0=0  may  be  obtained  via  the  method  of  "equivalent" 

edge  currents  which  indirectly  make  use  of  Ds  as  indicated  in  [22]. 

h 

The  equivalent  electric  and  magnetic  edge  currents  Ij  and  Mj, 
respectively  are  positioned  on  the  circular  edge  or  the  rim  at  the  open 
end  (z=0)  of  the  semi -inf inite  cylinder,  and  they  are  given  by  [22] 

I  (<*')  =  -y  ft*  ns('k«4'1  ;80)  (fi.<j>')  _  .  (35) 

1  0  /iT  sin30  nm 


and 


Mjf*') 


.7  fWH  Ph  ( ;gp)  (h1^1) 
"°/jF  sin30  rim 


(36) 


where  4>*  is  the  unit  tangent  vector  to  the  circular  edge  or  rim  (at  z=0) 

in  the  usual  cylindrical  coordinates  at  $>',  and  ,  ip,  6',  8  and 

0  0 

(Tj  )  .  ♦  change  as  a  function  of  position  on  the  rim.  The  coordinates 
rim 

of  any  point  on  the  rim  are  (p'=a,  ,  z'=0);  the  primes  denote  the 


source  coordinates  associated  with  the  equivalent  sources  Ii  and  Mj. 


The  quantities  ^  and  <J/'  in  [)s  are  interpreted  as  the  usual  projected 

h 

incident  and  diffracted  ray  angles  which  are  obtained  by  projecting  the 
incident  and  diffracted  rays  on  the  rim,  and  6'0  is  the  acute  angle 
between  the  incident  ray  and  the  edge  tangent  <j>'.  The  incident  magnetic 
field  TT1  is  related  to  F*  via  TT’=-Y0RxF'  in  which  Z0  is  the  free  space 
impedance,  and  Yo=(^o)”^*  When  the  above  equivalent  currents  are 
incorporated  into  the  usual  radiation  integrals,  one  obtains  the  fields 
radiated  by  these  currents  as 


(37) 


where  k  denotes  the  free  space  wave  number  and  p'=a  p1.  The  quantities 

A  A  A 

' ,  p '  and  R  are  shown  in  Figure  11. 

For  ka  sin  e»l,  Eq.  (37)  can  be  evaluated  asymptotically  via  the 

method  of  stationary  phase.  This  procedure  yields  the  stationary  points 

to  be  at  Qi  and  Q2  on  the  rim  (see  Figure  10),  and  the  corresponding 

stationary  phase  solution  agrees  exactly  with  Eg.  (28),  as  it  should. 

While  the  integral  in  Eq.  (37)  can  be  evaluated  exactly  in  closed  form 

only  for  0=0,  it  can  be  evaluated  approximately  but  accurately  in 

closed  form  for  0  small  as  follows. 

As  a  first  step  in  the  approximate  evaluation  of  the  integral  in 

Eq.  (37)  for  0  small,  one  may  fix  the  0S  in  and  0^  in  to  their 

values  close  to  the  caustic  direction  (0=0)  because,  the  diffraction 

coefficients  0S  are  slowly  varying  with  respect  to  0.  Thus,  for  0  very 
h 

small,  the  integral  in  Eq.  (37)  is  given  approximately  by 
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_bsl  I —  -jkR  2*  .  A  „  ,  .  j2ka$in0cos$' 

E  =  /®1  a  L__}n  (9-0)/d^' [RxRxij)1  )(E1.cos<ti'-E1cos9sin<t>' )]e 

4 n  J  jk  R  $  o  <P  9 

2  n  a  ^  j2kasin9cos<{>' 

-  D.  ( 9=0 )  /d 4> '  [  (Rx 4> * )  (E ^cos ©si n <}> '  +En_cos <i> 1 )  ]e  •  (38) 


The  above  integral  can  now  be  evaluated  in  closed  form  to  yield 


bsl  .  _  ,  .  -jkR 

Efln  =  e  a){Dh(fl=0)J  (2kasinO)-Oc(9=O)cos?0  * 

oo  0  hv  i'  s  ?kasin8  J 


bsl  JT  -jkR 

E  =  eU^  e  a )  {D  ( 0-O)J '  (2kasinO)+D.  ( 9=Q)cos20  Jl(^asin9)|  e 


?kasin9 


It  is  clear  that  the  above  result  which  is  valid  at  and  near  the  caustic* 

i.e.,  for  9  at  and  near  zero  can  blend  into  the  two-point  result  of  Eq. 

(28)  which  is  valid  for  9  not  close  to  zero  (i.e.,  far  from  the  caustic 

direction)  if  the  ns(o=0)  in  the  above  result  for  small  0  is  now  allowed  to 
h 

be  a  function  of  9.  Thus,  one  may  heuristical ly  replace  I)s(9=0)  in  Eqs. 

h 

(39)-(40)  by  D$(°)  to  yield 
h 


bsl  ,  „  ,  ,  -jkR 

E g  =  E1  4  l(l-JL-)  j  (2kasin9)-(l+_L_)cos20  e 

00  9  7  c.oso  V  cos9  2kasino  1  H? 


bsl  , 

E.0  T-Cl4  )  d '  f ?kasi n 0l  +  (  1-  j.  kos^a  Jl(?^s1nO)|  .  ± 

*°  ♦  7  coso  V  coso  JkasTno  1 


The  above  result  for  E$o  corresponding  to  the  singly  diffracted 

9o 

contribution  to  the  backscattered  electric  field  is  now  valid  for  0 

bsl 

small  and  large.  Thus,  E*o  in  Eqs.  (41  )-(42)  is  "uniform"  in  the  sense 

00 

that  it  remains  valid  at  and  near  the  caustic  where  the  ordinary  ray 
solution  of  Eq.  (28)  fails,  and  it  reduces  uniformly  to  Eq.  (28)  far 
from  the  caustic  where  the  latter  expression  is  indeed  valid  and 
accurate.  Furthermore,  it  can  be  shown  to  remain  quite  accurate  even 
for  intermediate  values  of  8  by  comparing  it  with  the  exact  Wiener-Hopf 
solution  when  only  the  first  order  terms  are  retained  in  the  asymptotic 
approximation  of  the  latter  solution. 

For  ka  small,  higher  order  (or  multiply)  diffracted  rays  can 
contribute  significantly  to  the  backscatter  from  the  open  front  end. 

The  expressions  for  these  higher  order  diffracted  field  contributions 
are  discussed  next. 

(ii)  Higher  Order  (or  Multiply)  Diffracted  Ray  Contributions  to  [Sn] 

Let  E^sm  denote  the  electric  field  associated  with  the  rays 
backscattered  from  just  the  open  front  end  (at  z=0)  after  undergoing 
multiple  diffractions  across  the  aperture  corresponding  to  the  open 
front  end.  It  is  convenient  to  express  T^s,n  as 

Ibsm  =  Ebsm  "  £bsm  J  .  (43) 

4>o  ^  9o 

In  determining  T?)Sm  for  backscatter  directions  which  are  not  close  to 
the  axial  (9=0)  direction,  only  the  doubly  diffracted  ray  contributions 


need  to  be  included  for  accuracy  when  ka  is  not  sufficiently  large, 
because  the  contribution  from  all  the  other  higher  order  multiply 
diffracted  rays  in  general  becomes  sufficiently  weak  in  this  case  as  to 
be  ignorable.  It  is  noted  that  m=2  in  the  doubly  diffracted  ray 
contribution  to  I^sm,  i.e.,  is  the  doubly  diffracted  ray 

contribution  to  'P>sm.  If  ka  is  sufficiently  large,  then  even  the  doubly 
diffracted  ray  contributions  become  insignificant  in  comparison  to  the 
singly  diffracted  ray  contributions  for  the  off-axis  (or  9  not  close  to 
zero)  backscatter  directions.  It  is  noted  that  the  entire  rim 
contributes  to  the  diffraction  at  axial  backscatter  in  contrast  to  just 
a  few  isolated  points  of  diffraction  on  the  rim  which  contribute  to  the 
backscatter  in  the  off-axis  case,  because  in  the  case  of  on-axis 
illumination  the  cylinder  axis  (0=0)  is  a  caustic  of  the  rays  diffracted 
by  the  rim.  Consequently,  an  equivalent  rim  current  solution  must  be 
obtained  for  o  at  and  near  zero  in  order  to  express  "P’s?  at  and  close  to 
the  axial  caustic  as  done  for  the  first  order  or  singly  diffracted 
contribution  near  the  caustic.  The  equivalent  current  caustic 
solution  for  will  be  shown  to  generally  yield  numerical  values 
which  are  approximately  the  same  as  those  obtained  from  the  off-axis 
solution  for  TThs‘)  (which  is  based  on  a  few  isolated  points  of 
diffraction  on  the  rim)  in  their  common  domain  of  overlap  in  0.  Thus 
one  could  switch  from  the  caustic  solution  to  the  off-axis  solution  for 
in  their  domain  of  overlap.  For  the  special  case  on  axis  (9=0),  a 
separate  solution  for  T^sm  -jS  also  provided  (in  addition  to  the  caustic 
solution  for  mentioned  earlier)  which  includes  not  only  doubly,  but 


also  triply  and  all  the  higher-order  diffracted  rays.  The  latter 
solution  E^s'n  for  0=0  which  is  admittedly  more  accurate  than  the  caustic 
solution  for  at  6=0,  especially  when  lea  is  not  sufficiently  large, 
is  also  more  cumbersome  to  use  in  practical  applications.  Furthermore, 
even  when  ka  is  not  sufficiently  large,  is  usually  within  one  or 

two  dBs  of  the  value  predicted  by  ft>sm#  Consequently,  in  the  MSM  based 
analysis  for  the  problem  configurations  in  Figures  1-4,  only  the  second 
order  solution  will  be  employed  for  both,  0  near  and  far  from  zero 

to  yield  the  higher-order  diffraction  contributions  to  the  total 
backscattered  field 

In  the  off-axis  backscatter  case,  it  is  seen  via  Figure  12  that  in 
general,  a  set  of  four  rays  contribute  to  the  backscattered  field  after 
undergoing  double  diffraction.  The  total  doubly  diffracted 
backscattered  field  contribution  can  be  expressed  as 

bs2  il2r  i 2 1 r  i 34r  i43r 

E  4>0  3  E<j>  +  E  +  E<j,  +  E  .  (44) 

00  0  0  0  0 

where  the  superscript,  i-l-2-r,  for  instance,  denotes  that  the  incident 
field  (-i-)  is  diffracted  from  point  Qi(-l-)  to  0?(-2-),  and  then  to 
the  far  field  (-r-),  and  vice  versa.  It  is  found  that  points  O3  and 

Q4  are  symmetrical ly  located  with  respect  to  the  xz-plane,  and  the 
cylindrical  coordinates  of  (p=a,  z= 0),  and  Q4(p=a,  <&=$4, 

z  =0 )  are  obtained  by  letting 
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i  .  e3  =  d  .  e3 

A  A  A  A 

d  •  6*  “  “i  •  G4 


to  satisfy  the  generalized  Fermat's  principle  at  points  O3  and  Q/t. 

A  A 

Here,  i  and  d  are  the  unit  vectors  of  the  incident  direction  and  the 
diffracted  ray  direction  across  the  aperture.  e3  and  are  the  unit 
edge  vectors  at  O3  and  Q4,  respectively,  i,  d,  03,  and  e4  are  given  by 


-x  sin  e  -  z  cos  0 
x  cos  a  +  y  sin  a 
-x  sin  $3  +  y  cos  $3 

A  A 

-x  sin  <t>4  +  y  cos  <t>4 


tan  a  = 


sin  $4  -  sin  $3 

cos  4>  -  cos  <}> 

4  3 


Thus,  by  substituting  Eqs.  (47 ) - ( 51 )  into  Eqs.  (45)-(46),  one  obtains 


i>3  =  J.  +  tan”^(sin0)  • 


and  <J>4  =  -4>3 


It  can  be  shown  that 


ilZr  i 2 lr  134r  143r 

^  ~  E  A  *  E  A  “Ea 


which  is  to  be  expected  from  reciprocity.  Therefore,  Fq.  (44)  becomes 


bs2 

E*o 

00 


il2r  i 34r 
2  E 4,  +2  F* 


6  0 


(55) 


where 


il2r  i  -j2ka  cos^l  -jkR 

E  =  f  a  ej _ 2  e  .  (56) 

^7  irkasin0  fos2o  R 


il2r  i  -j2ka  sin2i.  -jkR 

F  -  p  a  e  2  e  .  (57) 

0  0  7  ikasTnS  TT 


:;34r»-E;(0^in20^coS20)eJ2kasi"9COS*3- 


i  34 

_£ _ e 


s(pi34+s) 


~jksrW  e"jkR. 

*  P  R 

(58) 


and, 

TT 

-J'T 

D  = _ i? _ [sec(2!.  -  S)  ±  sec  (I.  +  5)]  ;  (60a) 

l  2/2"iik  sin  S'  *  7  *  ^ 

h  o 


sine' 


1 

/l+sinEo 


t  «  tan"1  (— SjolM 
/1-sin4^ 


(60b) 


Figure  13.  The  doubly  diffracted  ray  for  the  equivalent  current 
solution.  Here,  is  not  necessarily  equal  to 


-3/2 

P134  =  ._-a  ;  p34r  =  -a  llgsjnig  (l+sin2e)  •  (60d) 

(Usin2o)3/2  2sin^9 

Near  the  caustic  region  (0+0),  Eqs.  (56 ) - ( 59 )  break  down  and  the 
equivalent  current  formulation  is  employed  to  correct  the  results  as 
follows.  Let  (a,  0)  and  (a,  v'»  0)  be  the  cylindrical  coordinates 

of  the  first  and  second  diffraction  points  on  the  ri in  as  illustrated  in 
Figure  13.  To  satisfy  the  generalized  Fermat's  principle  at  point  (1), 
it  can  be  found  that 

A  A  /S 

d  =  x  cos  6  +  y  sin  6  .  (61) 

where 

0  =  ir  +  -  sin-1  (si^sin^1 )  .  (62) 


Thus,  v'  and  <t> '  are  related  by 

v '  =  tt  +  4> 1  -  2sin"l  (sin3sin<|>* )  .  (63) 

which  is  illustrated  in  Figure  14.  The  doubly  diffracted  field  in  the 
backscatter  direction  can  be  expressed  in  terms  of  the  equivalent 
current  form  as 


£bs2 


„  JkZo 


2n 

[/  dv1  a^xRxI^v1  +  YQRxMj,v‘)e 


jk  p*R 


(64) 


where 


I2(V) 


_y  fsr.  Ds  (  >  gp) 

°/jF  s(nS0  rim 


(65) 
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©0=xsiny-y  cos  y 


(70a) 


y  =  *  +  -  sin"l(sin8sin<j>‘ )  . 


(70b) 


i  12  i2 

Since  v'  .  ^  =0,  one  only  needs  to  evaluate  EgQ  and  Hpo  .  Thus,  one 
obtains  [1] 


Eil2  =  -  E1  0  (~~e  e-jks 

0o  VsMp?sT 


H112  =  -  H1  D.  I  P _ _  e‘jks 

0o  0o  h/^T 


where 


El  =  [El  si n 4> 1  cos  6  -  E1  (sin9  sin6  +  cost'  cos0  cos6)]eJk 
po  ? 


j  k  a  s  i  n  Ocos  4> ' 


(73a) 


Hfl  3  (H0  n 4*'  cos  5  -  H^.  (sin9  sin6  +  cos<J>'  cos©  cosS)  ]e^kasin0COS'*>  • 
°  (73b) 


tan  6  =  sin  0  cos 
cos  0 


(74) 


Also, 

p  =  a  sin280 _  ,  (75a) 

sin  0  cos  sin  Bq 

and 


s  =  2a  sin  80  .  (75b) 


The  Ds  in  Eqs.  (71)-(72)  is  given  by 
h 


n$ 

h 


-U 

-e 


2/2irk  sin 


[sec( )  +  sec  (jfetil)  ] 


(76a) 


with 


=  n  +  tan"1  (sinecos^1)  ,  ^  «  *  ;  (76b) 

COSO  2 

sinSQ  =  A  Tin 2 0  s i n^4> '  .  (76c) 

It  is  noted  that  the  integrations  in  Eq.  (64)  can  be  carried  out 
numerically  only  for  o<30°.  When  0>30°,  a  complication  arises  in  that 
more  than  one  singly  diffracted  ray  can  strike  the  second  diffraction 
point  (?)  as  may  be  readily  seen  from  Figure  14.  In  the  equivalent 
current  solution,  this  proximity  of  the  ray  fields  at  the  point  (?) 
leads  to  a  singular  result.  On  the  other  hand,  0>3O°  is  usually  a 
sufficient  condition  to  meet  the  criterion  ka  sin  0»1,  so  that  the 
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equivalent  current  solution  of  Eq.  (64)  asymptotically  reduces  to  the 
four-point  result  in  Eqs.  ( 56 ) - ( 59 )  which  is  valid  far  from  the  axis 
(at  0=0).  Furthermore,  the  four  point  douhly  diffracted  ray 
contribution  remains  bounded  and  accurate  even  in  the  case  for  which  the 
equivalent  current  solution  becomes  singular.  Thus,  in  practical 
computations,  it  is  recommended  that  the  equivalent  current  solution  for 
Et>s2  which  remains  valid  at  and  near  the  axial  caustic  (at  0=0)  be 
switched  to  the  off-axis  4-point  solution  for  whenever  0>30°. 

Although  the  doubly  diffracted  ray  contribution  to  the 
backscattered  electric  field  I^s?  has  been  evaluated  above  for  both,  the 
off-axis  or  0*0  case  and  also  for  the  paraxial  case  where  0  is  at  or 
close  to  zero  (and  0=0  is  the  caustic  direction),  a  separate 
solution  will  also  be  obtained  below  for  the  general  mth-order  multiply 
diffracted  field  Tl5sm  in  the  special  case  of  on-axis  incidence  ( i . e. , 
axial  backscatter)  as  mentioned  above.  The  reason  for  developing  an 
expression  for  E*ism  for  the  special  case  of  axial  (0=0)  backscatter  is 
that  it  would  enable  one  to  estimate  the  importance  of  the  triple  and 
higher  order  multiple  diffraction  effects  which  are  not  present  in  the 
expression  for  the  doubly  diffracted  field  "E^s?.  F^s'nwill  be  developed 
below  only  for  the  case  of  axial  {0=0)  backscatter  because  it  is 
tractable  in  this  case.  Such  an  mth-order  result  for  rt5sm  which  would 
be  valid  for  the  off-axis  (or  0*0)  case  cannot  be  obtained  in  any 
tractable  fashion.  Furthermore,  multiple  diffraction  effects  beyond  the 
doubly  diffracted  ones  are  in  general  quite  weak  for  the  off-axis  (or 
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0*0)  case  and  hence  there  is,  in  general,  no  need  to  consider  the 
mth-order  solution  (for  m>2)  in  the  latter  case. 

The  general  mth-order  solution  IT’S"!  for  the  9=0  (axial 
backscatter)  case  will  be  constructed  in  terms  of  the  previously 
described  equivalent  rim  current  concept  as  follows.  First  the  singly 
diffracted  contri buti on  "P^sl ( 6=0)  will  be  listed  to  be  followed  by  the 
expressions  which  will  be  developed  for  the  doubly  and  the  triply 
diffracted  contributions  Tbs2(e=o)  and  P)S^(0=0).  A  sequence  will 
then  be  shown  to  emerge  from  ( 6=0) ,  r^s^(0=0),  and  ^S^(0=O) 
which  can  be  directly  generalized  toT^s^(0=O),  r*5S^(0=0),  etc.,  and 
hence  to  E^sm(0=0)  via  the  process  of  induction.  One  begins  with  the 
first-order  (or  singly  diffracted)  contribution  E^sl(9=0)  which  is 
developed  as  follows.  Referring  to  Figure  15,  the  incident  field  is 
given  for  the  case  of  axial  incidence  by 

V  =  -  x  E0  eJkz  .  (77) 


Since  i|/  =  m  in  Eqs.  (35)-(36),  only  the  equivalent  electric  current 
contributes  to  the  axially  backscattered  field;  i.e.,  the  equivalent 
magnetic  current  gives  a  vanishing  contribution.  Thus,  as  in  Eq.  (37), 


Fhsl(e=0) 


=  jk7.0 

"Tn” 


-jkR 

.  ^ 

R 


2w  jkp'«R 

I  e 

o 


(RxRxI  )ad<t>' 


(78) 


wi  t  It 
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R  =  2 


p‘  =  a(x  cos  <!>'  +  y  sin  <j>‘ )  , 


=  -x  sin  <t>'  +  y  cos 


It  follows  that 


**  i  kR  2n  -  *■ 

Psl(0=O)  =  jk?o  •  e  J  (xsindi'  -ycos$>')  I  d$'  .  (80) 

IT  R  o  1 


where 


l:J1  o  (v  =  4»  =  if )  • 


/Bit  e 


j  Sir)  i p‘ 

kZ0 


at  2*0 


Thus,  one  obtains 


E  (6,0)  -x  £„.*•! 


when  0=0 


This  result  chocks  with  that  in  Eqs.  (41 ) - (42 ) ,  as  it  should,  when  the 
latter  is  specialized  to  o=0. 
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m*. _ 


One  can  extend  the  above  approach  to  obtain  the  axially 
backscattered  field  due  to  the  2nd-order  (or  doubly  diffracted)  ray 
interaction  as 


_bs2  -jkR  2"  ^ 

E  (9=0)=  Jkl  •  z  •  e _  J  (X  sin  4>  -  y  cos  4>  )  I _ d ■(> '  • 

4tt  0  R  o  '■ 

(83) 


in  which  only  the  equivalent  electric  current  I?  is  non-zero  (the 
equivalent  magnetic  current  contribution  vanishes  again).  Specifically, 


I 


•  ^  a  ^) 


(84) 


and 


E2  ’  Eo  •  sin  ♦'  •  °s(*'  =  T  ’  *  =  *’  /m 7T?a)  e 

p  =  -a  . 


-2jka 


(85) 


Therefore, 


2E  ,  .  -j  (2ka+JL) 

l0  =  -  ii£  •  sin  4  _ L  •  e  ^ 

2  kZ0  /URa 


(86) 


By  comparing  I2  in  Eq.  (36)  with  1^  in  Eq.  (81),  one  can  readily  obtain 


Ps2  (6=0)  =  u2  •  Psl  (9=0)  • 


(37) 


where  U£  is  given  by 


-j(2ka  -  1) 
,  4 


/  uka 


(88) 


Likewise,  one  may  obtain  the  axially  backscattered  field  r'1s^(0=O)  due 
to  the  Vd-order  (or  triply  diffracted)  ray  interaction  via  the 
equivalent  current  concept  as 


(0=0)  =  •  1  •  ® _  /  (x  sin  'J>'  -  y  cos  4> ' )  I  ~d«t> '  • 

4  tt  0  R  o  5 


-jkR  2/  A 


(89) 


If'  i  h 


/v  1  I -  J  ' 

ill!l  •  n  (*'«»,*  ■  w)  •  fee  7  at  *=0  .  (9°) 

7.  s  7  \  V 


As  shown  in  Figure  16,  it  is  noted  that  the  ray  singly  diffracted  from 
any  point  on  the  rim  (edge)  traverses  across  the  aperture  to  be 
diffracted  again  from  a  diametrical ly  opposite  point  on  the  rim,  thereby 
producing  a  doubly  diffracted  ray.  Likewise,  the  doubly  diffracted  ray 
wnic'n  traverses  the  aperture  can  in  turn  be  diffracted  once  again  from  a 
diametrically  opposite  point  on  the  edge  to  produce  the  triply 
diffracted  ray,  and  so  on.  Now,  the  doubly  diffracted  ray  which 
traverses  the  aperture  (to  produce  the  triply  diffracted  ray)  will  lie 
on  the  reflection  shadow  boundary  of  the  singly  diffracted  ray  (which 
produces  the  doubly  diffracted  ray).  Due  to  the  fact  that  the  doubly 
diffracted  ray  lies  on  the  reflection  shadow  boundary,  the  field 
associated  with  that  doubly  diffracted  ray  in  the  aperture  is  non-ray 
optical.  Thus,  one  cannot  use  the  edge  diffraction  coefficient  which  is 
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valid  only  for  an  incident  ray  optical  field  to  directly  calculate  the 
field  of  the  triply  diffracted  ray.  However,  it  is  possible  in  this 
special  case  to  decompose  the  non-ray  optical  doubly  diffracted  field  in 
the  aperture  into  two  distinct  ray  optical  components.  The  diffraction 
of  each  of  these  two  ray  optical  components  can  then  be  directly 
calculated  using  the  edge  diffraction  coefficient  which  is  valid  only 
for  a  ray  optical  incident  field.  Let  the  two  ray  optical  components  of 
the  non-ray  optical  doubly  diffracted  ray  field  along  the  aperture 
(i.e.,  along  the  reflection  boundary),  be  denoted  by  and  £3^2. 

respectively  when  the  field  is  evaluated  in  the  limit  as  the  observation 
point  moves  to  the  aperture  from  the  shadowed  side  of  the  reflection 
boundary.  Let  the  total  original  non-ray  optical  field  be  denoted  by 
E3.  Ihus, 

£3 1 1  i  m  =  ^3,1  +  .  (91) 

from  shadow  side 

If  the  limit  of  the  total  non-ray  optical  field  value  E3  in  the 
aperture  was  evaluated  from  the  lit  side  of  the  reflection  boundary,  one 
would  obtain 

_  i  -i  -i  _ ( go )  i 

^3 1 1  i in  »  E3.1  -  f3,2  +  •  (92) 

from  lit  side 


where  E3  is  the  geometrical  optics  field  "reflected"  from  the  edge 
illuminated  by  the  singly  diffracted  field  in  the  aperture.  This  term 
must  be  included  in  the  total  field  calculation  since  it  is  non-zero  on 
the  lit  side  of  its  reflection  boundary.  Furthermore,  it  can  be  shown 


03) 


so  that 


-  a  -i 

=  E3  J  =  e3  . 

lim  from  lim  from 

shadow  side  lit  side 


as  expected  since  the  field  in  the  aperture  after  the  second  interaction 
must  be  continuous  across  the  aperture.  The  ray  optical  components 
F3>1  and  t\t2  are  explicitly  given  by 


_i  ^ 

=  E0  0. 
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r  r 

Here  p  is  obtained  by  [1]  to  yield  p=-a.  pi  and  P2  are  the  two 
principal  radi i  of  curvature  of  the  reflected  wavefront  associated  with 

E^°)’ ;  they  are  given  as  follows  Cl]. 


rc  cose1 


2  cos  01  +  1 


where  rc  and  rc  are  the  two  principal  radii  of  curvature  on  the  surface. 
In  this  case,  e’^0,  r.=-a,  r^=»  ,  i^=a,  and  t^a.  Consequently, 


r  r  0 

Pl  =  "a  '*  ^2  =  2a 


One  thus  obtains 


t3  =  U3  E? 


(100) 


(101) 


where  E2  is  as  in  Eq.  (85)  and  U3  is  given  by 


-j2ka  -j" 

u^  -  JL  e  (1  +  e  ^  j  • 

2  /?  /2irka 


Hue  to  the  similarities  in  the  expressions  between  I2  and  I3  (see 
Eqs.  (84)  and  (90)),  one  readily  obtains 


(102) 


_bs3  hs2 

E  (0-0)  =  113  E  (0=0)  . 


(103) 


Combining  Eqs.  (82),  (87),  and  (103),  one  can  obtain  the  on-axis 
backscattered  field  T^(0=O)  up  to  and  including  the  triply  diffracted 
rays  as  follows. 

s  2  bsm  -  -jkR 

t  (9-0)  «  l  E  (N)  =  x  E  ae _  {l+u.,(ltuj)  •  (184) 

o  m=l  o?  R  1  <1  3M 

Therefore,  from  Eqs.  (4)  and  (104), 

S90  |  0=0  =  |  0=0  “  f  ( l+u 2 ( 1 +u 3 )  1  * 

or,  equivalently,  the  on-axis  RCS  from  the  open  end  of  a  semi-inf  ini te, 
perfectly-conducting  circular  pipe  is  given  by 


a 


1  +  U^(l  +  Ug) 


(106) 


The  above  results  for  P)S1(9=0),  El)S2(o=0),  and  s 3 ( o - 0 )  can  be 

extended  to  obtain  the  mth-order  diffracted  ray  contribution  to  the 

axial  backscattered  field  E^s  as 
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where  F,n  is  obtained  by  generalizing  the  procedure  used  in  obtaining 
F3  of  Eq.  (94).  Each  successive  bounce  (or  diffraction)  of  a  ray 
optical  field  component  in  the  aperture  gives  rise  to  a  non-ray  optical 
field  which  in  turn  can  be  decomposed  into  two  ray  optical  components  as 
in  Eq.  (91).  Furthermore,  each  of  the  latter  two  ray  optical  components 
generates  two  additional  ray  optical  field  components  upon  diffraction, 
and  so  on.  Finally,  as  shown  in  Appendix  B,  the  total  on-axis 
backscattered  field  which  includes  all  m  multiple  interactions,  where 
m-*-»,  is  found  to  be 
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and 
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Also,  the  on-axis  RCS  of  the  open  front  end  (z=0)  is  obtained  from 


0  1 6=0  =  *a 


1  +  iu  c 
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(113) 


Numerical  results  for  the  off-axis  RCS  (normalized  to  ita^)  of  just 
the  rim,  or  the  open  front  end  (at  z=0)  of  the  configurations  in  Figures 
1-4,  is  illustrated  in  Figures  17-18  as  a  function  of  o  for  fixed 
values  of  ka.  These  calculations  are  based  on  Egs.  (41)-(42), 


( 55 ) - ( 59 ) , 

and  ( 64 ) - ( 76 ) . 
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Some  of  these  results  are  compared  with  corresponding  results  based  on 
the  exact  Wiener-Hopf  solution  in  Figure  13.  The  agreement  between 
these  two  independent  results  is  seen  to  be  good.  Likewise,  numerical 
results  for  the  on-axis  RCS  are  shown  in  Figure  19,  as  a  function  of 
a/X  (X  is  the  wavelength).  A  result  based  on  the  Wiener-Hopf  solution 
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i-T> 


1  - I  si  ORDER  UTD 


(a) 


Figure  17.  Off-axis  RCS  (normalized  to  ita^)  of  a  semi¬ 
infinite,  hollow,  perfectly  conducting  circular 
with  ka=7.?61.  The  equivalent  current  technique 
is  employed  for  (K30°. 
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Figure  17.  (continued) 
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Mqure  19.  On-axis  backscattered  RCS  (normalized  to  na^)  from  a 
-infinite,  hollow,  perfectly-conducting  circular 
cylinder  with  various  terms  Included  in  UTD 
calculations. 


SR 


(c)  3rd-order 


Figure  19.  (continued) 


for  the  on-axis  illumination  is  also  compared  with  the  present  on-axis 
solution  which  includes  multiply  diffracted  rays  of  all  orders  (see  Eq. 
(113)).  Again  the  agreement  between  the  two  solutons  is  seen  to  be 
good.  It  is  possible  that  the  inclusion  of  slope  diffraction  effects  on 
the  multiply  diffracted  fields  could  improve  the  accuracy  of  the  present 
solution  for  smaller  ka  even  further  in  the  case  of  on-axis 
il lumi nation. 

B.  Ray  Analysis  for  [S12] 

The  elements  of  the  scattering  matrix  [S12]  describe  the  manner  in 
which  the  modal  fields  propagating  within  the  waveguide  region  radiate 
into  the  exterior  free  space  region  from  the  open  front  end  (at  z=0)  of 
the  semi-infinite  circular  pipe  (waveguide),  as  depicted  in  Figure  7. 

The  elements  of  [S12]  are  found  here  via  the  UTD  together  with  its 
modi f ications  which  are  required  within  caustic  regions.  However,  in 
order  to  use  the  UTD  ray  technique  systematically  in  this  case,  it  is 
necessary  to  convert  each  propagating  waveguide  modal  field  into  an 
equivalent  ray  optical  field  at  the  rim  (edge)  associated  with  the  open 
front  end  (at  z=0).  This  procedure  is  in  contrast  to  the  UTD  ray 
technique  employed  by  Narasimhan  [18]  for  analyzing  the  same  problem  as 
the  one  being  considered  in  the  present  UTD  ray  analysis  of  [S12]. 
Narasimhan  [18]  does  not  convert  his  waveguide  modal  field  into  an 
equivalent  ray  optical  field  in  his  UTD  analysis.  Consequently,  the 
UTD,  which  requires  that  the  edge  diffraction  coefficient  be  employed 
only  for  a  "ray-optical"  incident  field,  cannot  be  strictly  applied  to  a 


"modal"  field  incident  at  the  edge  (or  riin).  On  the  other  hand,  it 
appears  that  if  the  propagation  constant  of  a  mode  approaches  the  free 
space  wave  number,  as  is  the  case  for  a  mode  which  is  very  far  from 
cutoff,  then  the  modal  field  behaves  almost  like  a  ray  optical  field  and 
Narasimhan's  [18]  procedure  may  work  reasonably  well  in  that  case. 
However,  Narasimhan 's  [18]  procedure  is  expected  to  become  inaccurate 
for  a  mode  even  moderately  near  cutoff;  whereas,  the  present  UTD 
procedure  which  is  described  below  will  not  suffer  from  the  limitation 
present  in  Narasimhan's  procedure  [18]. 

Table  1  indicates  the  various  electric  and  magnetic  field 
components  for  the  TEnm  and  TMnm  modes  in  a  circular  waveguide.  The 
propagation  constants  of  the  TEnm  and  the  TMnin  modes  are  given  by 

en!T1  =  /|<2  _  rPnml2  and  fyim  =  Jyl  _  (Pnm respectively.  Consider  a 
a  ‘  a 

typical  modal  field  component  whose  p  and  z  variation  is  given  by 


'  5  =  {Pnm  for  TEnm  mode 

nm  pnm  for  TMnm  mode 

< 

y  =  {S’nm  for  TEnm  mode 
l  nm  $nin  for  TMnn1  mode 


(117) 


One  may  decompose  Jn(Snir|p)  in  Eq.  (117)  as  Jn(x) 
thus,  Eq.  (117)  becomes 
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TADLE  1 


PROPERTIES  OF  MOOES  IN  CIRCULAR  WAVEGUIDES 
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TABLE  1 


(CONTINUED) 


TE  modes  TM  modes 


The  normalization  factors  Nnm  and  Nnm  yield  unit  power  in  each  mode. 
Here,  eon  =  2  for  n  =  0,  and  1  for  n  >  0. 
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(1) 

from  the  large  argument  asymptotic  form  of  H(2)(£  o)  and  for  real 

n  nm 

Ynm»  ^  is  easily  verified  that  the  first  term  in  Eq.  (118)  containing 
Hn  represents  a  conical  wavefront  propagating  away  from  the 

guide  center,  (i.e.,  away  from  the  guide  axis)  as  shown  in  Figure  20; 
whereas,  the  second  term  which  contains  ^(£nmp)  a  conical 
wavefront  (also  for  real  ynm)  propagating  toward  the  guide  center  (or 
collapsing  on  the  guide  axis).  Clearly,  the  field  of  the  conical  wave 
which  propagates  away  from  the  guide  axis,  i.e.,  which  propagates  toward 
the  r!m  at  the  open  front  end  (at  z=0),  constitutes  the  ray  optical  field 
which  is  incident  at  the  rim  (edge).  This  incident  field  undergoes 
diffraction  at  the  circular  edge  and  thus  contributes  to  the  field 
radiated  into  the  exterior  region  by  the  modal  field  within  the 
waveguide.  Far  from  the  waveguide  axis  (0=0),  one  may  use  the  UTI)  to 
calculate  the  field  radiated  from  the  open  end  in  terms  of  the  two 
diffracted  rays  shown  in  Figure  21.  In  accordance  with  the  generalized 
Fermat's  principle,  these  two  rays  originate  at  the  two  points  of 
diffraction  on  the  circular  rim  (at  z=0)  which  correspond  to  the 
intersection  of  the  rim  with  the  plane  defined  by  the  far  zone 
observation  point  and  the  axis  of  the  open-ended  circular  waveguide. 

Tne  first  order,  or  the  singly  diffracted  field  associated  with  the  two 
rays  diffracted  from  the  rim  is  directly  produced  by  the  incident  ray 
optical  field  which  is  incident  there.  A  part  of  this  singly  diffracted 
field  also  initiates  multiply  diffracted  rays  which  bounce  back  and 
forth  across  the  aperture;  at  each  bounce,  a  part  of  the  energy  is 
again  radiated  to  the  far  zone  observation  point  along  the 
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aforementioned  same  two  ray  paths  after  being  multiply  diffracted  from 
the  same  two  points  of  diffraction  on  the  rim.  The  latter  contribution 
to  the  radiation  field  which  results  from  the  multiply  diffracted  rays 
may  be  included  for  improved  accuracy  when  ka  is  not  sufficiently  large. 
If  ka  is  sufficiently  large,  then  the  field  of  only  the  singly 
diffracted  rays  remains  significant;  thus,  the  singly  diffracted  fields 
are  the  dominant  contributors  to  the  far  zone  radiated  fields,  when  6 
is  not  close  to  zero.  The  UTD  electric  field  Tr  which  is  associated 
with  these  two  singly  diffracted  rays  is  given  by 
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{Note  that  in  this  case*i; 
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where  the  terms  within  the  square  brackets  in  fcqs.  (12D)-(123) 
UTD  edge  di f tract i or.  coefficients  multiplied  by  a  constant;  in 
particular. 
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(122) 


(123) 


are  the 


(124a) 


(124b) 


(1?5) 


S'  =  cos-1(e'  /k)  for  TE  modes 
nm  nm  '  nm 

S'  =  cos-1 (8/k)  for  TM  modes 
nm  1  nm  '  nm 

It  is  noted  that  the  half  cone  angles  of  the  conical  wavefront 

associated  with  the  rays  incident  on  the  rim  (at  z=0)  as  shown  in 

) 

Figures  20-21  are  denoted  by  Snm  and  5nin  for  the  TEnm  and  the  TMnm 

I 

modes,  respectively.  These  angles  Snm  and  5nm  appear  in 

I 

Eqs.  (120)-(123).  Note  that  0={Snm  or  5nm>  corresponds  to  the 

incident  and  reflected  shadow  boundary  directions  associated  with  the 

incident  conical  wa-e  at  the  edge  as  shown  in  Figure  20.  The  above 

~  „  9T5 

Eqs.  ( 120 ) - ( 123 )  are  presented  for  2kLcos?(  2  )>3  where  L=a/sin£. 

I 

Consequently,  these  equations  need  to  be  modified  if  0+Smn  or  6nm 

by  including  the  "F"  functions  in  Ds  as  indicated  in  Eq.  (29).  Since 

h 

0=0  is  a  caustic  direction  for  the  diffracted  rays,  the  two  point  first 
order  diffraction  solution  for  the  radiated  fields  in  Eqs.  (120)-(123) 
becomes  singular  at  9=0  and  therefore  cannot  be  employed  at  and  near 
that  direction.  The  equivalent  ring  current  concept  used  previously  for 
evaluating  the  fields  near  the  caustic  direction  (0=0)  cannot  be 
employed  in  this  case  because  the  incident  conical  wavefront  which 
undergoes  diffraction  at  the  rim  edge  (at  z=0)  gives  rise  to  a 
non-cylindrical  wavefront  in  the  immediate  vicinity  of  that  edge.  It  is 
noted  that  the  diffracted  wavefront  in  the  immediate  vicinity  of  the 
diffracting  edge  must  be  cylindrical  in  order  for  the  equivalent  ring 
current  concept  to  be  valid.  The  reasons  for  the  existence  of  a 
non-cylindrical  diffracted  wave  near  the  circular  edge  in  this  radiation 
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problem  are  two-fold,  first,  the  total  field  in  the  immediate 
vicinity  of  the  diffracting  edge  consists  not  only  of  the  edge 
diffracted  field  which  may  possess  a  cylindrical  wavefront,  but  also 
consists  of  the  incident  and  reflected  conical  wavefronts  associated 
with  the  incident  and  reflected  rays  which  exit  directly  from  the 
aperture.  The  latter  incident  and  reflected  conical  rays  represent  the 
ray  equivalent  of  the  modal  field  which  impinges  on  the  edge  (at  z=0) 
from  within  the  waveguide  region  as  described  earlier.  Secondly,  the 
wavefront  of  the  edge  diffracted  field  in  the  immediate  vicinity  of  the 
diffracting  edge  may  also  become  non-cyl i ndri cal  if  it  is  produced  by  an 
incident  conical  wave  corresponding  to  a  waveguide  mode  far  from  cutoff 
(i.e.,  when  6mn  and  6nm  are  close  to  zero)  because  in  this  case  an 
observation  point  in  the  close  vicinity  of  the  diffracting  edge  (at  z=0) 
can  lie  within  the  incident  and  reflection  shadow  boundary  transition 
regions  associated  with  the  incident  and  reflected  conical  rays  at  the 
edge.  Therefore,  a  procedure  different  from  the  equivalent  current 
concept  must  be  employed  to  find  the  field  near  6=0,  which  is 
radiated  from  the  open-ended  waveguide. 

A  procedure  is  described  below  for  obtaining  T  which  remains 
valid  for  o  not  only  at  and  near  zero,  but  also  for  0  far  from  zero, 
furthermore,  for  o  far  from  zero,  this  result  reduces  to  the  first 
o'-'dor  (singly  diffracted)  two  point  UTO  diffraction  solution  for  ^ 
obtained  earlier  in  Fqs.  (119)  -(123)  where  the  latter  is  indeed 
expected  to  be  aci. urate.  This  procedure  is  based  on  a  modification  of 
Ufimtsev's  Physical  Theory  of  Diffraction  [24];  it  is  described  next. 


The  radiated  field  E  can  be  expressed  in  terms  of  the 


Stratton-Chu  radiation  integral  [25]  over  the  equivalent  surface 

currents  that  exist  on  the  surface^  which  tightly  encapsulates  the 

semi-infinite  open-ended  circular  waveguide.  Let,«a  denote  the 

circular  disc  shaped  aperture  surface  at  z= 0  and  /&c  denote  the  surface 

_r 

formed  by  the  exterior  walls  of  the  waveguide.  Thus,  E  is  given  by 
[25] 


_r  -jkR  -  -  _  -  _  . j k R  - R 

E  =  Jjfo  e _  J  J  dS'  [RxRxJ  +  Y  RxM  ]e  •  (126) 

TT  R  A  =  £  $  +  ^  c  eq  0  eq 

where  R1  is  a  vector  from  the  origin  (0)  to  any  point  on/ :  .  Since  one 
is  concerned  with  the  use  of  rays  in  this  study,  it  is  convenient  to 
replace  the  equivalent  electric  and  magnetic  current  densities  and 
Meq  by  their  asymptotic  approximations  valid  for  large  ka.  Hence,  to  a 
first  order  of  asymptotic  approximation. 


_  I  .jnodal  d 
^eq  '  /sd  ~  ^eq  +  ^ 


K 


eq  1  #a  "eq 


jnodal  d 
Mpn  +  1] 


'eq  1  j 


(127) 

(128) 

(129) 


It  is  noted  that  Men  |  x  =0  because  is  a  perfectly  conducting 
'  A  c  *  C 

exterior  wall  of  the  waveguide.  In  Eqs.  (127)  and  (128), 


_modal  *  _inc  _modal  i  nc 

^eq  =  z  x  Kmodal  1  ^eq  =  wnodal  x  z  •  (130) 


_inc 


where  Hin0Cja  1  ’s  the  magnetic  field  associated  with  the  waveguide  modes 


which  impinge  upon  the  circular  edge  (at  z=0)  from  within  the  waveguide. 

inc  inc 

Likewise,  Em0dal  is  the  electric  field  associated  with  Hmodal*  0n 

the  other  hand,  the  current  densities  "3^  and  l3!^  in  Eqs.  (127 ) - ( 129 ) 


deq 


odal 


denote  the  corrections  to  Jpn  and  Mpn  that  result  from 

inc  _inc 


jnodal 

^eq 


the  effects  of  diffraction  of  the  incident  modal  fields  C^modal ,  Hmodal ) 
at  the  circular  rim  (edge).  The  and  TT  will  be  described 


subsequently  in  more  detail. 

It  is  obvious  from  the  above  definition  in  Eq.  (130)  that 

_modal 

Jeq  =0  on  js c.  Consequently,  the  result  in  Eq.  (126)  can  be 
decomposed  into 


f  -  tp„ 


(131) 


where 


.ikZ0  e'jkR 

"Tit - R 


/ /dS 
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-  »  _modal 
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_modal 


+  Y  RxM 
o  eq 


jkp' *R 
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(132) 


since  "S'-  p'  on  /8a,  and, 
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S  .v7  -jkR  «  -  _d  -  _d  jkp  *R 

Eu  =  e-R_  J  /  dS '  [RxRxJ  +  YqRxM  [e 

d 


-J'kR  ,  ,  .  .  _d  jkR ' -R 

+  e  I  I  dS'  [RxRxJ  ]e 

4  TT  R  A  * 


(133) 


Obviously,  the  term  Tpg  represents  the  vector  Ki rchhoff-Huygen ‘s 

aperture  integral  for  field  radiated  by  the  open  end  of  the  waveguide  in 

which  one  employs  the  common  approximation  that  the  field  in  the 

aperture  at  the  open  end  consists  of  only  the  waveguide  modal  fields 

which  are  incident  (on  the  aperture)  from  within.  The  diffraction 

current  densities  7  and  in  Eq.  (133)  for  e(j  provide  the  correction 

r  _r 

to  the  incident  mode  current  approximation  in  Tpo.  The  integral  for  Epo 
in  Eq.  (132)  has  been  essentially  evaluated  by  Chu  [2]  and  is  also 
described  by  Silver  [24].  The  result  can  be  rewritten  in  a  more 
convenient  form  for  the  present  analysis  as 


/  =  (eS, 

PO 


POe 


e  +  E 


P04> 


♦) 
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where 
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!t  is  further  convenient  to  identify  the  terms  within  the  square 
brackets  of  Fqs.  (134)-(137)  as  follows. 


$  ,  for  TE  modes 

5  _  j  run  nm 

6  »  for  TM  modes 

nm  nm 


nP0  nP0 

where  us  and  are  the  soft  and  hard  physical  optics  (PO)  edge 

diffraction  coefficients.  Far  from  the  axis  (9=0);  i.e.,  for  large  ka 

•  r 

sin  9,  and  9  not  close  to  (6nm  or  Sr,m),  the  above  results  for  Ypo 

in  Eqs.  (134 ) - ( 137 )  reduce  to  the  two  point  diffraction  form  as  in 
,  PO 

Eqs.  ( 120 )- ( 122 ) ,  but  with  Ds  appearing  in  place  of  Ds  as  one  might 

h  h 

_r 

have  anticipated.  It  will  be  shown  below  that  only  when  the  Eu 

contribution  of  Eq.  (133)  is  added  to  the  Tpo  of  Eqs.  (134)- ( 137 )  to 

obtain  "E  as  in  Eq.  (126),  and  the  result  for  ^  is  subsequently 

asymptotically  approximated  for  large  ka  sin  9  and  9  not  close  to 
1 

{Snm  or  5nm},  then  the  resulting  asymptotic  form  agrees  with  the  two 

point  UTD  diffraction  solutions  in  Eqs.  (120 )- (122 )  which  is  in  terms  of 

Ds.  Returning  next  to  the  evaluation  of  Ty  in  Eq.  (133),  it  is  noted 
h  d  d 

that  J  and  M  in  that  equation  represent  the  current  densities 
associated  with  the  edge  diffracted  field  components.  The  surface 
integrals  over  &  a  and  &  c  must  be  evaluated  asymptotical ly  along 
the  coordinate  transverse  to  the  circular  rim  thereby  reducing  those 
surface  Integrals  to  a  line  Integral  over  the  circular  rim  (edge)  at 
z=0.  A  stationary  phase  evaluation  of  that  line  integral  (or  the  rim 
Integral)  for  ka  sin  9  large  yields  a  closed  form  result  for  r[j  . 

However  that  closed  form  result  for  large  ka  sin  0  is  not  valid  at  and 
near  the  caustic  direction  (9=0).  Clearly,  the  rim  integral  therefore 
cannot  be  evaluated  via  stationary  phase  at  and  near  the  caustic; 
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instead,  it  must  be  evaluated  numerically  for  9  at  and  near  zero. 

While  the  reduction  of  the  surface  integral  for  Ty  to  the  line  or  rim 

integral  is  conceptually  straightforward.  It  is  quite  tedious  since 
d 

and  "ft  must  be  expressed  in  the  proper  coordinate  systems  to  carry  out 
the  asymptotic  end  point  evaluation  of  the  Integral  along  the  coordinate 
transverse  to  the  rim.  Of  course,  and  are  known  from  the  UTD 
edge  diffracted  fields  in  the  problem  of  plane  wave  diffraction  by  a 
half  plane  [1].  These  fields  are  also  identical  to  those  obtained  by 
Sommerfeld  in  his  solution  to  the  problem  of  plane  wave  diffraction  by  a 
half  plane  [26].  In  particular. 


_d 
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*  _d 
z  x  H 


on  ^  a 
on  ^c 


Md 
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X  z 


on 


A 


(140) 


in  which  f  and  H  are  the  edge  diffracted  electric  and  magnetic 

fields  produced  by  the  diffraction  of  the  incident  conical  wave  (which 

is  locally  plane  for  large  ka)  at  the  edge  (see  Figure  21).  A  far  more 

convenient  alternative  to  the  Uflmtsev  based  formulation  for  t  in 

Eq.  (131)  which  is  given  In  terms  of  the  surface  integral  representation 

for  'Ey  in  Eq.  (133)  is  to  use  the  concept  of  equivalent  Ufimtsev  edge 

r 

currents  to  obtain  'Ey  directly  In  terms  of  a  rim  Integral  rather  than 
the  surface  integral.  This  modification  of  the  Uflmtsev  based 
formulation  in  Eq.  (133)  is  described  below. 


86 


Consider  the  problem  of  the  plane  wave  scattering  (and  diffraction) 

by  a  half  plane  in  Figure  22  for  which  the  exact  solution  is  well  known. 

It  is  not  difficult  to  see  in  this  case  that  when  the  incident, 

reflected,  and  diffracted  fields  present  in  that  exact  solution  are 

employed  to  calculate  the  equivalent  electric  and  magnetic  current 

densities  in  the  Stratton-Chu  Integral  [2,25]  for  the  electric  field 

4  radiated  (or  scattered)  by  the  half  plane,  then  one  should  again 

recover  the  original  exact  solution  for  the  fields  scattered  (and 

diffracted)  by  that  half  plane.  The  Stratton-Chu  formulation  [2,25]  for 
r 

Thp  is  the  same  as  in  Eq.  (126),  and  the  contours  >8a  and  ,&c  for  the 
present  half  plane  case  are  as  shown  in  Figure  22.  Note  that/ta  of 
Figure  22  locally  simulates  the  aperture/8a  of  the  circular  waveguide 
(see  Eq.  (126)).  Again,  the  field  "dip  can  be  decomposed  into 


EpO(hp) 


(141) 


as  in  Eq.  (131).  The  equivalent  electric  and  magnetic  current  densities 

_  _  _r 

Jeq  and  Meq  that  appear  in  the  Stratton-Chu  integrals  for  Epo(hp)  and 

^U(hp)  are  given  by 
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A 
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CD  INCIDENT 


figure  22.  Diffraction  mechanism  from  a  perfectly-conducting 
half  plane. 
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where  I”1 ,  F*",  and  are  the  Incident,  reflected  and  diffracted 
components,  respectively,  of  the  electric  field  present  in  the  exact 
solution  to  half  plane  diffraction  problem  of  Figure  22  and  likewise  the 
TT^ ,  TTr ,  and  are  the  magnetic  field  components  associated  with  F  ,  FP, 
and  F**.  It  can  be  shown  that  an  evaluation  of  the  Fpo(hp)  integral  in 
terms  of  7eq  and  Feq  of  Eqs.  (142)-(144)  yields 


Fpo(hp)  =  | -it)  +  F^d^'  | -ti)  +  Epo 


(145) 


where  <|7,  =  0O, 


1,00 
1/2  ,  5  -  0 
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(146) 


twi  -r' 


=  PO  -jkR 
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(147) 


in  which  6^®(ir±0  ®0;7)  is  the  uniform  physical  optics  edge  diffraction 

coefficient  in  dyadic  form  and  it  corresponds  to  the  receiver  location 
at  x  >  0,  This  dyadic  fi**®  can  be  expressed  in  terms  of  the 

~P0 

scalar  diffraction  coefficient  Os  in  a  manner  similar  to  that  for  D  as 

~P0  h 

indicated  in  [lj.  The  0$  is  explicitly  given  by 

h 


39 
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~P0  -  j 11 

Os  (*,*';!)=  -g  1  |tan(»-*l)Ff2kttcos2(.^V)]Ttan(^)Ff2kRcos2^»^^»,)]} 
h  ‘  2/2ik  2  2  2  2 

(148) 

When  2k3  cosf>(^li_j>3,  the  "F"  functions  in  Eq.  (148)  are 
^  ~PO  PO 

essentially  unity,  then  Ds  (ip  , xp;  I.)  reduces  to  Ds  ( 4»* ;  1)  defined 

h  2  h  2 

earlier  in  Eq.  (139);  i.e.. 
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~P0  PO 
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(150) 


if  2k R  cos^AlA  >3. 


Also,  an  evaluation  of  E-u(hp)  in  terms  of  the  7eq  and  fteq  of 
Eqs.  ( 142 ) - ( 144 )  must  yield 


^U(hp)  =  . 


051) 


where 


_jj  i  *u  -jkR 

Eu  =  T  •  H  (*±9, 00;ir)  .  £ -  .  (152) 

at  0  7  /W 

=U 

The  0  may  be  referred  to  as  a  "Ufimtsev"  half  plane  edge  diffraction 

=U 


coefficient  in  dyadic  form.  This  dyadic  D  can  also  be  expressed  in 

ft  -u 

terms  of  the  {^grd}  scalar  diffraction  coefficient  Os  ,  where 


U 

Ds(  ^  ;*)  = 
h  7 


-i* 

e  JT 
2/T^k 


sec(jMll)-tan(JM[l^  +  ^sec(y  j-tanfj'+i'1 


(153) 


PO  U 

It  is  noted  that  in  contrast  to  Ds  ,  the  Ds  is  valid  at  and  near  the 

h  h 

incident  and  reflection  shadow  boundaries  in  the  problem  of  Figure  22; 

these  shadow  boundaries  occur  at  |  |  =  n.  Consequently,  the 

_d  _d 

wavefront  associated  with  Ey  and  Hy  is  always  cylindrical  and  one  can 
always  view  the  field  ^(hp)  as  being  produced  by  an  equivalent  edge 
(line)  current  referred  to  as  an  equivalent  Ufimtsev  edge  current;  in 
particular,  one  may  write 


e-M 

w 


_1 
if  E 


“  i 

=  y  E  only; 


and 


e-jkR 

~7T 


(154) 


(155) 


if  H  =  y  H  only.  Thus,  the  equivalent  electric  and  magnetic  Ufimtsev 
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half  plane  edge  currents  Iy  and  Mg,  respectively  are  given  by 


(156) 


(157) 


The  above  Itj  and  Mg  are  obtained  by  a  direct  comparison  of  Eq.  (152), 

with  Eqs.  (154)-(155),  respectively.  These  currents  Iy  and  My  which 

are  valid  for  normal  incidence  on  the  half  plane  can  also  be  directly 

generali?ed  to  oblique  angles  of  incidence  and  diffraction  on  the  edge 

[23],  Next,  comparing  the  half-plane  geometry  in  Figure  22  with  the 

local  geometry  at  any  point  on  the  circular  edge  of  the  open-ended 

circular  waveguide  as  shown  in  Figure  21,  it  is  not  difficult  to  see 

that  one  can  define  a  e0  ={5n'm  or  6nm}  at  any  point  on  the  circular 

edge  (at  z=0)  so  that  the  field  Iy  in  Eq.  (133)  can  now  be  viewed  as 

being  produced  by  the  equivalent  electric  and  magnetic  Ufimtsev  half 

plane  edge  currents  Iy  and  My.  Thus,  using  the  equivalent  current 

concept,  one  may  directly  employ  the  Stratton-Chu  radiation  integral  for 
_^r 

Ey  in  terms  of  these  rim  currents  Iy  and  My,  to  obtain 
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(158) 
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with  £  *  d£ '  =  ^'adV  for  a  circular  rim,  and  'p,  =  p,a.  Eq.  (158)  is  in 

contrast  to  the  surface  integral  in  Eq.  (133)  which  can  be  reduced  to  a 

riin  integral  only  after  asymptoti cal ly  evaluating  the  integral  along  the 

coordinate  orthogonal  to  the  rim.  Thus,  according  to  the  proposed 

_r 

modification  of  the  Ufimtsev  based  formulation  for  obtaining  E  ,  one 
may  write 


_r 
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ik7  r 

+  JkZo  e _  £ 
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(159) 


where  Ipo  is  still  the  same  as  in  Eq.  (132).  Only  the  Eg  term  in 
Eq.  (133)  which  occurs  in  Eq.  (131)  is  now  replaced  by  the  line  integral 
around  the  edge  as  in  Eq.  (158)  above.  Specifically,  the  equivalent 
Ufimtsev  edge  currents  and  My  in  Eqs.  ( 156 ) - ( 157 )  which  are 
modified  so  as  to  make  them  valid  for  oblique  angles  of  incidence  and 
diffraction  on  the  edge,  without  the  restriction  of  generalized  Fermat's 
principle  for  edges  [23]  are  now  defined  by 


I 


U 


and 


1 

sin0o 


M 


U 


Z  1 

°/jF  si  n  60 


(160) 


o*n 


93 


As  shown  in  Figure  23,  ;<0  and  30  are  the  acute  angles  made  by  the 
incident  and  diffracted  rays  with  the  tangent  (l‘)  to  the  edge  at  the 
point  of  edge  diffraction,  respectively.  In  the  present  case  of  conical 
wave  incidence  at  the  edge,  the  axis  (or  caustic)  of  the  conical  wave 

I  t 

coincides  with  the  waveguide  axis  and  hence  30  =  v/2  or  sin  30  -1  in 
Eqs.  (160 )- ( 161 ) ;  on  the  other  hand,  30  is  given  by 

si nb0  *Jl  -  sin^O  sin^1  .  (162) 

The  above  Iy  and  My  in  F.qs.  (160)  — (161 )  reduce  to  those  given 
previously  in  Eqs.  ( 1 56 ) - ( 157 )  for  normal  incidence  on  the  edge  in  which 

T 

30  =  n/2  and  li0  =  n/2. 

It  is  important  to  note  that  the  result  in  £q.  (169)  involving  the 
equivalent  Ufimtsev  edge  current  represents  a  useful  modification  of  the 
original  Ufimtsev  formulation  in  Eq.  (131 )- ( 133 ) .  Furthermore,  this 
modified  result  in  Eq.  (159)  is  quite  general  in  that  it  can  be  directly 
extended  to  other  antenna  and  scattering  problem  conf i gurations 
involving  edges.  Basically,  the  modified  result  in  Eq.  (159)  implies 
that  the  radiation  or  scattering  can  be  calculated  according  tc  a 
physical  optics  approximation  and  be  corrected  by  the  equivalent 
iJfimtsev  edge  current  integral  which  is  then  added  to  the  physical 
optics  contribution  of  Eq.  (159).  The  result  in  Eq.  (159)  is  to  he 

used  for  aii  0  (both  near  and  far  from  zero)  and  it  is  expected  to  be 

_r 

accurate  as  long  as  Epy  is  not  vanishingly  small.  For  0  -*•  n/2,  a 
doubly  diffracted  ray  becomes  important;  its  contribution  can  be 
calculated  separately  via  the  UTO  and  it  is  simply  added  to  the  solution 
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Figure  23.  diffraction  from  a  curved  edge 


for  /  in  Eq.  (159).  For  large  ka  sin  e,  the  physical  optics 

contribution  in  Eq.  (159)  together  with  Eqs.  (134 ) - ( 138 )  can  be 

approximated  such  that  the  Bessel  functions  involving  ka  sin  0  and  the 

Hankel  functions  therein  are  replaced  by  their  large  argument  asymptotic 

forms.  When  this  asymptotic  approximation  is  made,  the  physical  optics 

result  has  the  same  form  as  the  two  point  UTD  diffraction  solution  of 

Eqs.  ( 120 ) — ( 122 )  except  that  it  now  contains  the  physical  optics  edge 

J>0 

diffraction  coefficients  us  Instead  of  the  UTD  edge  diffraction 

h 

coefficients  Ds  of  Eqs.  ( 120 )-( 122 )  as  mentioned  earlier.  Furthermore, 
h 

the  rim  integral  in  Eq.  (159)  can  also  be  asymptotically  evaluated  for 

large  ka  sin  e  via  the  method  of  stationary  phase  to  yield  a  result 

which  is  again  like  the  two  point  UTD  diffraction  solution  except  that 

U 

it  contains  the  Ufimtsev  ede  diffraction  coefficients  Ds  instead  of 

h 


Ds  .  However,  the  sum  of  the  asymptotic  approximations,  for  large 
h 

ka  sin  Q,  of  the  physical  optics  and  the  equivalent  Ufimtsev  rim  current 
integrals  agrees  with  the  two  point  UTD  diffraction  solution  in 
Eqs.  ( 120 ) - ( 1 22 )  because 


P°  U  , 

Ds  Ds  +  Ds  .  (163) 

h  h  h 

It  may  be  remarked  that  while  the  two  point  UTD  diffraction  solution  in 

Eqs.  { 120 ) - ( 122 )  is  accurate  for  ka  sin  0  large,  the  solution  in 
Eq.  (159)  is  valid  even  when  ka  sin  9  is  not  large;  furthermore,  it  is 
expected  to  be  more  accurate  than  the  UTD  solution  for  o  ♦  Snm  or 

I 

6nin  even  if  the  "F"  functions  present  in  Ds  of  Eq.  (29)  are  included 

h 

in  the  latter.  One  notes  that  the  Bessel  functions  Jn  (ka  sin  6)  or 
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Jn(ka  sin  8)  In  the  physical  optics  solutions  of  Eqs.  (135)-(1 38 )  for 

» 

the  TMnm  and  TEnm  cases  vanish  at  0  =  «nm  and  0  =  6nm,  respectively; 

PO 

thereby  keeping  the  solutions  bounded  even  when  Ds  ♦  »  at  8  +  Snm 

»  h  PO 

or  6nin.  Clearly  the  Bessel  functions  which  keep  Ds  bounded  are 

h 

to  be  viewed  as  Incident  and  reflection  shadow  boundary  transition 

functions  which  yield  a  result  for  Tpg  that  is  uniformly 

valid  across  the  incident  and  reflection  shadow  boundaries  associated 
1 

with  0  =  {5nm  or  6nm}.  It  appears  that  the  Bessel -type  transition 
functions  are  more  accurate  than  the  Fresnel  integral  ("F"-type) 
transition  functions  of  the  UTD  at  least  for  the  present  waveguide 
problem  and  especially  for  ka  sin  8  small.  Finally,  it  may  be  remarked 
that  the  two  point  UTD  diffraction  solution  valid  for  large  ka  sin  e  can 
in  some  problems  be  generalized  heuristical ly  to  be  valid  at  and  near 
the  caustic  at  8=0  by  recognizing  the  two-point  solution  to  be  an 
asymptotic  approximation  of  some  functional  form  which  yields  a  bounded 
result  at  the  caustic.  However,  such  a  heuristic  procedure  for 
generalizing  the  two-point  UTD  diffraction  solution  to  the  caustic 
region  is  risky  in  the  present  open-ended  circular  waveguide  case  and  it 
can  lead  to  an  erroneous  result.  In  conclusion  therefore,  the  result 
for  in  Eq.  (159)  is  to  be  employed  for  all  0. 

A  direct  application  of  either  the  UTD  result  of  Eq.  (119)  (for  o 
not  close  to  zero)  or  even  the  use  of  the  expression  for  /  in  Eq. 

(159)  will  generally  fall  to  predict  the  correct  value  of  if  any 
component  of  Tpo  vanishes  at  all  aspects  (see  Eq.  (134)  for  Ipo). 

Such  a  situation  arises,  for  example,  in  the  prediction  of  the  <fr 
component  of  £po  f°r  the  radiation  by  an  open-ended  circular  waveguide 
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r 

excited  In  the  TMnm  mode.  In  the  latter  case,  the  term  Epo <(>  is 

identically  zero  In  Eq.  (134).  From  a  practical  point  of  view,  even 
r 

though  Epo^  =0,  its  overall  effect  in  the  RCS  calculation  is  small  in 
comparison  to  the  other  scattering  effects  which  dominate;  hence  an 
unavai 1 abi 1 ity  at  the  present  time  of  an  accurate  estimate 
based  on  ray  methods  and  its  modifications  for  T** •  $  radiated  from  the 
open  end  of  a  circular  waveguide  excited  by  a  TMnm  mode,  is  not  a 
serious  problem  in  the  present  RCS  study.  Some  typical  results  for  TTr 
based  on  Eq.  (159)  are  shown  In  Figures  24-25  where  they  are  seen  to 
compare  well  with  those  based  on  the  exact  Wiener-Hopf  solution  [12], 

A  A  r 

provided  the  0  and  <J>  components  of  Epg  are  non-zero. 

C.  Ray  Analysis  for  [Spi] 

The  scattering  matrix  [S2i]  describes  the  transformation  or 

coupling  of  the  incident  plane  wave  field  into  the  waveguide  modes  as 

illustrated  in  Figure  6.  This  scattering  matrix  [Spi]  was  defined 

earlier  in  Eq.  (11).  The  problem  of  determining  [S21]  is  the  reciprocal 

of  the  problem  of  determining  [S12]  (see  Figures  6-7).  Thus,  a 

knowledge  of  [S12],  or  more  specifically  S0nm,  $$nm,  S0nm,  and  S^, 

and  the  use  of  the  reciprocity  theorem  for  electromagnetic  fields  allows 

one  to  calculate  the  elements  sinm,  S^,  and  of  [S21], 

and  to  hence  calculate  [S2 1 ]  in  terms  of  Thus,  it  can  be  shown 

that  the  elements  of  [S21]  are  simply  related  to  the  elements  of 
_  T 

[S12J  *  *1ere  T  denotes  the  transpose  matrix  operator.  The  reciprocity 
theorem  is  applied  as  follows. 
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Figure  24.  (continued) 
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First  consider  the  geometric  plane  defined  by  z=-L  inside  the 
semi -inf inita  pipe,  and  let  S0  denote  the  area  of  this  waveguide  cross 
section  at  z=-L  as  shown  in  Figure  26.  Then  let  S  denote  the  surface 
area  which  tightly  encapsulates  the  complete  outer  (exterior)  surface  of 
the  semi-infinite  pipe,  and  also  a  portion  of  the  inner  walls  of  this 
pipe  up  to  the  distance  z=-L  within  the  pipe.  The  surface  S  does  not 
include  the  plane  at  z=-L  within  the  pipe.  Let  £  denote  the  sphere 
at  infinity  which  surrounds  the  semi-infinite  pipe  such  that  S  on  the 
exterior  or  outer  wall  of  the  semi-infinite  pipe  is  connected  to  the 
surface  r.  at  z  -*■  -«>  as  in  Figure  26.  Next  consider  the  following 
two  cases.  In  the  first  case,  the  semi-inf i nite  circular 
waveguide  is  excited  from  within  by  the  modal  fields  with  amplitudes 
Amn  and  Rnn,  which  propagate  in  the  z  direction.  However,  in  the 


second  case  the  same  geometry  is  excited  by  an  external  electric  field 
_i  ~  i  -  i 

t  =  0E0  +  which  for  convenience  is  assumed  to  be  produced  by  an 
electric  test  source  TIe  at  Rp  exterior  to  and  in  the  far  zone  of  the 


semi -infinite  pipe. 


It  is  assumed  that  the  spherical  wave  factor 


e-jkRp 

Rp 


in  E  is  suppressed  so  that 
the  field  generated  by  "J0. 


■  represents  a  locally  plane  wave  part  of 
-et  Je  =  (3E0  a-  jf-Ep|  where 


f(|R-Rp|'j  refers  to  Dirac's  delta  function.  The  modal  fields  with 

+  _r  r  ->  r  - 

amplitudes  Anm  and  Bnin  radiate  the  fields  Enm  =Eenm0  +E(j,nm<i>  from  the 

open  end  of  the  pipe  in  the  first  case.  Likewise  the  field  T  which 


is  locally  plane  at  the  open  end  of  the  pipe  launches  the  modal  fields 

with  amplitudes  A"m  and  B~  within  the  semi -infinite  pipe  in  the  second 
case.  As  shown  in  detail  in  Appendix  C,  one  can  apply  the  reciprocity 


theorem  to  the  pair  of  ‘he' ds  c.,,ri  and  f  within  the  region  V  bounded  by 
S,  $0  and  51,  and  tr->n  enforce  the  radiation  condition  on  Z,  as  well  as 
the  boundary  condition  on  S  (e.g.,  tangent i a 1  electric  field  /anishes  on 
the  perfectly  conducting  surface  S).  This  yields  tne  required  modal 

amplitudes  A-  and  3~  in  terms  of  cr  (related  to  [S.„])  as  follows: 
nm  nn  nm 


A- 
m  n 


(9e: 


:  Rp  j.Rp 


jmp 


2A+ 

nm 


I! 


nm 


h1 

nm 


2  ds 


4r 


(164) 


run 


( 9E  o  +  6TM(Rp)Rp.eJkRp 

z  ds 


23* 

nm 


f  f 


e  x  h 
nm  nm 


4n 

jmu 


(165) 


Mere,  the  subscripts  TK  and  TM  in  Cqs.  (164 ) - ( 165 )  denote  that  6p£  is 

*■  J"  + 

radiated  by  Ann  alone,  an j  ljm  is  radiated  by  Bnn  alone,  respectively. 

_r  _r 

ine  radiated  fields  f<;r  Epp  and  Eprvj  have  been  shown  in  Section  (3). 

Because  of  the  normal i zat ion  factor  N  and  N*  introduced  in  Table  1, 

nm  nm 

it  is  noticed  the* 


_( • )  _( ■ ) 

1  //  en;n'  <  h,.r>l  •  z  ds  =  1  .  (166) 

?  S0 


Tne  elements  of  Tspi]  may  be  easily  obtained,  by  inspection,  from 
T>|s.  (ll)-(l?)  and  ( 1 64 ) - ( 165 )  respectively. 


■  ■  I’'--*-:-  - 
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Figure  27  illustrates  an  example  for  the  coupling  into  the  circular 

waveguide  due  to  a  unit  plane  wave  incident  on  axis  with  electric  field 

polarized  in  the  y  direction.  It  is  noted  that  only  TE]m  and  modes 

yield  non-vanishing  radiated  fields  in  the  axial  (0=0)  direction,  thus 

only  TEim  and  TMim  modes  can  be  coupled  into  the  circular  waveguide  for 

the  configuration  of  Figure  27.  The  total  transverse  electric  fields 

evaluated  on  x  (or  y)  axis  are  then  reduced  to  only  E*  (or  Ep),  because 

all  the  coupled  fields  of  Figure  27  possess  a  rcos  -  type  of 

^sin  <y 

variations.  Figure  28  illustrates  the  total  transverse  electric  fields 
on  x  (or  y)  axis  at  various  z=-L  inside  the  circular  waveguide  (with 
ka=35)  due  to  a  unit  plane  wave  incident  on  axis  as  shown  in  Figure  27. 

D.  Ray  Analysis  for  [S22] 

The  scattering  matrix  [S22]  Is  a  modal  reflection  coefficient 

matrix  which  is  associated  with  the  interaction  shown  in  Figure  7  which 

occurs  at  the  open  front  end  (or  z=0).  This  matrix  [S??]  was  defined 

earlier  in  Eqs.  (13)- ( 14 ).  The  order  of  the  sub-matrix  elements  of 
hh  or  ee  he  or  eh 

[S22J;  i.e.,  of  [Rnm  ;  pq  3.  *nd  [Rni„  ;  pq  ],  are  found  by  employing  the 
UTD  in  the  following  manner.  As  before,  the  equivalent  ray  cone 

4-  4- 

corresponding  to  the  modes  with  amplitudes  Apq  and  Rpq  which  are 
incident  on  the  circular  edge  at  the  opening  of  the  semi-infinite  guide 

I 

make  an  angle  5pq  and  <Spq,  with  the  walls,  respectively  (see  Figure 
29).  Then  from  UTD  considerations,  each  of  these  rays  diffracts  hack 
into  the  waveguide  region  at  any  angle  ik  with  the  walls.  Except  for  a 
mode  near  cutoff,  the  angle  i|/will  be  far  from  any  incident  or 


109 


Figure  23.  Total  coupled  transverse  (to  z)  electric  field  along 
x  or  y  at  z=-l  for  a  unit  plane  wave  incident  on-axis 


I 

Figure  29.  Reflection  from  the  open  end  of  a  circular  waveguide. 


3 


reflection  shadow  boundary  that  is  associated  with  the  cone  of  incident 
rays  at  the  edge.  One  may  thus  define  equivalent  magnetic  and  magnetic 
dipole  edge  currents  of  strength  M  and  Mg,  respectively,  on  the  circular 
rim  (or  edge)  which  generate  the  same  diffracted  fields  within  the 
waveguide,  but  near  the  edge,  as  those  produced  by  the  incident  ray 
cone.  Unlike  the  equivalent  magnetic  and  electric  edge  currents  shown 
in  Section  (B)  which  radiate  into  free  space,  the  equivalent  magnetic 
and  magnetic  dipole  currents  mentioned  above  radiate  within  the  guide; 
i.e.,  in  the  presence  of  the  perfectly-conducting  waveguide  walls, 
f orisequently  the  equivalent  electric  edge  current  for  the  exterior  field 
calculation  is  now  replaced  by  the  equivalent  magnetic  dipole  edge 
current  Mg  for  calculating  the  field  in  the  interior  waveguide  region. 
Therefore,  care  must  be  exercised  in  defining  these  equivalent  currents 
wnich  radiate  into  the  waveguide  region.  It  can  be  shown  that  the 
general  form  of  M  is  the  same  as  in  Eq.  (36)  except  that  a  factor  of 
1/2  must  be  included  in  the  definition  of  this  current  and  also  for  Mg 
to  account  properly  for  the  presence  of  the  perfectly-conducting 
waveguide  walls  [28,20]  when  calculating  the  field  in  the  interior 
waveguide  region.  Specifically,  the  equivalent  magnetic  current  is 
given  by. 


_  (  -  „i nc 


>  )  •  •  D^fih  4/'  ;  •  (167) 


modal 

and  the  equivalent  magnetic  dipole  current  is  given  by  [27,28], 
_  ,.inc 


m.U.i  )  ■>  i  ;r 


I  I  V»  „  I - 

.  •  A  )  •  1  .  /81*  • 

modal  ?  >  7 


1  6q  =  ?i  .  (168) 


sin 
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where  4',  =  l$pq  or  Spq.  Clearly  one  requires  4<'  =  ipq  or  fipq  in 
Eqs.  (167)-(168)  corresponding  to  the  angle  made  hy  the  incident  ray 
cone  with  the  waveguide  walls.  Then  the  fields  radiated  into  the 
waveguide  region  by  the  above  rim  currents  may  be  expressed  in  terms  of 
modes  in  the  usual  manner  by  simply  invoking  the  reciprocity  theorem  to 
find  which  modes  are  excited  by  the  source  [20],  i.e.,  by  the  rim 

current.  Since  one  is  interested  in  calculating  the  amplitudes  A* 

or  of  the  reflected  TE _  or  TM„  modes,  respectively,  it  therefore 

follows  from  physical  considerations  that  one  must  require  4/  -  6n,n  or 

<5nm  in  Eqs.  ( 167 ) - ( 168 )  for  these  cases,  based  on  the  aforementioned 
ray-mode  equivalence.  Thus,  one  must  require  that  i|>  take  on  the  proper 
value  corresponding  to  the  ray-cone  angle  6nm  or  <s'1in  which  is 
associated  with  the  reflected  waveguide  mode  when  employing  reciprocity 
to  find  the  strength  of  the  reflected  modes  generated  by  V  and  Mq  of 
Eqs.  (167 ) - (168 )  via  [20].  Without  presenting  any  details,  one  obtains 

the  required  amplitudes  A-  and  B"  via  the  above  procedure  from  which 

the  elements  of  [S22]  can  be  directly  extracted  (by  inspection)  via 
Eq.  (13).  In  particular, 


p  -  ri|n[hnm-^(^nm»  ^pg)  *  ^znm»^  (l<inm»  fyq)  ]dt 

nmlP(l  2 f f  e'  k  h'  •  z  ds 


2 //  e'  x  h'  •  z  ds 
^  nm  nm 

0 


(169) 


- 

rim[hnin.^('5nm»  ®nq)  ]d?-  , 


nm;pq  2 j  j  e_  x  h  •  z  ds 
•  ‘  nm  nm 


(170) 


(171) 


4> 


p  |~i|n[hni[,.^( |Snm> jdi 

'nm;pq  2fje 

So 


x  h 
nm  nm 


z  ds 


h»» 


<6 


p  -  f i 10  jdnjn.M( 6njni  Spg  j  +  ^znm«^  ( 'Snm »  6pq )  1 .  (172) 

n,n : PC1  2/ /  e‘  x  h*  •  z  ds 

nm  nm 

bo 

It  is  important  to  note  that  the  integrals  over  the  rim  in 

Eqs.  (169)-(172)  can  be  evaluated  in  closed  form;  likewise  the 

integrals  over  the  guide  cross  section  S0  (at  z=0)  can  also  be 

evaluated  in  closed  form  (see  Eg.  (166)).  Emptying  the  orthogonality 

properties  of  the  waveguide  modes,  it  can  be  easily  seen  that  the 

numerators  of  Eqs.  (169)-(172)  vanish  when  orn.  Thus,  the  elements  of 

[$2p]  vanish  when  p*n.  Tor  the  dominant  mode  regime,  it  is  necessary  to 

modify  the  above  equivalent  rim  currents  to  include  the  effects  of  rays 

multiply  diffracted  across  the  aperture.  Taking  cognizance  of  the  fact 

that  fhp  multiple  di ff ractions  occur  along  the  reflection  boundaries  ot 

the  interacting  rays,  it  follows  that  multiple  diffraction  calculations 

must  V  treated  carefully  by  decomposing  the  non-ray  optical  fields  near 

tne  reflection  boundaries  into  their  ray  optical  components  as  done 

previously  in  the  calculation  of  [Sn]  at  0=0.  Following  a  similar 

procedure  us  in  'Tppenclix  8  to  include  all  orders  of  rays  multiply 

diffracted  across  the  rim,  one  can  show  that  the  final  expression  for 
rhh 

Knmjnq  which  denotes  the  reflection  coefficient  of  the  TEnm  mode 
that  is  reflected  back  into  the  waveguide  when  a  TEnq  mode  is  incident 
at  /-il,  and  whim  is  based  on  Eg.  (160),  is  given  by 
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hh 


Rnm;nq  “  Ro t 1+u 2 ^ ^nm *  6nq^ 


D 


l+( -1 )nA-D 


] 


+  AoRo  [1~up(^nm’  Vi^  *  - - -  1 

00  c  nm  nq  i+(_i)nAC 


where 


R  =  Nnq  J  n2 _  cos  6* 

o  WI  n?  -  p'  2  J  (p*  )  mi 

nm  rnm  n  v*nm' 


(173) 


[sec  f6nm  -  6nq)  +  sec  (*mn  +  6nq ) ] 


(174) 


.  .2  ,  r  H(2)  '  (p*  ) 

A  =  -jOli)  tan  6 '  tan  6*  tan  6rt(l  tan  nm  ,..n  _ n.cL  • 

oJn  nq  nm  2  2  4  ( 2  V  D  ‘  i 

n  u  nq J  (175) 

n  1  'j(2ka--J-)  ,  /17C, 

u9(i|>>4'  )  =  (-1)  *  e  ^  (sec  <J>  +  sec  ^  )  •  (176) 

2/itka 


n=l 


n  =  l 


(177) 

(178) 


and 


A  =  - 


.-■l(aa-j) 

4/uka 


B  -  -  1  e“j2ka  • 

7 


(179) 

(180) 


ee 

Similarly,  the  final  expression  for  Rnin.nq,  which  denotes  the 
reflection  coefficient  for  the  TMnin  mode  reflected  back  into  the 
waveguide  when  a  TMnq  mode  is  incident  at  z=Q,  is  given  by 
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(181) 


oe 

R  -  R  •  [l  +  u~(5  ,  5„  )  •  H  .  .  1 

rrn;nq  o  1  nm’  nq  1  +  („i)nA.|) 


where 


H(^)(p  )  .  . 

-  _i_  n  nP  sin  5nq 

Nnm  4ka  J '  (p  j~  6nm 
n'Kmn' 


COS  6n,, 


nm 


[sec  [K‘m  ~  °nq)  +  sec  f6nm  +  6ng  ,] 


(182) 


Also,  the  final  expression  for  Rnm;r)q,  which  corresponds  to  the 
reflection  coefficient  for  the  TMnm  mode  which  is  reflected  back  into 
the  waveguide  when  a  TEnq  mode  is  incident  at  i- 0,  is  given  by 


eh 

R  -■  +-  R, 

nm ; nq  i 


u2^nm’ 


6  )  ♦  _ !?_ _ 

nq  i+( -i )nA*n 


(183) 


where 


r,  --  ,Niiri  in/° 

*  fi  4Fa 

run 


«i2,(p 


p  j 

nm  n 


nq )  cos  ^nq 
(Pnm)  cos  6nm 


+  sec 


I 


(184) 


Inn  upper  sign  "r"  in  Fq.  (183)  corresponds  to  the  incident  TEnq  mode  in 
the  upper  row,  arid  reflected  TMnm  mode  in  the  lower  row  of  Table  1.  Tne 
lower  sign  in  Eq.  (183)  corresponds  to  the  incident  TEnm  mode  in  the 
lower  row,  and  reflected  TMnm  mode  in  the  upper  row  of  Table  1. 
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Likewise,  the  final  expression  for  Rnm.nq,  which  corresponds  to 
the  reflection  coefficient  for  the  TEnm  mode  which  is  reflected  back 
into  the  waveguide  when  a  TMnq  mode  is  incident  at  z=0,  is  given  by 


he  ~ 

D  =  +  D 

Knm;nq  ~  *1 


[l+u2( 


D 

l+(-l )nAD 


] 


±  A  R,  [l-u~(5  ,  6)  •  _ G _  ] 

o  1  2  nm  nq  i+(_i)nAC 


where 


(185) 


.(2) 


R,  =  J_  _  -n-  ...  ,  Hn  [Png)  sin  5 
N  4Z  n^  -  p'  2  J  (p  ) 


nq 


nm  n'rnm' 


[sec  (S™  "  6nq)  +  sec  (°nm  *  °nq)] 


'  $nm  +  5ri 


(186) 


H(2)  (p  )  . 

An  =  -j  .n  (!i . tan  «'  cot  6  tan  V  tan  i™. 

0  h(2)  (p  )  nm  nq  T1  — 

n  '-KnqJ 


(187) 


The  upper  sign  "+"  in  Eq.  (185)  corresponds  to  the  incident  TMnq  inode  in 

the  upper  row,  and  reflected  TEnm  mode  in  the  lower  row  of  Table  1.  The 

lower  sign  in  Eq.  (185)  corresponds  to  the  incident  TMnq  mode  in  the 

lower  row,  and  reflected  TEnm  mode  in  the  upper  row  of  Table  1. 

hh  hh  ee  ee 

Plots  of  the  magnitudes  of  R j_|_ ;  1 1 ,  ^01 ; 01  *  Rn ;  1 1 »  and  Roi;01  in 
Eqs.  ( 1 73)- ( 182 )  are  shown  in  Figures  30-31,  respectively,  as  a  function 
of  ka.  These  UTD  based  results  are  also  compared  with  those  based  on 


*b 


(b) 


Figure  30.  Magnitude  of  the  self  reflection  coefficient  for  a  TEnm 
mode  incident  on  an  open-ended  circular  waveguide. 
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Weinstein's  Wiener-Hopf  solution  to  this  problem.  One  observes  from 
these  figures  that  the  agreement  between  the  UTD  and  the  exact 
Wiener-Hopf  solution  for  the  reflection  coefficients  is  good.  It  is 
noted  that  Feisen  and  Yee  [19]  have  obtained  the  modal  reflection 
coefficients  associated  with  acoustic  modes  reflected  from  the  open  end 
of  a  semi-infinite  acoustically  rigid  pipe.  They  [19]  also  employed  ray 
methods,  and  they  used  a  Poisson  sum  technique  to  convert  the  reflected 
ray  cones  into  waveguide  modes.  On  the  other  hand,  the  ray  to  mode 
conversion  for  the  more  general  electromagnetic  case  accomplished  in  the 
present  work  involves  simply  invoking  the  reciprocity  theorem  to 
determine  which  modes  are  excited  by  a  given  source  inside  the 
waveguide;  in  particular,  the  sources  employed  in  this  work  are  the 
equivalent  rim  currents.  While  our  procedure  may  be  somewhat  less 
rigorous  than  the  one  in  [19]  (which  employs  a  Poisson  sum  technique), 
it  appears  to  be  simpler  and  yields  results  which  are  also  quite 
accurate. 

F.  Ray  Analysis  for  [S{j] 

The  scattering  matrix  [SbL  like  [S22L  is  also  a  reflection 
coefficient  matrix  which  is  associated  with  the  discontinuity  at  the 
back  end,  or  at  z=-L,  as  shown  in  Figure  8.  Clearly,  the  elements  of 
[Sb]  depend  on  the  exact  nature  of  che  discontinuity  at  z=-L;  as 
observed  from  Figures  1-4,  the  nature  of  the  discontinuity  is  different 
for  the  different  terminations  at  z=-L  in  those  figures.  Presently,  the 
following  terminations  or  discontinuities  at  z=-L  are  considered; 
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namely,  an  open  circuit  or  open  end  at  z=-L  for  a  finite  pipe,  a  planar 
dielectric  Interface  at  z=-L  In  a  semi-infinite  pipe,  a  disc-blade 
termination  in  a  semi-infinite  pipe,  and  a  short  circuit  or  closed  end 
at  z=-L  for  the  semi -infinite  as  well  as  the  finite  length  pipe, 
respectively. 

a)  [S5]  for  an  Open  Termination  at  z=-L 

When  a  finite  pipe  is  open  at  both  ends,  it  is  apparent  that  [S^] 
is  the  same  as  [S22L  Thus, 

C$b]  =  C$22]  •  (188) 

b)  [Sb]  for  a  Closed  Termination  at  z=-L 

As  shown  in  Figure  3,  there  exists  an  electric  field  which  is 
incident  on  the  termination  from  the  left  side  or  the  region  (z>-l); 

that  field  is  denoted  here  by  E£,  where  E£  may  represent  a  TE  or  TM 
electric  modal  field.  Then  the  total  electric  field  E  in  the  region 
z>-L  can  be  expressed  as 

OD 

?  =  E-  +  l  c*  E*  for  z  >  -L  .  (189) 

k  n=l  "  n 

+ 

where  cn  (n*l,2,...,  »  )  denotes  the  reflection  coefficient  for  the 
+ 

modal  field  Tn  which  Is  reflected  from  the  termination  at  z=-L  when  the 
modal  field  E£  is  incident  there.  One  needs  to  solve  for  the 


--WTwr-tffiiiiirii’ir  t — 
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coefficients  c*  in  Eq.  (189)  for  obtaining  the  elements  of  [$b]  in  this 
case  as  follows.  On  the  perfectly  conducting  termination  at  z=-L,  the 

following  boundary  condition  must  be  satisfied;  namely, 

z  x  r  I  =  0  .  (190) 

(z  =  -L 


+ 

Now  the  unknowns  cn  can  be  determined  quite  easily  from 

Eqs.  (189 ) - ( 190 )  together  with  the  orthogonality  properties  of  the 

waveguide  modes;  thus, 


+  f-1  n=k 

cn  =  1  0  otherwise 


(191) 


It  follows  that 


[Sb]  ^  -[I]  ,  (192) 

where  [I]  is  an  identity  matrix  of  order  «  x 

c)  [S^]  for  a  Planar  Termination  at  z=-l 

As  shown  in  Figure  32,  the  guide  is  loaded  with  a  dielectric 
material  of  relative  permittivity  er  in  the  region  z<-L.  The 
dielectric  interface  at  z=-L  is  planar,  and  perpendicular  to  the  guide 

axis.  Consider  a  modal  electric  field  Z~  incident  upon  this 
termination  at  z=-L  from  the  region  z<-L,  then  the  total  electric  field 
E  in  the  region  z^-L  within  the  guide  can  be  expressed  as. 
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Figure  32.  Backscatterlng  froi 
with  a  dielectric 
er  In  the  region  z 


E.  +  I  c+  E  + 
k  n=l  n  n 


l  l  c-  E- 
n=l  n  n 


-L  <  z  <  0 


z  <  -L  . 


(193) 


likewise  the  total  magnetic  field  H  corresponding  to  E  may  be  expressed 


H.  +  l  c+ 

k  n-1  n  n 


U  c*  H- 
n  =  l  n  n 


-L  <  z  <  0  . 


z  <  -L  . 


(194) 


Here  cn(or  cn)  denotes  the  modal  reflection  (or  transmission) 
coefficient  at  the  junction  z=-L.  It  is  noted  that,  for  -L<z<0, 


En  =  (en  *  enz)  e 


_±  -  —  +J  Y  Z 

Hn  =  cn  (±  hn  +  hnz)  e  nm 


(195) 


(196) 


where  en  (or  TTn ) ,  and  enz  (or  Fnz)  denote  the  transverse  and 
longitudinal  electric  (or  magnetic)  modal  fields,  respectively,  and  can 

^  I 

oe  found  in  Table  1.  ynm  can  be  8nm  or  8nm  to  denote  the  TE  or  TM 
modes.  For  z<-L  (in  the  dielectric  region). 


•  =  c  (e  --  e  )  e 

n  n  '  n  nz' 


''J'Tnm2 


Mn  ■  cn  <*  hn  *  h„z>  e 


T jfnmz 


(197) 


(198) 


where  the  expressions  for  en,  hn,  enz,  and  hnz  can  be  found  from 
'able  1  except  that  the  free  space  permittivity  e  has  to  be  replaced  by 
e  er.  Therefore,  k2  are  replaced  by  k2er,  and  &nm  (or  0nm)  by 
(or  f£m),  where 


e'  =  /k2e  -  for  TE  modes. 

nm  /  r  v  '  nm 


fk^e^  -  (Pfim)^  for  modes. 


(199) 


Then,  by  enforcing  the  boundary  conditions  at  the  dielectric  interface 
which  require  that  z  xT  and  z  x  TT  are  continuous  at  z=-L,  and  employing 
the  orthogonality  properties  of  the  waveguide  modes,  one  obtains. 


~hh  „i 

n  3  Pnm  ~  Pnm 

nm;nm  —  Z7~ 

&  +0 
nm  nm 


(200) 


d  3  pnm  *  er®nm 
nm;nm  ~Z 

6  +  €r0 
nm  r  nm 


(201) 


All  the  other  elements  of  [St,]  vanish  identically.  Thus,  [Sb]  is 
reduced  to  a  diagonal  matrix,  i.e.. 


[R  ]  [  o  ] 

nm;nm 


[SK1  = 


[  0  ) 


[£*e  ) 

nm;nm 


«  x  ® 


(202) 


d)  [Sb]  for  a  Disc-Blade  Termination  at  z=-L 


As  shown  in  Figure  1,  the  planar  disc-blade  discontinuity  comprises 
of  N  blades,  placed  with  an  azimuthal  periodicity  around  a  central  disc. 
The  blades  of  length  c  and  the  disc  of  radius  b  are  assumed  to  be 
perfectly  conducting  and  to  lie  in  the  plane  z=-L.  Each  blade  spans  an 
angle  4»,  and  the  distance  (a-c)  between  the  blades  and  the  cylinder  is 
assumed  to  be  small.  For  a  sufficiently  large  number  of  blades,  the 
surface  current  density  7b  induced  on  the  composite  disc-blade  geometry 
may  be  given  to  a  first  approximation  by 

3|,  ■  2  z  x  ,  on  the  disc  and  blades  only,  (203) 

in  which  TT  1 s  the  sum  of  all  the  propagating  TEnm  and  TMnm  modes 
which  are  incident  on  the  disc-blade  discontinuity.  An  expression  for 
"Hg" is  provided  in  Eq.  (7).  The  approximation  in  Eq.  (203)  corresponds 
to  the  geometrical  optics  approximation  for  the  current.  The  fields 
generated  by  "db  propagate  in  the  +z  direction;  for  propagation  in  the 

A 

+z  direction,  these  fields  constitute  the  reflected  waveguide  modes  in 
the  region  0>z>-L,  and  for  the  -z  directed  propagation  case,  they 
constitute  the  transmitted  waveguide  modes  in  the  region  z<-l.  Once 
again,  as  in  the  calculation  of  the  elements  of  [S22J»  the  reflected  and 
transmitted  waveguide  modes  generated  by  7b  may  be  readily  calculated 
by  invoking  the  reciprocity  theorem  to  determine  which  modes  will  be 
excited  by  a  given  source  within  the  guide  [20].  Thus,  the  amplitudes 

+  f 

Anin  and  Bnm  of  the  TEnm  and  TMnm  modes,  respectively,  which  are 
reflected  from  the  disc-blade  discontinuity  are  found  to  be  [20] 


(204) 


nm 


and 


-//  7* 
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in  which  Sd  denotes  the  surface  area  of  the  disc-blade  geometry.  The 
elements  of  [Sb]ooX<»  can  be  directly  obtained  from  Eqs.  (204)-(205), 
which  are  found  as  follows. 
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(207) 


(208) 


and 


129 


Some  interesting  effects  are  observed  in  the  reflection  from  a 
rotational ly  symmetric  object  inside  a  circular  waveguide.  It  has  been 
shown  [34]  that  if  an  object  is  rotating  about  the  axis  of  the  waveguide 
and  has  periodic  characteristics  in  the  azimuthal  direction  with  period 
?i r/N  (N:  integer),  an  incident  TEpq  (or  TMpq)  mode  can  only  excite 
reflected  TEnm  (or  TMnm)  modes  where  n=p  or  n=JiN±p  [l:  integer). 
Furthermore,  only  the  reflected  modes  withe  n=AN±p,  but  not  n=p,  contain 
information  on  the  modulation  of  the  reflected  wave  due  to  the  rotation 


of  the  periodic  object.  Therefore,  if  the  disc-blade  model  in  Figure  1 
contains  a  large  number  of  fan  blades  (N),  as  is  usually  the  case  in  a 
true  jet  intake,  and  if  the  waveguide  is  not  big  enough  to  propagate  the 
high-order  modes  with  mode  number  n=*N±p,  then  one  will  generally  not  be 
able  to  observe  any  modulation  effect  in  the  reflected  wave. 

In  the  configuration  of  Figure  2  where  the  stator  and  the  rotor  are 
both  taken  into  account,  Eqs.  (206)-(209)  still  apply  except  that  the 
area  of  integration  Sq  appearing  in  the  numerators  is  now  a  function  of 
the  relative  position  of  the  stator  and  rotor.  When  the  geometry  of  the 
rotor-stator  combination  is  no  longer  symmetric,  all  possible  modes  can 
be  reflected  corresponding  to  an  Incident  mode.  However,  the  conversion 
of  an  incident  mode  into  a  different  reflected  mode  is  usually  small 
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compared  with  the  self  reflection  coefficient.  The  modulation  effects 
are  mostly  due  to  the  change  in  the  total  area  covered  by  the  rotor- 
stator  combination. 

Finally,  the  matrix  [P]  which  is  associated  with  the  phase  delay  in 
the  propagation  path  L  can  be  conveniently  written  as 


[P] 


e‘^01L  ..  0  0  -  0  - 

\ 

X 

0  e'jS01L  0 

x 

x- j  e  l 


o 


0  e  nm 
\ 


oo  X  “> 


(210) 
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SECTION  4 


NUMERICAL  RESULTS  AND  SUMMARY  FOR  PART  I 

Numerical  results  for  the  on-axis  RCS  of  the  intake  blade  geometry 
in  Figure  1,  which  are  based  on  the  MSM,  in  conjunction  with  the  UTD 
together  with  its  modi fications  at  caustics,  are  shown  in  Figure  33  as  a 
function  of  a/A  ,  with  the  length  L  (from  the  intake  opening  to  the 
blades)  being  kept  such  that  L=10a.  All  of  the  RCS  plots  here  are 
normalized  to  wa^  .  The  number  of  blades,  N=40  for  the  plot  in  Figure 
33.  The  mean  RCS  level  in  Figure  33  is  basically  controlled  by  the 
blade  scattering,  and  the  ripple  structure  results  from  the  interference 
between  the  blade  backscatter  and  the  backscatter  from  the  intake 
opening.  The  effect  of  evanescent  (non-propagating)  modes  in  the  intake 
is  ignored  in  these  calculations  because  it  is  expected  to  be 
negligible  as  long  as  L  does  not  become  extremely  small  (or  approach 
zero).  The  scattering  matrices  in  the  MSM  become  fairly  simple  in  this 
case  because  first,  they  are  of  finite  order  having  neglected  the 
evanescent  modes,  and  secondly,  there  are  only  a  few  elements  involved 
for  the  range  of  a/A  values  being  considered,  since  only  the  first  few 
propagating  modes  exist  in  this  case.  It  is  interesting  to  compare  the 
intake-blade  RCS  results  of  Figure  33  with  the  RCS  of  two  other  related 
geometries  as  shown  in  Figures  3  and  4,  respectively.  Thus,  Figure  34 
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Figure  33.  On-axis  RCS  (normalized  to  na^)  of  the  simplified  jet 
intake  model  with  N  blades  attached  to  a  circular  disc 
as  shown  in  Figure  1. 
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Figure  34.  On-axis  RCS  (normalized  to  na^)  of  a  hollow,  perfectly 
conducting  finite  length  circular  cylinder  open  at  the 
front  end  and  closed  (shorted)  at  the  back  end. 
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Figure  35.  On-axis  RCS  (non  1ized  to  wa^)  of  a  hollow,  perfectly 
conducting  finite  length  circular  cylinder  open  at 
both  ends. 
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illustrates  the  on-axis  RCS  of  a  hollow,  perfectly-conducting  cylinder 
of  finite  length  L  witn  the  front  end  open  (as  in  the  intake  case)  and 
with  the  back  end  closed  (or  shorted)  with  a  perfectly-conducting  planar 
disk.  Likewise,  Figure  35  illustrates  the  on-axis  RCS  of  the  same 
finite  length  hollow  cylinder  as  in  Figure  4  except  that  it  is  now  open 
at  both  ends.  The  RCS  in  Figures  34-35  have  also  been  calculated  via 
the  MSM  used  in  conjunction  with  the  UTD.  Clearly,  as  one  would  expect, 
fhe  mean  RCS  level  is  higher  for  the  finite  cylinder  with  a  closed 
.’snorted)  oack  end  (Figure  3)  than  for  the  same  cylinder  with  both  ends 
open  (Figure  4).  In  the  case  of  the  finite  cylinder  with  a  closed  back 
end  (Figure  3),  the  mean  RCS  level  is  controlled  by  the  backscatter  from 
the  closed  end,  and  the  ripple  structure  again  results  primarily  from 
the  interference  between  the  backscatter  from  the  closed  back  end  and 
the  backscatter  from  the  open  front  end.  In  contrast,  the  mean  RCS  in 
Figure  35  of  the  finite  length  hollow  cylinder  open  at  both  ends  is 
controlled  primarily  by  the  backscatter  from  the  open  front  end. 

However,  the  ripple  structure  is  primarily  the  result  of  interference 
between  the  backscatter  from  the  open  front  end  and  the  backscatter 
resulting  from  the  wave  coupling  from  the  internal  to  the  external 
regions  of  the  hollow  pipe,  and  vice  versa. 

It  is  important  to  note  that  only  the  TEim  and  TMjm  modes  inside  a 
circular  waveguide  yield  non-vanishing  radiated  fields  in  the  axial 
(0-0)  direction.  Therefore,  from  the  foregoing  discussions  of  the 
scattering  matrix  [S^j]  in  Section  3,  only  the  TE i,n  and  TMim  modes  can 
be  induced  inside  the  circular  waveguide  of  Figures  1-3  when  a  plane 
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wave  is  incident  in  the  axial  (-z)  direction.  Furthermore,  the  TMjm 
modes  yield  very  low-level  fields  in  the  direction  6=0  and  the  mode 
conversion  between  the  TEim  and  TMim  modes  are  small  in  the 
discontinuity  at  z=-L;  for  practical  purposes,  only  the  TE modes  need 
to  be  considered  in  the  calculations  of  the  on-axis  backscattered  fields 
of  Figures  1-3.  On  the  other  hand,  the  TMim  modes  may  become  important 
for  calculating  the  RCS  for  the  problem  in  Figure  4  where  these  modes 
are  coupled  from  the  exterior  region  to  the  interior  at  z=-L  (and  vice 
versa).  These  TMim  modes  have  been  Ignored  in  the  calculations  leading 
to  Figure  35  and  could  possibly  account  for  only  a  part  of  the  almost  3 
dB  shift  between  the  results  based  on  (MSM-UTD  with  modifications)  and 

the  independent  moment-method  calculations  [35]  in  Figure  35  for  the 

\ 

following  reason.  It  is  expected  that  the  mean  RCS  level  in  Figure  35 
should  correspond  to  the  field  backscattered  from  the  open  front  end  as 
accurately  predicted  via  the  (MSM-UTD  with  modifications)  as  shown 
earlier  in  Figure  19;  consequently,  the  effect  of  the  TMjp,  modes  would 
have  to  virtually  be  almost  in  phase  with  the  backscatter  from  the  open 
end  to  yield  a  better  agreement  with  the  moment  method  solution  by 
providing  a  constant  shift  in  the  mean  RCS  level,  but  such  an  in-phase 
behavior  of  the  TMim  modes  may  not  actually  take  place  over  the  entire 
range  of  a/A  values  considered  in  Figure  35  and  will  therefore  not 
completely  explain  the  puzzling  approximate  3  dB  shift  between  the  two 
independent  calculations. 

The  on-axis  solutions  of  Figures  33-34  can  Indeed  be  extended  to  a 
very  large  a/A.  As  an  example,  which  can  be  readily  seen  from  Appendix 
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A,  only  11  modes  (i.e.,  TEj i  through  TE 1 , 1 1 ) *  instead  of  all  the  320 
modes,  need  be  included  when  ka=35. 2187. 

Some  additional  interesting  preliminary  results  of  the  work  are 
shown  in  Figures  36-41.  Figures  36-37  indicate  the  UTD-MSM  based  RCS 
calculations  for  a  finite  or  semi -infinite  hollow  metallic  circular 
cylinder  open  at  the  front  end  and  closed  (shorted)  at  the  termination. 
These  UTD-MSM  based  calculations  are  shown  to  agree  very  well  with 
experimental  results  obtained  from  elsewhere  [13,29]. 

The  computational  speed  can  be  examined  from  Table  2. 

TABLE  2 

NUMBER  OF  MOOES  INCLUDED  VERSUS  CPU  TIME  REQUIRED 


No.  of  modes  included 

Approximate  CPU  time  required  on  a  VAX-11/780 

computer  system  to  generate  the  numerical 

values  of  the  RCS  for  the  problems  in 

Figures  1-4  and  for  a  given  aspect  angle  (e) 

5 

0.13  sec 

10 

0.57  sec 

15 

1.50  sec 

20 

3.01  sec 

25 

6.66  sec 

13(3 


Because  the  required  CPU  time  increases  sharply  as  one  increases  the 
number  of  modes,  some  measures  must  be  taken  to  extend  the  solution  for 
large  values  ka.  To  obtain  the  RCS  at  any  desired  angle,  it  is  possible 
for  the  sake  of  computational  efficiency  to  discard  the  modes  which 
radiate  weakly  in  that  direction.  A  weak  radiated  field  also  means  a 
weak  induced  field  in  the  waveguide  via  reciprocity.  It  is  noted  that 
the  distance  L  between  the  open  end  and  the  termination  of  the  waveguide 
is  so  chosen  that  all  the  evanescent  modes  can  be  negligible. 

Figures  39-41  indicate  the  UTO-MSM  based  calculations  for  the 
modulation  envelope  of  both  the  the  off-axis  and  the  on-axis  RCS  for  the 
configuration  in  Figure  1,  but  with  two  identical  sets  of  planar  blades 
which  are  physically  located  at  the  same  position  inside  the  circular 
intake  duct.  One  set  of  these  blades  is  kept  fixed  as  a  stator,  while 
the  other  set  is  allowed  to  rotate  (i.e.,  it  forms  a  rotor).  The  RCS 
modulation  envelope  then  indicates  the  extent  of  the  change  in  RCS 
resulting  from  the  motion  of  the  rotor  relative  to  the  stator.  Another 
example  is  shown  in  Figure  2  to  illustrate  the  RCS  of  a  semi -inf inite 
hollow  circular  cylinder  with  a  rotor-stator  combination  inside  in  which 
the  stator  and  rotor  blades  are  non-identical.  This  configuration  which 
is  shown  in  Figure  2  resembles  most  jet  intake  geometries.  Consider  the 
case  where  the  stator  and  the  rotor  are  assumed  to  have  60  and  30 
blades,  and  the  various  angular  widths  for  the  stator  and  the  rotor  are 
2°  and  8°,  respectively.  The  stator  blades  are  connected  to  the  intake 
wall,  the  length  of  the  rotor  blades  is  c=0.9a,  and  the  radius  of  the 
center  disk  is  assumed  to  he  b=0.2a.  Since  the  relative  position  of  the 
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Figure  36.  RCS  (normal ized  to  ira^)  of  a  finite  length  hollow 
metallic  circular  cylinder  open  at  one  end  but 
closed  at  the  other  end. 
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Figure  37.  RCS  (normalized  to  *a^)  of  a  finite  length  hollow 
metallic  circular  cylinder  open  at  one  end  but 
closed  at  the  other  end. 
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RCS  (normalized  to  na2)  of  a  semi -infinite  holl 
metallic  circular  cylinder  open  at  one  end  but 
closed  at  the  other  end. 


,RCS  MODULATION 
ENVELOPE  FOR 
ka  ®  3.8,  etc. 


Figure  40.  RCS  (normalized  to  it a^)  modulation  from  a  jet 
Intake  model . 
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Figure  41.  RCS  (normalized  to  ira?)  modulation  from  a  jet 
Intake  model. 
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rotor  blades  with  reference  to  the  stator  blades  is  a  periodic  function 

of  time  (t)  with  a  period  T,  the  reflection  coefficient  matrix  [S5] 

associated  with  that  composite  rotor-stator  blade  structure 

bS 

discontinuity,  and  also,  the  total  backscattered  field  t  has  the  same 

periodic  property.  The  angular  frequencies  present  in  the  backscattered 
_bS 

electric  field  E  are  obtained  in  this  case.  Frequencies  other  than 
those  corresponding  to  the  one  associated  with  the  incident  wave  will  be 
present  in  because  of  blade  rotation.  In  particular,  these 

_bs 

additional  angular  frequencies  in  E  resulting  from  blade  rotation 
(relative  to  the  stator)  comprise  the  so-called  modulation  frequencies, 
and  they  give  rise  to  the  RCS  modulation  frequency  spectrum.  A  short 
analysis  leading  to  the  determination  of  the  RCS  modulation  frequencies 
is  provided  next  as  in  [30].  Let 


^.1 


(t)  e 


(211) 


where  «j0  is  the  angular  frequency  of  the  incident  plane  wave  and  z(t) 
is  a  complex  function  of  t  with  period  T.  It  can  be  seen  that 


T=6 1  •  (212) 

tor 

where  wr  is  the  angular  velocity  with  which  the  blades  rotate. 

z(t)  is  a  periodic  function;  therefore.  It  can  be  represented  by 
a  summation  of  its  Fourier  series, 

«°  j  2  irmt  /T 

t(t)  =  I  T{m)  e  .  (213) 

m=-«° 
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If  | F (m) ! -0,  for  |m|>N0/2,  then  Eq.  (213)  can  be  adequately 
approximated  by 


No/2-1 

j2Timt/T 

T(t)  -  l  T[m)  e 
m=-N0/2 


N0-l  j2ir(m-N0/2)t/T 

*  1  T(m-H0/2)  e 

m=0 

Let 

_  j*N0t/T  _ 
r(t)  e  =  Ej(t) 

and 


(214) 


(215) 


T(m-N0/2)  =  Fx(m)  .  (216) 

Then 

No-l  j2nmt/T 

=l(t)  =  l  Fi(m)  e  .  (217) 

m=0 

The  Inverse  transform  Is 

N0-l  -j2nmt/T 

Fi(m)  «  1  l  Ii(t)  e  .  (218) 

tNo/T=0 

The  ^i(m)  In  Eq.  (218)  can  be  obtained  via  the  Fast  Fourier  Transform 
(FFT)  algorithm  applied  to  the  RHS  of  Eq.  (218).  From  Eqs.  (21 1 )- (213) , 
we  have 


•  J  (<*v>+2irm/T)t 

■  I  F(m)  e 

m=-« 


(2H) 
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The  power  density  spectrum  of  Eq.  (219)  is  defined  by 

T0  ®  j[(w-a)o)-2iTm/T]t  2 

S(oj)  =  1  /  l  T(m)  e  dt 

7T^  -T0  m=- 


=  2  l  7(m)  sinC(ui-uio)-2TTm/T]To 
m=-»  (ui-u)o)-2Ttm/T 


•  l  F* (m 1 )  sin[(ai-mo)-2T!Til/T]T0 

in'=-«  (u)-w0)-2irm,/T 


(220) 


'o  *  "» 


S(u)  =  2n  l  | F (m )  | ^  5[ ( aj-ajo) -2-mn/T]  . 


(221) 


Thus,  the  angular  frequencies  or  spectral  lines  of  t  are  located  at 


=  ojq  +  2irm/T  =  wq  +  mu)p  ;  m=0,  ±1,  ±2, 


(222) 


J)S 

Clearly  mu)p  represent  the  modulation  frequency  spectrum  for  E  .  The 
power  density  spectral  lines  for  the  backscattered  field  at  0=40°  (see 
Figure  2  for  o)  are  shown  in  Figure  42  where  it  is  seen  that  the 
results  for  the  o  and  <j>  polarized  returns  are  almost  identical  for  the 
special  case  of  0=40°. 

Additional  results  for  the  modulation  frequency  spectra  are  shown 
for  various  combinations  of  M  rotor  blades  and  N  stator  blades,  and  for 
different  0,  and  also  different  ka  values  in  Figures  43-46.  The 
dimensions  of  the  rotor  and  stator  blades  as  well  as  of  the  central  disc 
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are  also  shown  In  these  figures.  The  various  angular  widths  for  the 
stator  vane  and  the  rotor  blade  are  all  kept  at  2°  and  8°,  respectively. 
Note  that  the  period  T  depends  on  the  rotor-stator  combination,  and 
except  for  Figure  42,  is  not  necessarily  given  by  Eq.  (212).  It  is  seen 
from  these  RCS  modulation  spectra  that  their  amplitudes  in  certain  cases 
(particularly  for  certain  combinations  of  the  number  of  rotor  and  stator 
blades)  can  be  significant.  Clearly,  a  more  realistic  rotor-stator 
model  would  be  worth  investigating  to  obtain  a  more  realistic  RCS 
modulation  spectrum  of  jet  inlet  configurations.  This  task  is  not 
simple,  but  it  is  hoped  that  based  on  some  of  the  results  of  the  initial 
study  on  the  RCS  of  a  simplied  jet  inlet  model  reported  here,  it  is  a 
topic  which  is  certainly  worth  pursuing  in  the  future.  The  accuracy  of 
the  RCS  analysis  of  the  simplified  models  chosen  in  Figures  1-4  is  based 
in  turn  on  the  good  accuracy  obtained  in  the  DTD  ray  solutions  developed 
here  for  the  elements  of  the  various  scattering  matrices  in  the  MSM 
based  formulation  of  these  problems.  The  good  accuracy  of  the  UTD  ray 
solutions  and  their  modifications  at  caustics  is  established  in  Section 
3  where  the  results  based  on  the  UTO  are  compared  against  results 
based  on  available  exact  (but  more  cumbersome)  solutions  in  some  cases 
and  based  on  measurements  in  other  cases.  It  is  also  noted  that  the  MSM 
has  proved  to  be  quite  accurate  wherever  it  has  been  employed  elsewhere 
for  solving  various  other  scattering  problems  [10];  consequently  the 
present  UTD-MSM  analysis  is  also  expected  to  yield  the  same  type  of 
accuracy  in  dealing  with  the  problems  in  Figures  1-4.  An  advantage  of 
the  present  UTO  ray  representation  for  the  elements  of  the  scattering 
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Figure  42.  Power  density  frequency  spectrum  of  the  polarized 
backscattered  electric  field  for  the  problem  in 
Figure  2  with  0=40°,  ka=3.8,  number  of  stator  vanes=60, 
number  of  rotor  blades~30,  disc  radius  b=0.2af  and  the 
blade  length  c=0.9a. 
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(b)  4>-po1ar1zation 


Figure  42.  (continued) 


Power  density  frequency  spectrum  of  the  backscattered 
electric  field  for  the  problem  in  Figure  2  with 
9=0°,  ka=3.8,  number  of  stator  vanes  =40,  number 
of  rotor  blades=40,  and  disc  radius,  b=0.2a.  The 
blade  length,  c*0.9a. 


-CO  -50  -'50  -30  -20  -10 


o 


WP  =  T 


I  i  l  l 

-4  -3  -2  -1 


i  i  i 

2  3  4 


(  CJ  -  <jJq  )  /  CUp 


(a)  9-polarlzation 


Figure  44.  Power  density  frequency  spectrum  of  the  backscattered 
electric  field  for  the  problem  in  Figure  2  with 
0=0°,  ka=3.8,  number  of  stator  vanes  =20,  number 
of  rotor  blades*40,  and  disc  radius,  b=0.2a.  The 
blade  length,  c*0.9a. 
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igure  45.  Power  density  frequency  spectrum  of  the  backscattered 
electric  field  for  the  problem  in  Figure  2  with 
0=0°,  ka=5,  number  of  stator  vanes*20,  number 
of  rotor  bladesa40,  and  disc  radius,  b=0.2a.  The 
blade  length,  c=0.9a. 
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Figure  46.  Power  density  frequency  spectrum  of  the  backscattered 
electric  field  for  the  problem  in  Figure  2  with 
o=0°,  ka=5,  number  of  stator  vanes*8,  number  of 
rotor  blades=40,  and  disc  radius,  b»0.2a.  The 
blade  length,  c>=0.9a. 
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matrices  required  in  the  MSM  is  that  they  are  not  only  simple  (in 
comparison  to  those  available  from  the  exact  Wiener-Hopf  solution),  and 
quite  accurate,  but  they  also  yield  a  physical  insight  into  the  dominant 
scattering  mechanisms  present  in  the  problems  of  Figures  1-4.  The 
latter  aspect  of  the  UTD  ray  analysis  offers  the  potential  for 
controlling  the  RCS. 

While  it  would  be  worthwhile  extending  the  present  work  to  deal 
with  more  realistic  stator-rotor  models,  it  is  also  worthwhile  as  a  part 
of  future  research  to  study  the  effects  on  the  RCS  of  changing  the 
cross-sectional  shapes  of  the  inlet  duct,  e.g.,  from  the  present 
circular  form  to  rectangular,  skewed  rectangular,  semi-circular,  or 
elliptical  forms.  In  addition,  it  is  also  of  interest  to  consider  the 
continuous  changes  in  the  duct  cross-section  along  the  length  of  the 
inlet  duct,  e.g.,  from  rectangular  to  circular,  or  a  semi-circular  to  a 
circular  one  with  a  bend.  The  latter  cross-sectional  shapes  and  changes 
in  the  cross-section  along  the  length  of  the  inlet  duct  commonly  occur 
in  practice.  A  related  study,  which  has  been  under  investigation  by 
Volakis  [36],  is  on  estimating  the  RCS  of  a  jet  intake  mounted  on  a 
convex  surface.  All  of  the  above  mentioned  problems  are  difficult, 
challenging,  and  worthy  of  future  Investigation  in  the  area  of  RCS 
studies. 
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PART  II 

CANONICAL  RESPONSE  WAVEFORMS  OF  FINITE  AND  OPEN 
LOADED  CIRCULAR  WAVEGUIDES 

Chun-Yue  Lai 
David  L.  Moffatt 
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SECTION  1 


PRELIMINARY  DISCUSSION  FOR  PART  II 

The  interaction  of  electromagnetic  waves  with  the  jet  intake  and 
exhaust  configurations  on  modern  aerospace  vehicles  are  important 
scattering  mechanisms.  While  both  cross  section  modification  effects 
have  oeen  measured,  a  firm  understanding  of  the  scattering  mechanisms 
has  not  followed.  The  purpose  of  this  report  is  to  first  demonstrate 
the  simplicity  and  diagnostic  potential  of  the  canonical  response  wave 
forms  for  a  cavity  structure  and  then  to  postulate  an  approximate  time 
domain  derived  model  to  predict  the  scattering  from  circular  waveguide 
geometries  which  is  suitable  for  extension  to  more  complex  cavities. 

This  report  is  primarily  concerned  with  the  spectral  rather  than 
the  spatial  variation  of  the  electromagnetic  scattering  properties  of 
finite  and  semi-infinite  circular  waveguides  at  guide  diameters  less 
than  3.0  or  4.0  wavelengths.  Consequently,  the  Interrogating  signal 
waveforms  are  limited  to  axial  incidence. 
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In  the  second  section  of  this  report,  the  exact  low  frequency  and 
the  asymptotic  high  frequency  spectral  solutions  to  an  unloaded  semi- 
infinite  open  circular  waveguide  are  summarized.  The  analytical  impulse 
response  waveform  is  also  discussed. 

The  third  section  centers  upon  the  backscatter  from  finite  circular 
waveguides  with  open  or  short  rear  terminations.  The  analytical  model 
chosen  is  the  finite  circular  waveguide  shown  in  Figure  47.  The 
waveguide  has  a  diameter  2a  and  a  length  L.  The  waveguide  walls  are 
assumed  to  be  infinitely  thin  and  perfectly  conducting.  Consider  first 
the  case  of  a  plane  electromagnetic  wave  with  shock  type  time  dependence 
axially  incident  on  the  waveguide.  Regardless  of  the  observation  point, 
the  time-dependent  scattered  waveform  must  be  the  same  as  that  which 
would  be  obtained  from  an  open  circular  waveguide  of  the  same  dimension 
until  such  time  when  the  incident  wavefront  travels  to  the  rear  of  the 
guide  and  the  subsequent  perturbation  then  reaches  the  observer.  It  is 
clear  therefore  that  some  short  time  portion  of  the  canonical  response 
waveforms  of  the  finite  circular  waveguide  can  be  taken  from  those  for 
an  open  or  semi-infinite  circular  waveguide.  Simply  stated,  the 
excitation  cannot  anticipate  the  termination  of  the  waveguide.  Further¬ 
more,  the  evanescent  and  propagating  modes  initially  launched  at  the 
guide  rim  are  independent  of  the  rear  termination.  Therefore  the 
initial  incident  waveguide  modes  on  any  postulated  internal  structure 
(e.g.,  blades)  and  the  subsequent  reradiation  coupling  also  come  from  a 
solution  for  the  open  circular  waveguide.  This  brief  discussion 
illustrates  one  reason  why  a  time  domain  analysis,  where  a  single  real 
time-dependent  waveform  of  relatively  simple  shape  sums  up  the 
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scattering  properties  of  an  object  at  all  frequencies,  is  felt  to  be  a 
fundamental  approach. 

In  the  fourth  section  of  this  report,  the  effects  of  internal 
loading  on  open  circular  waveguides  are  studie'.  The  different  kinds  of 
loading  include  short,  disc  and  blade  structures  (2  and  40  blades). 

They  are  compared  to  the  case  of  a  hollow  semi -inf inite  circular  wave¬ 
guide.  This  involves  the  joining  of  a  high  frequency  and  a  low 
frequency  spectrum  through  a  rational  function  fit,  thereby  enabling  the 
time  domain  responses  to  be  found.  Though  there  is  no  exact  way  to 
prove  the  correctness  of  the  fit,  the  time  domain  waveforms  produced 
thereby  do  satisfy  both  physical  and  moment  conditions. 

In  the  fifth  section,  the  poles  of  the  open  circular  waveguide  are 
explored  and  some  preliminary  results  are  presented.  They  are  based 
both  on  results  in  the  previous  sections  and  on  a  proposed  characteris¬ 
tic  equation  of  the  open  circular  waveguide. 

An  appendix  discusses  the  implications  of  the  research  in  this 
report  on  the  recognition  by  active  radar  sensors  of  targets  with  cavity 
structures  with  resonant  dimension  apertures. 
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SECTION  2 


SEMI -INFINITE  OPEN  CIRCULAR  WAVEGUIDE 


The  formulas  and  results  presented  in  this  section  are  partially 
abstracted  from  a  report  by  Johnson  and  Moffatt  [12]  and  from  a  paper 
[41]  by  the  same  authors.  A  relatively  complete  list  of  significant 
references  to  the  Wiener-Hopf  solution  is  given  in  the  report.  The  open 
circular  waveguide  and  the  coordinate  system  are  shown  in  Figure  48. 

With  the  unfortunate  e"iwt  time  convention,  the  incident  and  scattered 
fields  are  given  respectively  as 


e'  .  jLii  oe;  ♦  . 


(223) 


Es  =  e_"1KL  (0ES.  +  $ESJ  » 


(224) 


where 


(225) 
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Figure  48.  Coordinate  system  for  Wiener-Hopf  solution  to 
semi-inf inite  cylinder. 
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The  exact  defining  integrals  for  the  Wiener-Hopf  factorization 
functions  (L+(k),  M+(k)  in  Equation  (230)*)  are  given  in  Johnson  and 
Moffatt  [12],  In  the  same  reference  these  functions  are  evaluated 
exactly  (diameter/wavelength  less  than  3.0)  using  numerical  integration. 

For  axial  incidence  (9j=0,n=l)  the  infinite  summations  for  the 
polarization  matrix  coefficients  disappear  and 


_  2iacos(4>s)Ji(kasin0s)  L+(k) 

sin(9s)l+(kcos9s)[(2kaM+(k))2-(L+(k))2] 


(232) 


s  I  _  4i ka  cos (<t>s)Ji(kasin9s)  M+(k) 

♦♦  'Q  =Q  D+cos9s]M+(kcos9s)[(2kaM+(k))2-(L+(k))2]  ,  (233) 


i4as1n($s)Ji(kas1n0s)L+(k) 

sin(9s)L+(kcoses)[(2kaM+(k))2-(L+(k))2] 


(234) 


4ik2a3sin(<t>s)Ji(kas1n9s)(M+(k))2 

[1+cos  9s]M+(kcos  9$ )L+(k )[(2kaM+(k ) )?-(L+(k ) )2] 


(235) 


*  more  properly  written  as  L+(k,n),  M+(k,n), 


For  axial  backscatter  (9s=o) 


S°6  !9(.es.o’  L,.es=o'  ' 

(236) 

and,  of  course,  the  cross -p il ari zed  terms  are  zero.  Bowman  [8]  has 
given  an  asymptotic  result  for  the  case  of  axial  backscatter  which 
yields  an  impulse  response  approximation  through  inverse  Laplace 
transform 
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(237) 


where  t0='Vc  is  the  transit  time  for  the  guide  radius  and  the  spatial 
propagation  and  delay  have  been  suppressed.  Beyond  the  reflected 
impulse  the  response  consists  of  alternating  (sign)  paired  causal  and 
non-causal  contributions  delayed  by  integer  multiples  of  the  transit 
time  for  the  guide  diameter.  A  sketch  of  the  response  of  Equation 
(237)  i$  shown  in  Figure  49.  While  the  peaks  predicted  are 
singularities,  the  attenuation  is  severe  because  of  the  weighting  and 
only  the  first  few  terms  contribute  significantly.  A  comparison  of  the 
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IMPULSE  RESPONSE 


impulse  response  of  the  open  circular  waveguide  predicted  using  Fourier 
synthesis  of  the  exact  results  in  Equation  (236)  and  Bowman's  [8] 
asymptotic  approximation  Is  reproduced  from  the  paper  by  Johnson  and 
Moffatt  [41]  in  Figure  50.  Axial  backscatter  cross  sections  were 
also  compared  and  are  reproduced  in  Figure  51.  The  asymptotic 
frequency  domain  results  of  Bowman  (1970)  were  also  obtained  by  Chuang, 
Liang  and  Lee  [11]*.  The  impulse  response  comparisons  in  Figure  50 
show  that  the  major  failure  of  the  asymptotic  waveform  occurs  between 
delays  of  2  to  6  guide  radii  transit  times.  This  translates  into  a 
major  cross  section  error  of  roughly  8.0  dB  at  a  guide  diameter  of  0.59 
wavelengths  which  is  the  cutoff  wavelength  of  the  first  propagating  mode 
( TE 1 1 ) .  Progressively,  as  the  guide  diameter  increases,  the  peaks  of  the 
axial  backscatter  cross  section  occur  at  essentially  the  TE  mode  cutoffs 
(TEji,  TEpi,  TEfli,  etc.)  and  the  cross  section  nulls  are  slightly  in 
advance  (lower  ^/\)  of  the  cutoffs  for  the  TM  modes  (TMoi,  TM^,  TM21 
etc.).  From  Figure  51,  the  asymptotic  result  is  a  reasonable 
approximation  for  guide  diameters  greater  than  one  and  an  excellent  one 
for  guide  diameters  greater  than  two. 

The  simplicity  of  the  time  domain  waveforms  is  emphasized,  despite 
the  complexities  of  the  Wiener-Hopf  solution  the  time  domain  waveforms 
in  Figures  49  and  50  are  extremely  simple.  The  noncausality  of  the 
odd  order  summation  in  Equation  (237)  precludes  direct  utilization  of 
this  model  for  terminated  guides  and  guides  with  internal  structures. 


*The  exact  results  of  Chuang,  Liang  and  Lee  [11]  were  in  error,  as 
reported  in  a  corrections  letter,  Chuang  et  al.  [39]. 
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Figure  50.  Ir.erse  Laplace  transfer™  of  on-axls  backscatter, 
comparing  exact  numerical  results  with  asymptotic 
approximation. 


If  the  approximation 


oo 
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U)m+1[l-(-l)m]u(2(nt0-t)  = 

5/2  ZZZZI 

m  '  2/2mt0-t 


-a _  y 

m=l 


(j)  [l-(-l)m][u(t-2(m-l)t0)-u(t-2mtQ)]  j 

j—  - 

m  2/2mt0-t  (238) 


is  made  then  the  waveforms  are  identical  for  t  _>  2(m-l)t0  (m  is  odd)  but 
a  small  step  discontinuity  results  at  t  =  2(m-l)t0.  The  asymptotic 
waveforms  in  Figures  49  and  50  have  been  terminated  at  zero  but 
actually  are  nonzero  for  negative  time  and  hence  noncausal.  The 
addition  of  a  real  axis  pole,  with  high  damping,  to  remove  the  step 
discontinuity  in  Equation  (238)  would  appear  to  be  feasible.  Instead 
of  the  correction  of  Equation  (238),  we  have 


-  / .  \m+1n  (  /,  -«(t-2(m-l)t0) 

gjtl-za  l  «>  Cl-(-U  3  (l-e 
c  2ir  7*0  m=l 


3/2 

™  2/2'mTo~t 


1  Ju(t-2(m-l)to) 
-u(t-2mtQ)J  _ 


(239) 

where  gc(t)  replaces  the  odd  order  summation  in  Equation  (237).  The 
resulting  time  waveform  with  the  modification  of  Equation  (239)  is 
shown  in  Figure  52.  It  is  found  that  when  ato  is  of  the  order  of  10, 
the  asymptotic  results  In  the  frequency  domain  are  not  altered 
significantly. 

A  comparison  was  made  of  the  radar  cross  section  of  the  circular 
guide  as  given  by  Bowman  [8]  with  that  using  the  modification  of 
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Equation  (239).  The  radar  cross  section  according  to  Bowman  is 
(converted  to  eJwt  time  convention) 


J  l  (-j  )Vj2mka  ,2 

/irka  m=l  3/2 

m 


(240) 


The  radar  cross  section  with  the  modification  of  Equation  (239)  is 


-jw/4 

1  +  • _ 

/irka 


m  even 


(-j) 


me-j2mka 


(-j)m  g-j?mka 

m=1  m3/2  /Jtta 


'+  I 

m= 
m  odd 


1-Erfc(j^j2ka) 


-  e 


-j2mka-2at0  r 


/rr  /jka  +  ato 


1-Erf c(j/2jka+2at0) 
(241) 


In  Table  3,  these  asymptotic  forms  are  compared  to  the  exact  solution 
given  by  Johnson  and  Moffatt  [12].  It  is  evident  from  Table  2-1  that 
for  D/x  >  0.6,  the  cross  sections  corresponding  to  causal  and  noncausal 
time  waveforms  are  equally  good  approximations  of  the  exact  solution. 

Asymptotic  results  obtained  from  a  ray  optical  approximation, 
Bowman  [8],  are  identical  to  Equation  (237)  up  to  m=2  but  thereafter 
differ  in  the  weights  of  the  terms.  The  asymptotic  result  in  Figure  50 


177 


177 


TABLE  3 


COMPARISON  OF  RADAR  CROSS  SECTION  (dB)  OF  EXACT  SOLUTION 
(JOHNSON),  ASYMPTOTIC  SOLUTION  (BOWMAN)  AND  MODIFIED 
ASYMPTOTIC  SOLUTION  (Eqn.  2-19,  MOFFATT) 


D/X 

Johnson 

Bowman 

Moffatt 

O.lOO 

-9.458848 

o.no 

-8.520934 

0.120 

-7.646716 

0.130 

-6.824223 

0.140 

-6.045925 

0.150 

-5.305766 

0.160 

-4.595446 

0.170 

-3.913079 

0.180 

-3.255218 

0.190 

-2.617804 

0.200 

-2.001531 

0.210 

-1.408495 

0.220 

-0.833483 

0.230 

-0.276990 

0.240 

0.264939 

0.250 

0.792181 

0.260 

1.301408 

0.270 

1.795182 

0.280 

2.270724 

0.290 

2.731840 

0.300 

3.174628 

0.310  •„ 

3.627388 

0.320 

4.064384 

0.330 

4.490535 

0.340 

4.902714 

0.350 

5.301191 

0.360 

5.686905 

0.370 

6.056279 

0.380 

6.408279 

0.390 

6.744378 

0.400 

7.066416 

0.410 

7.288048 

0.420 

7.530067 

0.430 

7.786634 

0.440 

8.061458 

0.450 

8.346395 

0.460 

8.636331 

0.470 

8.927008 

0.4C0 

9.213129 

0.490 

9.493250 

0.500 

9.743763 

0.510 

9.976593 

0.520 

10.161530 

0.530 

10.306625 

0.540 

10.395186 

0.550 

10.413555 

0.560 

10.367692 

0.570 

10.254469 

0.580 

9.973226 

0.590 

8.092514 

0.600 

6.002594 

0.610 

4.758259 

3.071432 

2.880305 

0.620 

3.745943 

3.179154 

3.006446 

0.630 

2.918265 

3.244816 

3.090103 

0.640 

2.248425 

3.336831 

3.203558 

0.650 

1.711760 

3.476034 

3.363394 

0.660 

1.291267 

3.585063 

3.488144 

0.670 

0.976372 

3.642307 

3.560661 

0.680 

0.755607 

3.740452 

3.678051 

0.690 

0.615962 

3.907654 

3.066302 

0.700 

0.557735 

4.030650 

4.003542 

0.710 

5.867754 

4.074891 

4.060959 

0.720 

S.61S5S6 

4.203629 

4.211191 

0.730 

5.300941 

4.496467 

4.528727 

0.740 

4.915490 

4.701146 

4.745268 

0.750 

4.462694 

4.495789 

4.531564 

0.760 

3.960338 

3.819680 

3.829793 

0.770 

3.412798 

2.965487 

2.946289 

0.780 

2.842850 

2.349164 

2.317832 

0.790 

2.256176 

2.000432 

1.971990 

0.800 

1.679505 

1.635842 

1.605475 

0.810 

1.117669 

1.187165 

1.147456 

0.820 

0.591479 

0.819100 

0.775175 

0.830 

0.107829 

0.564319 

0.523153 

0.840 

-0.325362 

0.281201 

0.238399 

0.850 

-0.703480 

-0.053583 

-0.104108 

0.860 

-1.026303 

-0.324297 

-0.378683 

0.870 

-1.288107 

-0.521342 

-0.574454 

0.880 

-1.497142 

-0.749196 

-0.805479 

0.890 

-1.655850 

-1.009100 

-1.073821 

0.900 

-1.767546 

-1.208304 

-1.277504 

0.910 

-1.839869 

-1.354398 

-1.423762 

0.920 

-1.882918 

-1.531832 

-1.606125 

0.930 

-1.901502 

-1.726239 

-1.809954 

0.940 

-1.899772 

-1.861168 

-1.949681 

0.950 

-1.890381 

-1.957783 

-2.047384 

0.960 

-1.877861 

-2.005210 

-2.181114 

0.970 

-1.866457 

-2.217362 

-2.323117 

0.980 

-1.861093 

-2.292324 

-2.402434 

0.990 

-1.867049 

-2.341610 

-2.453027 

1.000 

-1.880707 

-2.420392 

-2.538843 

1.010 

-1.907  961 

-2.494346 

-2.622281 

1.020 

-1.937236 

-2.515522 

-2.646585 

1.030 

-1.976914 

-2.522350 

-2.654363 

1.040 

-2.021478 

-2.556740 

-2.695997 

1.050 

-2.065371 

-2.579477 

-2.727279 

1.060 

-2.106021 

-2.555760 

-2.705028 

1.070 

-2.140774 

-2.527534 

-2.G77213 

1.080 

-2.161213 

-2.524461 

-2.681547 

1.090 

-2.173288 

-2.505359 

-2.670043 

1.100 

-2.165415 

-2.447143 

-2.611728 

1.110 

-2.146071 

-2.392444 

-2.557278 

1.120 

-2.109116 

-2.360129 

-2.533267 

1.130 

-2.059807 

-2.309289 

-2.489650 

1.140 

-1.994752 

-2.226755 

-2.406308 

1.150 

-1.919178 

-2.153824 

-2.334631 

1.160 

-1.835327 

-2.100295 

-2.292448 

1.170 

-1.745910 

-2.027018 

-2.228523 

1.180 

-1.656677 

-1.928881 

-2.129956 

1.190 

-1.570702 

-1.844659 

-2.051133 

1.200 

-1.484781 

-1.776799 

-2.004253 

1.210 

-1.408654 

-1.688724 

-1.934084 

1.220 

-1.339279 

-1.501738 

-1.831867 

1.230 

-1.277167 

-1.491552 

-1.767274 

1.240 

-1.222956 

-1.414743 

-1.771731 

1.250 

-1.173344 

-1.317664 

-1.782492 

1.260 

-1.131993 

-1.206690 

-1.732754 

1.270 

-1.089466 

-1.114349 

-1.629171 

1.280 

-1.045954 

-1.032542 

-1.506572 

1.290 

-1.001394 

-0.930503 

-1.378139 

1.300 

-0.949340 

-0.818919 

-1.253332 

1.310 

-0.892004 

-0.727010 

-1.142442 

1.320 

-0.821347 

-0.642964 

-1.033542 

1.330 

-0.741740 

-0.538906 

-0.910714 

1.340 

-0.646153 

-0.428665 

-0.787844 

1.350 

-0.541254 

-0.338876 

-0.682887 
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1.360 

-0.426064 

1.370 

-0.301055 

1.380 

-0.168339 

1.390 

-0.035823 

1.400 

0.104695 

1.410 

0.237799 

1.420 

0.367332 

1.430 

0.484647 

1.440 

0.597126 

1.450 

0.692532 

1.460 

0.781571 

1.470 

0.859947 

1.480 

0. 930879 

1.490 

0.995187 

1.500 

1.055551 

1.510 

1.116309 

1.520 

1.181357 

1.530 

1.251361 

1.540 

1.343095 

1.550 

1.441111 

1.560 

1.557706 

1.570 

1.694245 

1.580 

1.843339 

1.590 

2.005405 

1.600 

2.185416 

1.610 

2.370679 

1.620 

2.560735 

1.630 

2.746525 

1.640 

'  2.916111 

1.650 

3.073452 

1.660 

3.206708 

1.670 

3.309040 

1.680 

3.373470 

1.690 

3.398627 

1.700 

3.303114 

1.710 

3.316839 

1.720 

3.205293 

1.730 

3.044869 

1.740 

2.833142 

1.750 

2.592613 

1.760 

2.303109 

1.770 

1.991521 

1.780 

1.658170 

1.790 

1.305361 

1.800 

0.957284 

1.810 

0.610778 

1.820 

0.279078 

1.830 

-0.032145 

1.840 

-0.311415 

1.850 

-0.558779 

1.060 

-0.772689 

1.870 

-0.951250 

1.880 

-1.097941 

1.890 

-1.204541 

1.900 

-1.290164 

1.970 

-1.342392 

1.920 

-1.369620 

1.930 

-1.383562 

1.940 

-1.391595 

1.950 

-1.385437 

1.960 

-1.376858 

1.970 

-1.372984 

1.980 

-1.373787 

-0.254430 

-0.58040 1 

-0.150045 

-0.461617 

-0.042492 

-0.343499 

0.044128 

-0.244776 

0.127807 

-0.147102 

0.231507 

-0.031908 

0.335697 

0.061024 

0.418518 

0.173686 

0.500756 

0.267246 

0.603329 

0.379283 

0.703902 

0.487138 

0.782561 

0.573963 

0.863047 

0.663914 

0.964523 

0.773385 

1.061533 

0.876530 

1.135011 

0.957706 

1.214536 

1.044577 

1.315421 

1.152379 

1.409195 

1.251371 

1.478916 

1.327015 

1.556167 

1.411466 

1.657583 

1.519011 

1.748774 

1.614656 

1.813689 

1.684713 

1.890157 

1.767725 

1.994319 

1.077441 

2.084119 

1.971007 

2.143686 

2.035024 

2.220860 

2.118327 

2.332511 

2.235130 

2.423400 

2.329006 

2.476077 

2.385360 

2.557715 

2.472930 

2.689892 

2.609956 

2.786655 

2.710307 

2.827770 

2.751536 

2.929719 

2.859467 

3.150998 

3.004904 

3.290674 

3.218311 

3.093302 

3.003269 

2.531862 

2.419261 

1.873926 

1.745675 

1.433071 

1.302936 

1.187997 

1 .060206 

0.921487 

0.791057 

0.604546 

0.470287 

0.359371 

0.225954 

0.19129S 

0.060154 

-0.002104 

-0.134483 

-0.224069 

-0.358561 

-0.393827 

-0.527132 

-0.516608 

-0.648159 

-0.664207 

-0.796644 

-0.827888 

-0.961013 

-0.946561 

-1.079222 

-1.033748 

-1.164748 

-1.144185 

-1.276220 

-1.262092 

-1.395434 

-1.339409 

-1.471423 

-1.395425 

-1.526003 

-1.473282 

-1.605004 

-1.552135 

-1.684992 

-1.594074 

-1.725301 

predicts  the  positions  of  similarities  and  discontinuity  quite  well  but 
fails  to  predict  the  damped  sinusoid  corresponding  to  the  large  error  in 
cross  section  at  the  cutoff  frequency  of  the  first  propagating  mode 
TEn. 

If  the  open  circular  waveguide  Is  shorted  internally  at  a  distance 
L  down  the  guide  then  each  propagating  mode  Is  reflected  as  the  same 
propagating  mode.  For  axial  Incidence  the  contribution  of  the  TEu  mode 
to  the  backscatter  is,  according  to  Johnson  and  Moffatt  [41], 

ES  =  _ei2aL8ueC0HE1  ,  (242) 

where  the  coupling  coefficient  is 

B  0  «  2i  (k+qii)M+(aii)  fl  ,  (243) 

<hi(1  --UV")  wf 

J11 

the  radiation  coefficient  Is 

I 

r  =  +(ka)2  M«ll)  1  fl 

*  L+(k)  2  .  <224> 

and 


°11  = 


(225) 


I 

The  J 1 1  is  the  first  (m»l)  zero  of  the  derivative  of  the  Bessel 


function,  i.e.,  1.84118 


SECTION  3 

FINITE  CIRCULAR  WAVEGUIDES  WITH  OPEN  OR  SHORT  REAR  TERMINATIONS 

The  axial  backscatter  by  the  finite  circular  waveguide  shown  in 
Figure  47  has  been  computed  using  a  moment  method  program  for  rotation- 
ally  symmetric  objects.*  Both  open  and  shorted  rear  termination  of  the 
waveguide  were  used.  The  normalized  axial  radar  cross  sections  of  both 
the  open  and  shorted  finite  guides,  a  thin  circular  disk  of  radius  a  and 
the  semi-infinite'  open  waveguide  (a  few  points)  are  shown  In  Figure  53. 
Several  features  should  be  noted.  First,  for  guide  diameters  less  than 
the  cutoff  wavelength  of  the  first  propagating  mode  (TEn)  the  radar 
cross  sections  of  the  open  and  shorted  cylinders  are  essentially  the 
same,  l.e.,  the  rear  termination  is  unimportant.  This  would  also  be 
true  for  any  structure  in  the  interior  of  the  guide  provided  they  are 
sufficiently  removed  from  the  front  of  the  guide  (evanescent  modes). 
Second,  the  disk  cross  section  in  this  same  range  (D/x  <  0.6)  is  a  good 

♦These  computations  were  made  by  Dr.  W.S.C.  Chang  of  the  ElectroScience 
Laboratory,  Department  of  Electrical  Engineering,  The  Ohio  State 
University,  Contract  No.  N00014-78-C-0049. 
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Figure  53.  Normalized  axial  radar  cross  section  of  finite 
hollow  cylinders  (open  and  shorted  at  far  end), 
a  circular  disk  and  a  semi-infinite  circular 
waveguide. 


approximation  for  the  average  of  the  finite  guide  cross  sections. 

Third,  the  disk  cross  section  is  a  reasonable  approximation  for  the  open 
waveguide  cross  section  in  this  same  range. 

Using  computations  spanning  a/X=0. 01 (0.01)1.20  the  axial  impulse 
response  waveforms  of  the  open  and  shorted  finite  cylinders  have  been 
generated  using  a  Fourier  synthesis  procedure,  Kennaugh  and  Moffatt 
[43].  These  response  waveforms  are  shown  in  Figure  54  for 
respectively  the  open,  54a,  and  shorted,  54b,  guides.  The  initial 
response,  impulse  plus  ringing  (0<t/to<20),  is  the  same  for  both  guides 
as  was  anticipated  in  the  Introduction.  This  response  is  also  the  same 
as  that  for  the  open  guide  in  Figure  50,  but  for  the  open  guide  much 
more  resolution  was  obtained  because  computations  to  a/x  of  6.0  were 
used.  The  negative  Impulse  at  the  origin  was  removed  in  these  and 
subsequent  figures.  The  basic  scattering  mechanisms  of  the  finite 
guides  can  be  enumerated,  some  of  which  are  isolated  and  identifiable  in 
the  waveforms  of  Figure  54.  Recall  that  both  the  finite  and 
semi-infinite  guides  have  an  Infinitely  thin  wall  thickness. 

1.  Specular  (frequency  Independent)  return  which  is  a  weighted 
impulse  at  Vto=0. 

2.  Rim  scattering  or  diffraction  back  and  forth  across  the  front 
rim.  This  Is  the  ringing  for  0<Vto<20.0. 

3.  Singly  diffracted  return  from  the  rear  of  the  guide.  This 
excitation  and  return  travels  at  approximately  the  free  space 
velocity  exterior  to  the  guide  and  is  clearly  identifiable  as 
the  spike  (bandlimlted  approximation  of  an  impulse)  at 
t/to*20.0. 
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Figure  54b.  On-axis  backscatter  Impulse  response  from  finite 
circular  cylinder  shorted  at  the  rear. 
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4.  Diffraction  across  the  rear  of  the  guide  or  a  doubly  diffracted 
contribution,  which  is  masked  with  other  effects  in  Figure  54. 

There  are,  of  course,  also  evanescent  waveguide  modes  for  both  guides 
which  are  unimportant  for  these  geometries.  The  above  four  mechanisms 
are  identical  for  both  guides. 

5.  Propagating  modes  are  coupled  into  the  guide  structure, 
reflected  at  the  termination  and  then  reradiated  at  the  front 
of  the  guide.  The  reflection  coefficient  at  the  rear  is 
obviously  different  for  the  open  and  shorted  guides.  Note, 
however,  that  there  is  no  mode  conversion  for  these 
terminations.  In  time,  these  mechanisms  must  occur  later  than 
t/to=20.0  because  of  the  mode  energy  velocity  in  the  guide. 

6.  For  the  open  finite  guide  another  scattering  mechanism  exists 
which  has  two  reciprocal  forms.  Coupling  into  the  guide  at  the 
front  as  waveguide  modes,  guide  propagation  and  then 
diffraction  at  the  rear  of  the  guide  and  free  space  propagation 
in  the  backscatter  direction.  Obviously  the  reverse  or 
reciprocal  mechanism  also  occurs. 

A  principal  advantage  of  the  time  domain  viewpoint  is  that  the 
above  scattering  mechanisms  can,  in  principle,  be  separately  removed 
from  the  total  response  waveform.  More  importantly,  approximate  impulse 
response  waveforms  for  these  and  similar  cavity  geometries  can  be 
constructed  by  simple  addition  of  basic  waveforms  with  appropriate 
delays.  The  impulse  at  the  waveform  origins  in  Figure  54  and  Figure 
50  is  simply  predicted  by  asymptotic  theory.  In  fact,  the  impulses 
shown  have  been  added  artificially.  The  weight  of  this  impulse.  Equation 
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(237),  was  subtracted  from  the  frequency  domain  data  before  Fourier 
synthesis  in  order  to  improve  details  of  the  waveforms  for  short  times. 
In  a  similar  fashion,  the  singly  diffracted  contribution  which  is 
easily  predicted  via  the  geometrical  theory  of  diffraction  (GTD)  as  * 
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can  be  removed. 

The  impulse  response  waveforms  of  the  finite  open  and  shorted 
guides  with  the  specular  contribution  removed  are  shown  in  Figure  3-3. 

It  remains,  in  this  time  span,  to  approximate  the  sinusoidal  delay  for 
0<t/to<20.0.  This  type  of  approximation  is  discussed  later. 

UTD  provides  excellent  asymptotic  frequency  domain  approximations 
for  the  doubly  diffracted  contribution.  Unfortunately  the  finite  guides 
have  a  cylindrical  caustic  on  the  guide  axis.  The  proper  constant  phase 
shift  correction  at  all  frequencies  for  this  caustic  results  in  a 
nonphysical  system  and  consequently  a  noncausal  response  in  the  time 
domain.  Thus  if  this  term  Is  subtracted  before  synthesis  a  distortion 
of  the  resulting  waveform  at  times  before  this  contribution  can  poscibly 
arrive  results.  This  is  illustrated  in  Figure  56  for  the  finite 
shorted  guide.  Comparing  Figure  56  and  Figure  55b  it  is  seen  that  a 
distortion  before  t/to=20.0  has  occurred  and  this  cannot  occur 


♦Modern  Unified  Theory  of  Diffraction  (UTD),  Pathak  and  Kouyoumjian 
offers  more  sophisticated  results  but  the  simpler  GTD  formula  is 
adequate  here. 


IMPULSE  RESPONSE 
SHORTED  FINITE  CYLINDER 

DIRECT  DIFFRACTION  FROM 
REAR  REMOVED 
DOUBLE  DIFFRACTION  FROM 
REAR  REMOVED 
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physically.  The  asymptotic  phenomenon  occurring  was  shown  in  Figure  49 
for  a  different  type  of  response.  Note  carefully  that  In  time  the 
doubly  diffracted  contribution  physically  begins  to  occur  at 
t/to=20.0+e.  According  to  the  asymptotic  theory  it  Is  identically  zero 
for  Vto>22.0.  It  is  precisely  this  caustic  correction  problem  which 
prevents  a  large  number  of  GTD  and  UTD  results  from  being  fully 
exploited  in  the  time  domain.  The  difficulty  has  been  known  for  some 
time,  Moffatt  [53],  but  a  simple  effective  solution  has  not  ensued. 

For  the  finite  shorted  guide  we  can  consider  removing  mechanism  5, 
the  propagating  mode  contributions,  by  utilizing  results  from  the 
solution  for  the  open  circular  waveguide  with  a  shorting  plate  as  shown 
in  Figure  57.  This  contribution  must  be  the  same  for  the  two 
structures.  The  axial  backscatter,  based  on  coupling  and  reradiation  of 
the  TEji  mode  and  neglecting  the  rim  scattering  was  given  in  Equation 
(242).  Radar  cross  section  calculations  of  the  open  waveguide  in 
Figure  57  which  include  the  TEj^  and  TMji  modes  and  the  rim 
contribution  have  been  made,  Johnson  and  Moffatt  [41],  Unfortunately 
the  propagating  mode  contributions  were  not  separated  and  the  time 
waveform  corresponding  to  Equation  (242)  will  require  some  repeated 
calculations.  Since  these  Involve  the  Wiener-Hopf  factorization 
functions  they  are  not  simple. 

An  extensive  table  of  values  of  the  factorization  functions  L+(a) 
and  M+(a)  computed  by  numerical  Integration  is  found  in  Johnson  and 
Moffatt  [12].  We  are  now  Interested  in  the  case  of  n=l  (backscatter). 

For  0.1<0/x<0.7,  the  L+  and  M+  functions  vary  smoothly  except  at  a  cusp 
In  the  imaginary  part  of  the  M+  function  at  rVx=0.5861.  An  interpolation 
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program  was  written  to  obtain  values  of  the  factorization  functions  at 
Inbetween  frequencies.  It  avoids  interpolation  across  the  cusp  which 
would  lead  to  erroneous  results.  The  results  for  l+(k)  and  M+(k)  for 
the  n=l  case  are  given  In  Table  4. 

Based  on  the  simplest  asymptotic  forms,  the  axial  backscattered 
field  is  given  by,  Johnson  and  Moffatt  [41] 

rs  =  f  -sa2  *  a(4+jil2)  1  „  '2sL/c 

L^I  UixVlJBs  VcJ  •  {W) 

The  time  waveform  corresponding  to  Equation  (3-2)  is 


F,(t)  - 


(t-2L) 

c 


u(t-2L)  . 

c 

(248) 


The  most  obvious  difficulty  with  Equation  (248)  is  the  delay--we  know 
that  it  should  be  slightly  greater  than  21 /c.  The  precise  inversion  of 
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TABLE  4 

RESULTS  FOR  L+(k)  AND  M+(k)  FOR  THE  n=l  CASE. 


L+(k>  H+(k) 


o/x 

r/kI 

Real 

Imag. 

Real 

Imag. 

.100 

1.0 

0.9595 

0.1825 

0.3033 

1.6039 

.110 

1.0 

0.9523 

0.1984 

0.3020 

1.4605 

.120 

1.0 

0.9446 

0.213b 

0.3007 

1.3411 

.130 

1.0 

0.9362 

0.2207 

0.2994 

1.2403 

.140 

1.0 

0.9273 

0.2433 

0.2981 

1.1541 

.150 

1.0 

0.9179 

0.2574 

0.2968 

1.0790 

•  ICO 

1.0 

0.9078 

0.2711 

0.2954 

1.0145 

.170 

1.0 

0.8972 

0.2843 

0.2941 

0.9572 

.180 

1.0 

0.8861 

0.2971 

0.2928 

0.9064 

.190 

1.0 

0.8743 

0.3095 

0.2915 

0.8611 

.200 

1.0 

0.8620 

0.3214 

0.2902 

0.8205 

.210 

1.0 

0.8510 

0.3324 

0.2884 

0.7829 

.220 

1.0 

0.8398 

0.3429 

C .2866 

0.7487 

.230 

1.0 

0.8284 

0.3528 

0.2848 

0.7175 

.240 

1.0 

0.8167 

0.3622 

0.2829 

0.6088 

.250 

1.0 

0.8048 

0.3711 

0.2810 

0.6623 

.200 

1.0 

0.7927 

0.3795 

0.27S1 

0.6379 

.270 

1.0 

0.7804 

0.3873 

0.2771 

0.6152 

.280 

1.0 

0.7678 

0.3946 

0.2752 

0.5942 

.290 

1.0 

0.7550 

0.4014 

0.2732 

0.5745 

.300 

1.0 

0.7420 

0.4077 

0.2712 

0.5561 

.310 

1.0 

0.7300 

0.4139 

0.2696 

0.5378 

.320 

1.0 

0.7181 

0.4196 

0.2679 

0.5205 

.330 

1.0 

0.7061 

0.4248 

0.2664 

0.5040 

.340 

V  1.0 

0.6941 

0.4296 

0 . 26  46 

0.4833 

.350 

1.0 

0.6822 

0.4340 

0.2632 

0.4733 

.360 

1.0 

0.6702 

0.4379 

0.2617 

0.4569 

.370 

1.0 

0.6583 

0.4413 

0.2602 

0.4451 

.380 

1.0 

0.6464 

0.4443 

0.2587 

0.4319 

.3  90 

1.0 

0.6344 

0.4469 

0.2573 

0.4192. 

.400 

1.0 

0.6225 

0.4490 

0.2556 

C.40CP 

.410 

1.0 

0.6112 

0.4513 

0.2530 

C . 3972 

.420 

1.0 

0.6001 

0.4533 

0.2509 

0.3870 

.430 

1.0 

0.5690 

0.4549 

0.2494 

0.3764 

.440 

1.0 

0.5781 

0.4563 

0.24C6 

0.3653 

.450 

1.0 

0.5673 

0.4572 

0.2483 

0.3538 

.460 

1.0 

0.5566 

0.4579 

0.2485 

0.3420 

.470 

1.0 

0.5460 

0.4582 

0.2493 

0.3296 

.480 

1.0 

0.5356 

0.4582 

0.2505 

0.3173 

.490 

1.0 

0.5253 

0.4579 

0.2521 

0.3044 

.500 

1.0 

0.5151 

0.4572 

0.2542 

0.2913 

.510 

1.0 

0.5053 

0.4568 

0.2566 

0.2779 

.520 

1.0 

0.4956 

0.4562 

0.2595 

0.2643 

.530 

1.0 

0.4861 

0.4553 

0.2626 

0.2504 

.540 

1.0 

0.4768 

0.4543 

0.2661 

0.2362 

.550 

1.0 

0.4676 

0.4530 

0.2699 

0.2219 

.500 

1.0 

0.4586 

0.4515 

0.2740 

0.2073 

.570 

1.0 

0.4498 

0.4497 

0.2783 

0.1926 

.500 

1.0 

0.4411 

0.4478 

0.2895 

0.1693 

.590 

1.0 

0.4326 

0.4456 

0.3190 

0.1553 

.600 

1.0 

0.4243 

0.4432 

0.3523 

C .  164  8 

.610 

1.0 

0.4162 

0.4410 

0.3720 

0.1720 

.620 

1.0 

0.4083 

0.4387 

0.3882 

0.1792 

.630 

1.0 

0.4005 

0.4362 

0.4011 

0.1665 

.640 

1.0 

0.3929 

0.4336 

0.4108 

0.1938 

.650 

1.0 

0.3855 

0.4308 

0.4175 

0.2011 

.660 

1.0 

0.3703 

0.4279 

0.4213 

0 . 2085 

.670 

1.0 

0.3713 

0.4249 

0.4223 

0.2159 

.600 

1.0 

0.3644 

0.4217 

0.4206 

G.2234 

.690 

1.0 

0.3577 

0.41C4 

0.4165 

0.2309 

.700 

1.0 

0.3512 

0.4149 

0.4099 

0.23  84 

This  waveform  would  have  to  be  convolved  with  the  inverse  transform  of 
the  b,9,C  u  expression  In  Equation  (242)  to  yield  the  true  Impulse 
response  of  the  coupled  and  reradiated  mode.  The  impulse  term  in 
Equation  (250)  is  a  precursor  Illustrating  that  even  for  guided  waves 
the  first  presence  of  the  signal  (perfect  detector)  arrives  with  the 
velocity  of  light  in  a  vacuum.  The  signal  front  buildup  is  dictated  by 
the  second  term  in  Equation  (250).  For  guided  waves  the  topic  of 
transients  was  discussed  by  Cerlllo  [37],  The  present  case  Is  further 
complicated  by  the  coupling  and  radiation  expressions,  which  are  also 
frequency  dependent. 

The  doubly  diffracted  contribution  in  the  time  domain  for  the 
finite  shorted  cylinder  is  given  by  (asymptotic  estimate) 


f  DOft>  =  -2  u(22t0-t) 

1  mr-t 

o 


(251) 


According  to  this  estimate  the  contribution  endures  only  for  t/toc22.0 
(t0-a/c).  The  essential  truth  of  this  can  be  seen  by  comparing  Figures 
54b  and  56.  For  to/t>22.0,  the  response  waveforms  are  identical. 
Therefore  for  t/to>22.0,  the  waveforms  in  Figures  54b,  55b  and  56  are 
the  coupling,  guide  propagation  and  radiation  of  basically  the  TEn  and 
TMu  circular  waveguide  modes.  At  the  smallest  wavelength  calculated, 
a/X=1.2,  there  are  actually  10  TE  and  6  TM  modes  which  can  propagate  but 
for  axial  backscatter,  only  the  TEn  and  TMji  modes  are  of  any  real 
consequences.  At  other  aspects  or  for  terminations  which  cause  mode 
conversion  this  will  no  longer  be  true.  We  postulate  therefore  that  on- 
axis  the  backscatter  contribution  from  coupling  and  radiation  consists 
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of  a  weighted  doublet  delayed  by  2L/(energy  velocity  of  TEn  mode)  and 
then  damped  sinusoid  contributions  which  originate  at  times  dictated  by 
the  energy  velocity  of  the  higher  order  modes.  Impulse  and  step 
response  waveforms  due  to  the  short  inside  the  semi-infinite  guide  are 
obtained  from  Inversion  of  Equation  (242)  for  0.6<°/X<2.6  in  steps  of 
D/X=0.01,  as  shown  in  Figure  60a  and  60b  respectively.  (Recall  that 
below  D/x=  0.5861,  the  cutoff  frequency  of  TEn  mode  of  The  circular 
waveguide,  there  is  no  return  from  the  short.)  These  time  waveforms 
agree  well  with  our  postulation. 

For  completeness  of  the  canonical  response  waveforms,  the  step  and 
ramp  response  waveforms  of  the  open  and  shorted  finite  circular 
waveguides  are  shown  in  Figures  58  and  59  respectively.  The  smoothing 
effected  by  integration  and,  for  the  ramp  response,  the  emphasis  on  low 
frequencies  is  very  evident. 


RAMP  RES  (  |00t 


RAMP  RESPONSE 
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SECTION  4 


INTERNALLY  LOADED  OPEN  CIRCULAR  WAVEGUIDES 

In  this  section,  we  study  the  effects  of  loading  inside  a  semi- 
infinite  circular  waveguide.  We  choose  this  waveguide  instead  of  the 
finite  guide  because  we  thus  eliminate  a  myraid  of  diffractions  from 
the  termination  of  the  guide.  Furthermore,  the  exact  low  frequency 
Wiener-Hopf  solution  is  available  for  the  semi -infinite  circular 
waveguide.  The  loading  is  placed  at  a  position  five  guide  diameters 
down  from  the  opening  of  the  guide. 

In  particular,  we  are  studying  the  time  domain  on-axis  impulse 
response  of  the  semi-infinite  circular  waveguide  with  or  without 
loading.  A  Fourier  synthesis  procedure  is  used  to  obtain  the  time 
domain  response.  The  frequency  spectrum  consists  of  three  regions:  1) 
low  frequency  (below  cutoff)  region;  2)  high  frequency  asymptotic 
region;  and  3)  a  resonance  region  between  the  high  and  low  frequency 
regions  which  will  be  approximated  by  a  rational  function  fit. 

The  low  frequency  region  was  calculated  for  ^/X=0.01  to  r)/X=0.58  in 
steps  of  0.01.  It  is  calculated  using  equations  and  data  on  the  Wiener- 
Hopf  factorization  functions  given  in  the  report  by  Johnson  and  Moffatt 
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[12].  It  is  therefore  the  exact  Wiener-Hopf  solution.  Since  the 
cutoff  wavelength  for  the  TEn  mode  of  the  circular  waveguide  occurs  at 
°/X=0.586,  this  low  frequency  region  is  identical  for  all  cases  of 
loading,  including  the  absence  of  loading. 

The  high  frequency  solution  is  computed  via  the  methods  of  physical 
optics  for  the  internal  load  and  the  geometrical  theory  of  diffraction.* 
It  extends  to  n/x=2.60.  A  detailed  discussion  is  given  by  Pathak  and 
Huang  [50]. 

For  the  region  between  the  high  and  low  frequency  regions,  a 
rational  function  of  the  form 

l 

F(jka)  =  n  0  An(Jka)n  (252) 

h  m 

1 +1  Bm(jka) 
m=l 


is  used  to  fit  points  in  the  high  and  low  frequency  regions. 

To  ensure,  in  part,  the  correctness  of  the  rational  function  fit, 
the  following  checks  are  made. 

1.  Since  the  load  is  placed  five  diameters  down  from  the  open  end 

of  the  guide,  the  impulse  response  for  0<t/to<20.0  =a^  and 

0  c 

all  cases  of  loading  should  be  the  same  as  the  case  of  no 
loading  in  the  same  time  period.  Furthermore,  the  impulse 
response  for  0<t/to<20.0  should  be  the  same  as  that  of  the 

♦These  computations  were  made  by  Mr.  C.C.  Huang  of  the  ElectroScience 
Laboratory,  Department  of  Electrical  Engineering,  The  Ohio  State 
University,  Contract  No.  F19628-80-C-0056. 
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finite  waveguide  discussed  in  Section  3. 

2.  The  impulse  response  and  step  response  must  satisfy  the  moment 
conditions. 


/  Fi(t)  dt  =  0 

0 

(253) 

00 

I  Fu(t)  dt  =  0 

(254) 

o 

where  Fj (t )  and  Fu (t. )  are  the  impulse  and  step  response 
respectively. 

It  must  be  understood  that  in  the  absence  of  experimental  data  it 
is  not  possible  to  completely  verify  the  calculated  results. 

Experimental  data  to  verify  the  calculations  would  require  either  an 
experimental  pulse  range  or  very  broadband  scattering  data  to  isolate 
the  rim  plus  load  return. 

In  the  following  sections,  (4-1)  to  (4-5),  results  for  the  cases  of 
Open  circular  waveguide,  shorted  open  circular  waveguide,  non-shorting 
disc  in  open  circular  waveguide  and  planar  blade  geometry  in  open 
circular  waveguide  (40  blades  and  2  blades)  are  presented.  These 
configurations  are  shown  In  Figure  61. 

All  the  frequency  data  have  spacings  of  0/X=0.0l  and  the  following 
notation  is  used: 

Fasy (s )dV*=0.64,  2.6]  means  the  asymptotic  frequency  response  for 
0.64<°/X<2.6  in  steps  of  °/X=0.01.  It  is  equivalent  to 
Fa$y(s)[u('Vx-0.64)-u(fl/X-2.61)]  where  u(^/X)  is  the  discrete  unit  step 
function. 
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4.1  Open  Circular  Waveguide 


This  Is  a  well  known  geometry  and  will  also  be  referred  to  as  the 
hollow  cylinder  from  here  on.  Frequency  data  used  are  given  by 


Fhollow{s)  CD/a=0-01*  2*63  =  F1ow(s)  [°/X=0.01,  0.58] 


+  Frat(s)  [°/X=0.59,  0.99]  +  Fa$y(s)  CD/X=1.0,  2.6]  . 

(255) 


Flow(s)  tfie  exac^  Wiener-Hopf  solution  given  by  Johnson  and  Moffatt 
[12].  Fasy(s)  [°/X=1.0,  2.6]  is  the  asymptotic  solution  given  by 
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Lastly,  Fra(;(s)[0/X=0.59,  0.99]  is  the  rational  function  fit  spanning 
the  two  solutions.  Figures  62a  and  62b  show  the  magnitude  and  phase 
respectively  of  Fh0n0w($)  up  to  °/X=1.4.  The  frequency  response  is  a 
slowly  varying  function  and  the  rational  function  fit  is  represented  by 
the  solid  curve  in  Figure  62.  The  impulse  and  step  responses  obtained 
by  Inversion  of  Equation  (255)  are  shown  in  Figures  63a  and  63b 
respectively.  Again  the  negative  impulse  at  the  origin  of  the  impulse 
response  curve  is  removed  by  adding  an  appropriate  positive  constant  to 
the  frequency  response.  This  also  applies  to  all  subsequent  plots. 

When  Figure  63a  Is  compared  with  the  0<t<20to  region  of  Figures  54  and 
55,  it  is  seen  that  agreement  is  excellent. 
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Figure  62b.  Frequency  spectrum  (phase)  of  semi-infinite 
hollow  cylinder. 
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4.2  Shorted  Open  Circular  Waveguide 

A  short  is  placed  5  diameters  down  from  the  open  end  of  the  semi- 
infinite  guide.  An  exact  solution  for  the  contributions  of  the  short  is 
given  by  Equation  (242),  which  we  shall  call  FSh(s).  Therefore,  the 
frequency  data  for  this  geometry  are  constructed  as  follows: 

Fshort(s)C°/*«0.01,2.6]  =  Fhollow(s)[°/A=0.01,2.6]  +  Fsh(s )[°/A=0.6,2.6]  . 

(257) 

Fhollow(s)  is  the  same  as  that  given  by  Equation  (255)  and  Fsh(s)  is  added 
to  ^hol 1 0W ( s )  term  term  f°r  0*6<D/A<2.6  .  Figures  64a  and  64b  show 
the  impulse  and  step  responses  of  the  shorted  open  circular  cylinder  via 
inversion  of  Equation  (257).  It  is  notable  that  for  0.0<t<20to(to=d/c) 
the  impulse  response  is  the  same  as  that  of  Figure  63a  while  for 
t/to>20.,  the  impulse  response  is  the  same  as  that  of  Figure  60a.  This 
suggests  that  we  can  deal  with  the  cylinder  and  the  loading  separately, 
thus  greatly  facilitating  the  computations. 

4.3  Non -Shorting  Disc  in  Open  Circular  Waveguide 

This  is  a  more  complicated  problem  since  the  frequency  response  is  a 
rapidly  varying  function  above  cutoff  of  the  TE j j  morfe  the  circular 
guide.  There  is  no  closed  form  expression  for  the  high  frequency  response 
and  an  asymptotic  solution  (Fdi (s))  is  obtained  numerically  using  GTD  and 
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the  physical  optic  approximation  for  the  disc.*  For  frequencies  below 
cutoff,  F i ow ( s )  in  Equation  (255)  is  used. 

The  frequency  data  used  are  as  follows: 

FdiSc(s)[0/x=0-0l.2-6]  =  F1ow(s)[°/X=0.01,0.58]  +  Frat(s)[0/X=0.59,0.63] 

+  Fd1(s)['Vx=0.64,2.6]  (258) 

The  rational  function  fit  is  needed  as  we  know  that  the  asymptotic 
solution  is  not  satisfactory  for  the  region  immediately  above  cutoff.  A 
number  of  points  from  F]ow(s)  and  Fdj(s)  are  used  as  input  to  a  program 
computing  the  coefficients  of  Equation  (252)  by  a  least  square  error 
algorithm.  This  is  a  trial  and  error  method  and  the  "correct"  solution 
must  satisfy  both  physical  constraints  and  the  moment  conditions. 

Figure  65  shows  the  impulse  response  of  an  "incorrect"  frequency 
solution.  The  waveform  for  0.<t/to<20  does  not  agree  with  the 
corresponding  waveform  for  the  hollow  cylinder.  Figure  63a. 

A  satisfactory  solution  is  given  in  Figure  66  and  Figure  67.  The 
circles  represent  the  solution  for  the  hollow  cylinder,  while  the 
triangles  represent  the  asymptotic  solution.  The  solid  line  of  Figure 
66  shows  how  the  rational  function  fit  behaves  inside  and  outside  the 
unknown  region  0.59<0/X<0.63.  The  squares  in  Figure  67  represent  the 
rational  function  fit.  The  corresponding  impulse  and  step  responses  are 


♦These  computations  are  made  by  Mr.  C.C.  Huang  of  the  ElectroScience 
Laboratory,  Department  of  Electrical  Engineering,  The  Ohio  State 
University  under  Contract  No.  F19628-80-C-0056. 
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Frequency  spectrum  (phase)  for  cylinder  loaded 
with  non-shorting  disc  showing  behavior  of 
rational  function  fit  (solid  line)  inside  and 
outside  the  unknown  region. 
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given  In  Figure  68a  and  68b  respectively.  Again  note  that  for 
0.<t/to<20.0  the  impulse  response  agrees  with  that  of  the  hollow 
cylinder  and  the  return  from  the  disc  starts  only  at  t=20to. 

In  general,  the  high  frequency  asymptotic  frequency  data  is 
excellent  for  ®/X  greater  than  1.0.  For  the  case  of  the  non-shorting 
disc,  we  have  used  the  asymptotic  solution  down  to  ^/X^O^;  Equation 
(258).  The  reason  that  this  is  possible  is  that  the  return  from  the 
disc  is  much  stronger  than  that  from  the  rim  at  these  frequencies.  This 
is  clear  when  we  examine  the  region  D/x>0.6  of  the  semi -infinite  hollow 
cylinder  (Figure  62a)  and  that  of  the  hollow  cylinder  loaded  with  a 
non-shorting  disc  (Figure  67a). 

In  target  identification,  we  would  be  interested  in  the  structure 
that  is  present  inside  the  cylinder.  In  our  present  study  the  loadings 
are  always  placed  five  diameters  down  the  guide  so  that  by  the  time 
t=20to,  the  ringing  from  the  rim  has  died  down  substantial ly.  Also  in 
the  case  of  the  non-shorting  disc  the  return  from  the  loading  starting 
at  t=20to  is  rather  strong.  When  these  two  conditions  are  not  met,  that 
is,  if  the  return  from  the  loading  occurs  when  the  ringing  from  the  rim 
has  not  died  down  enough  (due  to  position  of  loading  inside  cylinder)  or 
if  the  return  from  the  loading  Is  so  small  that  it  has  comparable 
magnitude  as  the  ringing  from  the  rim  at  the  instant,  the  waveform  of 
the  return  from  the  loading  is  masked.  Therefore,  it  is  of  special 
importance  that  the  reflection  from  the  loading  inside  the  cylinder  be 
Isolated. 
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Furthermore,  in  the  process  of  isolating  the  return  waveform  from 
the  loading,  we  can  partly  verify  the  correctness  of  the  frequency 
domain  solution.  A  valid  waveform  would  have  to  be  zero  for  0.<Vto<20. 
and  starts  at  t=20to. 

We  shall  first  obtain  the  frequency  spectrum  for  the  reflection  due 
only  to  the  non-shorting  disc  inside  the  guide  which  we  shall  call 
^disc  only(s)*  Note  that  below  °/X=0.5861  there  is  no  wave  propagating 
inside  the  guide  and  the  spectrum  is  identically  zero  in  this  region. 

F  disc  only ($ )[^/A-0.59,2.6]  = 

Fd1sc(s)C0/*»0.01,2.6]  -  Fho11ow(s)[0/X=0.01,2.6]  .  (259) 

The  impulse  response  obtained  via  inversion  of  Equation  (259)  is 
shown  in  Figure  69.  When  compared  with  Figure  68a,  we  see  that  this 
is  an  excellent  way  to  study  the  effects  of  loading.  It  can  be  applied 
to  experimental  situations  with  other  waveguide  geometries  and  loadings. 

4-4  Planar  Blade  Geometry  in  Open  Circular  Waveguide  (40  blades) 

The  geometry  of  the  loading  is  illustrated  in  Figure  61.  There 
are  forty  blades  around  a  small  circular  disc  positioned  5  diameters 
down  the  guide.  The  frequency  data  are  obtained  in  the  same  way  as  that 
of  the  non-shorting  disc  case  and  are  shown  in  Figure  70a  and  70b. 

The  impulse  response  and  the  step  response  are  shown  in  Figures  71a 
and  71b  respectively.  Note  that  they  are  very  similar  to  the  solution 
for  the  non-shorting  disc  (Figures  68a  and  68b).  This  is  attributed 
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Figure  71b.  Step  response  for  semi -inf Inite  cylinder  loaded 
with  a  40-blade  planar  geometry. 


231 


to  the  highly  symmetrical  property  of  the  geometry  and  the  large 
reflection  area  of  the  loading.  Thus  it  is  not  very  different  from  the 
non-shorting  disc  particularly  as  we  are  using  the  physical  optics 
approximation  for  the  loading. 

4-5  Planar  Blade  Geometry  in  Open  Circular  Waveguide  (2  blades) 

The  geometry  for  this  case  is  shown  in  Figure  61.  We  shall  study 
the  loading  with  the  blades  at  0°,  45°  and  90°  with  respect  to  the 
polarization  of  the  incident  electric  field.  The  frequency  spectrum  for 
the  0°  case  is  shown  in  Figures  72a  and  72b.  The  corresponding 
impulse  and  step  waveforms  are  shown  in  Figure  73a  and  73b.  It  is 
apparent  that  the  return  from  the  loading  is  not  strong  and  is  therefore 
masked  by  the  ringing  from  the  rim.  Here  is  an  example  where  we  would 
like  to  isolate  the  return  due  only  to  loading.  Following  the  procedure 
described  by  Equation  (259)  in  Section  4-3  we  obtained  the  impulse 
responses  for  the  0°,  45°  and  90°  cases  as  shown  in  Figure  74a,  74b 
and  74c  respectively. 

The  modulation  effect  is  primarily  evident  in  the  amplitude  of  the 
waveform.  We  have  used  frequency  data  up  to  (V\=2.6  and  thus  included 
altogether  18  TE  and  TM  modes.  However,  the  TEn  is  the  dominant  mode 
and  its  electric  field  lines  diagram  is  shown  in  Figure  75.  Therefore 
it  is  clear  that  for  the  0°  degree  case,  the  return  is  strongest.  Note 
that  the  area  of  the  small  circular  disc  at  the  center  is  0.09ira2  while 
the  area  of  the  two  blades  together  is  only  0. OlSna^.  Thus  the  modula¬ 
tion  effect  is  not  very  dramatic  since  the  small  disc  plays  a  leading 
role. 
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Figure  72a.  Frequency  spectrum  (magnitude)  for  cylinder 
loaded  with  a  2-blade  planar  geometry. 
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inside  cylidner  with  blade  at  45°  with  respect 
to  polarization  of  Incident  electric  field. 


Impulse  response  due  only  to  2-blade  planar 
geometry  inside  cylinder  with  blades  at  90° 
with  respect  to  polarization  of  incident 
electric  field. 


SECTION  5 


COMPLEX  NATURAL  RESONANCES  OF  OPEN  AND  FINITE  CIRCULAR  WAVEGUIDES 

The  complex  natural  resonances  of  a  scatterer  offer  a  unique  way  of 
identification,  Moffatt  [47].  They  are  also  necessary  for  the  make-up 
of  the  K-pulse,  Kennaugh  [44],  a  somewhat  different  target  recognition 
tool.  Although  the  importance  of  the  complex  natural  resonances  is 
recognized  by  many  researchers,  the  resonance  locations  have  been 
obtained  only  for  a  few  simple  geometries. 

In  this  section,  three  different  approaches  are  used  to  find  the 
complex  natural  resonances  of  the  open  circular  waveguide;  ei genanalysis 
of  time  domain  waveforms,  rational  function  fit  of  frequency  domain  data 
and  complex  roots  of  a  postulated  characteristic  equation. 

The  eigenanalysis  solution  of  min-max  problems  involving  squared 
error  is  described  extensively  in  the  literature,  Guillemin  [40].  It 
has  also  been  shown,  Moffatt  et  al.  [49],  to  be  less  sensitive  to 
noise  than  similar  procedures,  e.g.,  Prony's  method. 

Time  domain  impulse  response  waveforms  of  the  finite  cylinder 
(Figures  54a  and  54b)  are  used  In  the  present  analysis.  Since  we  note 
that  for  0.0<t/to<20.0  (t0=a/c,  a=radius  of  circular  guide,  c=speed  of 
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light),  the  impulse  responses  are  the  same  for  the  open  cylinder  and  the 
shorted  cylinder,  the  result  of  eigenanalysis  in  this  region  is  the  same 
for  both  cases  and  is  shown  in  Figure  76.  The  plot  is  in  the  complex 
ka  (normalized  s)  plane  and  it  is  evident  from  the  plot  that  the  real 
part  of  the  extracted  pole  locations  vary  over  some  range.  The  shaded 
area  in  the  plot  indicates  that  the  pole  can  occur  anywhere  in  the 
region.  This  is  a  common  shortcoming  of  the  eigenanalysis  methods. 

The  uncertainty  of  the  real  part  of  pole  locations  according  to  the 
eigenanalysis  method  (including  Prony's  Method)  is  attributed  to  the 
fact  that  the  exact  order  of  the  system  is  unknown  while  we  look  for  a 
certain  number  of  poles.  Also  it  is  known  that  the  sampling  interval 
affects  the  solution,  Moffatt  et  al .  [49]. 

Figure  77  and  73  show  results  of  eigenanalysis  of  the  region 
24.<Vto<40.  for  the  cases  of  a  cylinder  open  at  both  ends  and  a 
cylinder  with  rear  end  shorted  respectively. 

In  the  second  method,  the  moment  method  solution  for  the  hollow 
finite  guide  shown  in  Figure  53  is  used  in  a  rational  function  fit. 
Frequency  data  are  for  0.02<^/X<0.56  in  steps  of  0.02  (D  is  diameter  of 
guide,  X  is  wavelength).  As  we  noted  earlier,  the  cutoff  frequency  for 
the  TE i i  mode  of  the  circular  waveguide  is  n/X=0.586,  the  frequency  data 
includes  only  the  effect  of  the  rim  and  diffraction  at  the  rear  via 
propagation  on  the  outside  of  the  cylinder. 
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Figure  76.  Poles  extracted  from  impulse  response  from  rim 
finite  waveguide. 
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Figure  77.  Poles  extracted  from  impulse  response  from  rear 


of  finite  waveguide  with  both  ends  open. 


Figure  78.  Poles  extracted  from  impulse  response  from  rear  of 
finite  waveguide  with  rear  end  shorted. 


245 


Im  (ka) 


Thus  a  model  represented  by  the  rational  function 


M  -20ika  N 

l  An,(ka)m  +  e  l  Bn(ka)n 

F(ka)  =  m^o _ _ n=o _  ,  (260) 

1  +  I  Cj(ka)j 
j3l 

is  proposed.  The  coefficients  Am,  Bn  and  Cj  are  real  and  are  determined 
through  a  least  square  error  algorithm.  Then  the  expression  in  the 
denominator  of  Equation  (260)  is  solved  to  obtain  the  complex  natural 
resonances.  The  e-j20ka  term  represents  the  time  delay  of  contributions 
from  the  rear  of  guide.  Two  cases  are  presented  here  in  Table  5. 


TABLE  5.  COMPLEX  NATURAL  RESONANCES  EXTRACTED  USING  MODEL  OF 
EQUATION  (5-1)  AND  MOMENT  METHOD  SOLUTION  OF  BELOW 
CUTOFF  REGION  OF  FINITE  CIRCULAR  GUIDE. 


Set  A 


Set  B 


We  recall  that  the  resonance  locations  in  Figure  76  did  not 
include  the  diffraction  from  the  rear,  while  the  complex  natural 
resonances  in  Table  5  include  both  rim  scattering  and  diffraction  from 
the  rear  of  the  finite  guide.  Comparing  the  imaginary  parts  of  the 
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complex  natural  resonances  of  the  two  results,  it  appears  that  at 
frequencies  below  cutoff,  the  effect  of  rim  scattering  is  strongest  as 
the  inclusion  of  rear  diffraction  in  the  model  of  Equation  (260)  did  not 
alter  the  pole  locations  much. 

When  we  examine  the  denominator  of  each  component  in  the 
polarization  scattering  matrix  (Equations  226,  227,  228  and  229),  we 
note  that  all  contain  the  factor  l-fn2.  The  same  denominator  factor  is 
found  in  the  coefficients  which  define  the  waveguide  modes  coupled  into 
the  guide  and  also  appears  in  the  radiation  coupling  coefficients  for 
waveguide  modes.  We  postulate  therefore  that  the  characteristic 
equation  for  the  open  circular  waveguide  is 

1  -  fn2  =  (l+fn)d-fn)  -  0  (261) 

The  factored  form  of  Equation  (261)  is  interesting  in  that  it  appears  to 
define  even  and  odd  complex  natural  resonances  for  the  structure.  For 
axial  incidence,  the  infinite  summations  in  the  components  of  the 
scattering  matrix  reduce  to  the  n=l  term  and  for  axial  incidence  the 
complex  natural  resonances  excited  are  given  by 

(l+fl)(l-f l)  =  0  ,  (262) 

where 

fi  =  U(k)  .  (263) 

2ka  M+('k) 


247 


Asymptotic  forms  for  the  factorization  functions  (L+,  M+)  have  been 
given  by  Chuang  et  al .  [11]  as  (e"'’wtconvention) 


L+(k)  -s!l£  J  i  -  »  ,  ..  *  - 

/Wa  I  2  /irka  e-1"/4  8irka  e-i,r/2 


(264) 


and 


M+(k) 


3  eW4 


1  R'  ,  ♦  R1?  ^ 

/2ka  1  2  /ika  e"1*^  8nka  e“W2 


(265) 


where 


R  =  y  (-i)mei2mka 
m=i  mfn — 


(266) 


and 


to 

R'  =  l  (i)m  ei^mka  .  (267) 

m=l 


Thus  the  complex  natural  resonances  are  defined  by 


1  8itka  -  4/irVa  e111/4  R  -  iR? 

8  irk  a  -  4/ika  eW*  R'  -  iR'2 


±  1 


0 


(268) 


As  noted  by  Chuang  et  al .  [11],  the  infinite  summations  R  and  R' 
converge  relatively  rapidly  except  when  ka=3*/4  and  a  few  terms  are 
sufficient  to  represent  the  series.  However,  we  seek  those  complex 
values  of  ka  which  satisfy  Equation  (268)  and  each  terms  in  R  and  R' 
adds  additional  zeros.  This  obviously  complicates  any  search  procedure. 
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A  Cauchy  integration  program  developed  by  Singaraju  et  al .  [51] 
has  been  used  to  search  for  the  zeros  of  Equation  (268).  The  zeros  of 
1 -f i=0  and  1+f i=0  are  searched  separately  and  they  are  plotted  in  the 
complex  ka-plane  in  Figure  79  and  Figure  80  respectively.  Four  terms 
each  are  used  in  the  summations  of  R  and  R'.  The  solid  line  is  drawn  to 
show  the  cusp-like  trend  of  the  distribution  of  the  zeros.  Furthermore, 
the  cusps  have  a  period  of  it  and  suggests  a  simple  relationship  of  the 
wavelength  with  the  diameter.  As  is  expected  of  an  asymptotic 
expression,  the  locus  of  the  zeros  becomes  more  regular  as  the  frequency 
increases.  Due  to  the  complexity  of  the  characteristic  equation,  any 
solution  becomes  a  formidable  task. 

The  solution  to  Equation  (268)  has  been  attempted  with  different 
number  of  terms  in  R  and  R'.  A  particular  zero  appears  consistently  at 
around  -0.35+jl.l.  This  is  similar  to  the  pole  (set  8  of  Table  5) 
obtained  via  the  rational  function  model  of  the  below  cutoff  region  of 
the  frequency  spectrum.  It  Is  also  similar  to  the  dominant  pole  of  the 
wire  loop. 

Due  to  the  approximate  nature  of  the  equations  and  the  complexity 
of  the  problem,  the  solutions  are  only  approximate.  However,  it  is 
encouraging  to  note  that  there  is  good  agreement  between  the  three 
methods. 

The  question  of  a  characteristic  equation  and  attendant  complex 
natural  resonances  for  the  open  circular  waveguide  is  not  considered 
settled.  The  topic,  particularly  as  reexamination  of  the  Wiener-Hopf 
solution,  is  beyond  the  scope  of  this  report.  Exact  characteristic 


equations  are  known  for  only  two  scatterer  geometries,  spherical  and 
cylindrical.  The  vector  wave  equation  is  also  separable  for  the 
circular  disc  and  some  effort  toward  extraction  of  the  characteristic 
equation  has  been  made,  Mithouard  and  Hodge  [45].  It  would  now  appea 
that  the  open  circular  waveguide  can  possibly  also  be  treated. 


SECTION  6 


SUMMARY  FOR  PART  II 

Backscatter  calculations  for  the  low  resonance  region  of  loaded 
cavity  structures  (open,  loaded  circular  waveguides)  have  been  obtained 
by  combining  exact  (Wiener-Hopf )  solutions  at  low  frequencies  with 
asymptotic  estimates  (CTO)  at  high  frequencies.  The  spectral  solutions 
were  joined  using  rational  function  approximations  to  span  the  critical 
spectral  region  of  the  first  propagating  mode  for  the  guide.  A  proper 
fit  for  the  rational  function  model  was  obtained  by  satisfying  both 
physical  constraints  and  moments  of  the  corresponding  time  domain 
estimates.  The  basic  approach  is  not  new,  one  of  the  earliest 
utilizations  of  the  time  domain  concept  by  Kennaugh  and  Moffatt  [43]  was 
to  combine  Rayleigh  and  physical  optics  estimates  to  predict  resonance 
region  scattering.  In  this  case,  however,  exact  computations  were  used 
at  low  frequencies  and  the  proper  spectral  region  (below  the  cutoff 
frequency  of  the  first  propagating  mode)  was  known  a  priori. 

In  the  time  domain  the  asymptotic  estimates  could  not  be  utilized 
directly,  the  predicted  time  waveforms  being  noncausal.  It  has  been 
demonstrated  however  that  the  noncausality  of  these  solutions  can  be 
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removed  without  altering  the  correctness  of  the  asymptotic  model  in  its 
proper  spectral  range. 

For  an  open,  unloaded  circular  waveguide,  asymptotic  estimates 
proviae  a  reasonable  estimate  of  backscatteri ng  for  guide  diameters 
greater  than  1.0  wavelength  and  an  excellent  estimate  for  guide 
diameters  greater  than  2.0  wavelengths,  Johnson  and  Moffatt  [12],  It 
has  been  demonstrated  using  physical  constraints  and  moment  conditions 
on  the  time  domain  models  that  when  the  guide  is  substantual ly  loaded, 
i.e.,  a  short  or  multiple  blade  load  geometry,  there  the  asymptotic 
estimates  can  be  extended  to  lower  frequencies  (guide  diameters  of  0.7 
to  0.75  wavelengths).  The  reason  for  this  is  that  the  rim  scattering 
(unloaded  waveguide)  becomes  relatively  unimportant  compared  to  the 
contributions  from  the  load.  The  rim  still  must  be  accounted  for 
however  in  the  vicinity  of  the  cutoff  frequency  for  the  first  (TEn) 
propagating  mode. 

While  the  extension  has  not  been  made,  the  manipulations  of 
spectral  limited  analytical  solutions  for  the  case  of  a  open  loaded 
circular  waveguide  have  clearly  demonstrated  that  in  the  low  resonance 
region  much  more  realistic  jet  engine  geometries  could  be  handled.  The 
approach  would  be  via  moment  method  calculations  of  the  structure  at  low 
frequencies  which  will  not  support  propagating  modes  within  the  unloaded 
cavity  and  asymptotic  estimates  of  the  loaded  cavity  when  only  a  few 
propagating  modes  are  possible.  These  solutions  could  then  be  joined 
using  rational  function  estimates  to  span  the  low  resonance  region  of 
the  aperture  of  the  cavity.  As  in  the  cases  demonstrated  here,  both 


254 


physical  constraints  and  moment  conditions  can  be  applied  to  correct  the 
corresponding  time  waveforms. 


The  dominant  complex  natural  resonances  of  finite  and  open  circular 
waveguides  have  been  extracted  using  both  time  and  frequency  domain 
methods.  An  argument  was  given  illustrating  that  the  complex  natural 
resonance  string  associated  with  the  rim  scattering  (leading)  is 
identical  for  all  finite  circular  waveguides  and  is  independent  of  rear 
terminations  and  of  internal  loading.  This  was  demonstrated  for  the 
extracted  poles.  Therefore  we  can  conclude  that  target  recognition 
procedures  based  on  the  complex  natural  resonances  of  the  aperture  of 
the  jet  engine  intake  (see  Appendix  0)  would  be  independent  of  the  jet 
engine  itself  and  would  not  require  Interrogating  signals  with  spectral 
content  above  the  first  cutoff  frequency  for  the  intake. 

A  postulated  characteristic  equation  for  the  open  circular 
waveguide  has  been  suggested  based  on  the  exact  Wiener-Hopf  solution. 
First  estimates  of  the  complex  natural  resonances  corresponding  to  this 
characteristic  equation  have  been  obtained  using  asymptotic 
approximations  for  the  Wiener-Hopf  factorization  functions.  While  some 
rough  agreement  with  the  poles  extracted  via  time  and  frequency  domain 
methods  has  been  obtained,  additional  analysis  is  needed  to  properly 
Interpret  the  complex  natural  resonances  predicted  by  the  proposed 
characteristic  equation. 
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CONCLUSIONS  ANO  RECOMMENDATIONS 


The  following  conclusions  can  be  drawn  from  the  research  detailed 
in  this  report;  these  are  grouped  below  into  catagories  1  and  2 
corresponding  to  parts  1  and  2  of  the  report,  respectively. 

(la)  The  dominant  scattering  mechanisms  have  been  identified  and 
isolated  in  the  frequency  domain  analysis  of  the  problem  of 
electromagnetic  fields  backscattered  from  the  simplified  inlet 
geometries  of  Figures  1  and  2.  Junctions  pertaining  to  these  various 
dominant  scattering  mechanisms  are  defined,  and  the  scattering  from 
these  junctions  is  described  in  terms  of  appropriate  "generalized" 
scattering  matrices.  The  elements  of  the  generalized  scattering 
matrices  alluded  to  above  are  obtained  in  a  relatively  simple  form  via 
the  uniform  GTD  (or  UTD)  ray  analysis  together  with  it's  modifications 
which  are  required  along  ray  caustics.  The  multiple  scattering  method 
(MSM)  has  been  employed,  in  a  self  consistent  fashion,  to  sum  all  the 
interactions  between  the  scattering  junctions.  As  noted  above,  the 
scattering  junctions  are  characterized  by  the  generalized  scattering 
matrices  whose  elements  are  obtained  via  the  ray  method  (see  (lb)) 
together  with  it's  modifications  at  ray  caustics.  The  accuracy  of  the 
elements  of  the  generalized  scattering  matrices  has  been  established  by 
comparison  with  the  exact  (but  far  more  complicated)  Wiener  Hopf 
solution  for  these  elements,  and  also  by  comparison  with  available 
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measurements  in  some  cases.  Consequently,  it  is  concluded  that  the 
solution  to  the  complete  problem  of  the  electromagnetic  backscatter  from 
the  simple  inlet  model  of  Figures  1  and  2,  which  can  be  "built-up"  from 
the  use  of  the  generalized  scattering  matrices  via  the  self  consistent 
MSM  based  procedure,  is  also  expected  to  be  reasonably  accurate. 

(lb)  The  modulation  of  the  backscattered  electric  field  due  to  the 
blade  rotation  is  also  predicted  for  the  simple  inlet  model  in  Figures  1 
and  2.  It  is  concluded,  that  this  modulation  can  be  significant  in  some 
cases  for  the  simple  inlet  model  chosen  in  the  present  work. 

(2a)  The  diagnostic  and  interpretation  potential  of  the  canonical 
(impulse,  step,  and  ramp)  response  waveforms  of  cavity  and  loaded 
cavity-type  structures  have  been  vividly  demonstrated.  For  a  finite 
circular  waveguide  for  example  all  of  the  scattering  me-.hanisms  have 
been  recognized  and  in  many  cases  resolved  in  the  time  domain. 

2(b)  It  has  been  conclusively  demonstrated  that  scattering  data  for 
cavity  structures  in  the  resonance  region  (cavity  ci rcumference  of  the 
wavelength)  of  the  cavity  aperture  can  be  obtained  by  properly  combining 
low  frequency  calculations  and  high  frequency  asymptotic  estimates.  The 
approach  is  particularly  attractive  because  it  provides  a  means  to 
combine  low  frequency  data  where  the  only  fields  in  the  cavity  are 
evanescent  modes  with  solutions  provided  for  example  by  the  geometrical 
theory  of  diffraction.  Computations  at  the  low  frequency  end  are 
Inevitably  simpler  if  propagating  modes  in  the  structure  need  not  be 
considered. 
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2(c)  It  has  been  shown  that  the  non-causal  response  waveforms 
predicted  by  certain  asymptotic  results  for  finite  and  semi-infinite 
guides  can  be  corrected  without  altering  their  correct  form  at  high 
frequencies. 

2(d)  A  definite  procedure  using  rational  functions  to  join  the  low 
frequency  calculations  and  high  frequency  asymptotic  results  was 
established  for  the  case  of  loaded  and  unloaded  semi -infinite  circular 
waveguides.  The  rational  functions  are  simple  in  form  and  both  physical 
constraints  and  moment  conditions  on  the  time  response  waveforms  have 
been  given  to  assure  the  relative  correctness  of  the  rational 
functions. 

2(e)  A  form  for  the  characteristic  equation  for  a  semi-infinite 

\ 

circular  waveguide  has  been  postulated,  and  some  initial  estimates  of 
the  complex  natural  resonances  made  based  on  asymptotic  formulas.  The 
postulated  characteristic  equation  came  from  an  examination  of  the 
elements  of  the  polarization  scattering  matrix  for  the  exact  Wiener-Hopf 
solution,  and  also  from  an  examination  of  the  coefficients  for  the  modes 
coupled  into  the  guide. 

The  following  recommendations  are  made  for  future  studies  on  the 
subject  of  electromagnetic  backscatter  from  inlet  geometries. 

1(a)  It  is  recommended  that  the  present  frequency  domain  analysis 
be  extended  to  treat  non-circular  inlet  rim  shapes. 

1(b)  It  would  be  worth  analyzing  the  effect  of  a  hub  structure  in 
front  of  the  blades  in  the  simple  inlet  geometries  of  Figures  1  and  2. 

It  would  also  be  worthwhile  to  consider  more  realistic  blade  structures 
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(esg.,  those  which  include  a  pitch  on  the  blades)  in  the  inlet 
geometries  of  Figures  1  and  2. 

1(c)  It  is  recommended  that  the  present  analysis  be  extended  to 
treat  higher  frequencies  where  a  large  number  of  modes  can  be  excited 
within  the  inlet  duct.  Alternative  representations  for  the  high 
frequency  fields  within  the  duct  regions  should  be  sought  because  the 
use  of  a  large  number  of  modes,  or  a  large  number  of  equivalent  rays, 
would  lead  to  a  very  cumbersome  solution.  Also,  it  is  recommended  that 
the  effect  of  tapers  or  transitions  in  the  inlet  duct  shapes  be  studied 
since  actual  jet  inlets  do  not  maintain  a  constant  cross-section  within 
the  duct  regions. 

2(a)  It  is  recommended  that  the  rational  function  procedure  for 
obtaining  resonance  region  scattering  data  be  extended  to  finite  cavity 
structures  more  closely  modelling  actual  jet  engine  configurations.  In 
this  case  one  would  be  combining  low  frequency  electric  field  integral 
equation  calculations  with  high  frequency  asymptotic  estimates. 

2(b)  The  initial  research  on  the  characteristic  equation  and  the 
related  complex  natural  resonances  of  a  semi-infinite  circular  waveguide 
should  be  completed. 

All  of  these  recommended  problems  in  l(a)-2(b)  are  very  difficult 
and  challenging,  but  are  certainly  important. 
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APPENDIX  A 

SUMMARY  OF  CIRCULAR  WAVEGUIDE  MOOES 
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APPENDIX  B 


INCLUSION  OF  ALL  MULTIPLE  INTERACTIONS 
ACROSS  THE  APERTURE 

As  shown  in  Figure  16,  the  incident  plane  wave  field,  after  double 

_i 

diffraction,  gives  rise  to  a  non-ray  optical  field  E3  on  the  rim,  which 
was  described  by  Eqs.  (91)  and  (94). 


e3  3  E3,l  +  E3,2  •  (A-l) 

_i  _i  '  ' 

where  ^nd  £3^  are  the  two  ray  optical  components  given  in 

Eqs.  (95)-{96).  Now,  the  triply  diffracted  ray  which  traverses  the 

_i  _bs4 

aperture  to  produce  E4  (and  therefore  E  )  lies  again  on  the 

i 

reflection  boundary  of  the  doubly  diffracted  ray  (which  produces  E3  and 
bs3 

therefore  E  ).  Accordingly,  one  can  decompose  the  non-ray  optical 
i 

E4  into  four  ray  optical  components  as  follows. 


_i  i  1  i  2  i  1  i  2 
E4  =  r4>i  +  Uti  +  E4>2  +  E4,2 


_i  _i 

due  to  E3  ^ 1  due  to  E3 ^ 2 


(A-2) 
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_i  1  _i  2  _i 

Here,  £4^1  and  E4 #  1  are  the  two  ray  optical  fields  scattered  hy  £3^; 

i  1  _i  2  12 

[42  and  £4^2  are  the  two  ray  optical  fields  scattered  hy  E3  p. 

i  1  *  i  2 

Proceeding  the  same  as  in  Eqs.  (95)-(96),  one  can  write  I4 ^  ,  T4  1  , 
_i  1 

E4,2  ,  and  E4j2 


_i  1 

_i 

E4,l 

■  E3,l  ' 

_i  2 

_i 

E4,l 

a  1  F 

7  3,1 

_i  1 

_i 

E4,2 

=  E3,2  * 

_i  2 

_i 

E4,2 

:  1  f 

7  3,2 

-j2ka 


where 


2/tt 


r  r 
PI  P2 


-j2ka 


(pr+2a)(pr+2a) 
-  -j2ka 


r  r 
P3  P4 


(p^+2a)(p^+2a) 


r  r  r 

p=-a  ;  p^=-a  ,  P£=2a  ;  p^-a  ,  P^=4a 


_  I 

Thus,  £4  can  he  rewritten  as 

_1  _i  _i  _1 
E4  *  E4,l  +  E4,2  +  E4,3 


(4-3) 


(A— 4 ) 


(4-5) 


(4-6) 


(4-7) 


(4-8) 


(4-9) 


where 


_il  _i2 


r1  =  r  +  t 

4,1  4,1  4,1 


_i  1 


=  t,  *  o. 


.  1 


-j2ka  . 


_i  _  r-12 
%,2  '  E4,l 


_i 


•ie3,.  i-1’  •  i  •  yF 


r  P-j2ka 


<.3  • 


*  |  E3,2  (_1)  ‘  J  ‘  JT 


S-J‘2ka 


Proceeding  the  same,  one  is  able  to  obtain 
m-1  i 


E1  =  £  T  m>3  . 

m  n=i  m,n 


where 


-1  -1  J  rr  0-J2ka 

E.,l  ‘  Em-1  •  "s  •  J  •  /t?  ' 
-j(2ka-I.)' 


_1 


“m-1 


4/irka 


Lm-1 


and 
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/ 


(A-10) 


(A-ll) 


(A- 12 ) 


(A-13) 


(A-14) 


_1  _i  i - 

E  =  1  E  ,  ,  (-1)  •  j  •  /n-1 

m,n  7  m-l,n-l  v  J 


,-j2ka 


_1 


=  E. 


i-l, n-l  •  j®  •  2  ‘  " 


<  m-1 


(A— 15 ) 


with 


A  = 


4 /irk  a 


(A— 16 ) 


B  =  4.  e-j2ka 


(A-17) 


_ i  _i 

From  Eqs.  (107)- (108),  one  can  try  to  sum  up  £  Em  in  terms  of  £2  to 

m=2 

~  ® 

obtain  a  multiplication  factor  f,  then  [  Ft>sm  is  just  the  product  of 

m=2 

^52  anc|  f#  Now, 


«  4  00 

l  £  -  l 

m=2  m  m=2 

00 

m-1 

l 

n=l 

m-1 

_i 

E 

m,n 

r1 

“  i 

-  I 

I 

+  A  • 

l  r  , 

m=3 

n=2 

m,n 

m=3 

=  jB  • 

00 

l 

m-1 

l 

r1 

/  n  + 

m=2 

n=l 

m,n 

\ /  H+r 

(A-18) 


Thus, 
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(1-A)  I  F 

m=2  m 


®  m-1  i  _  i 

jb  .  i  if  •  [J~+r  . 

m=2  n=l  m»n  ^  n+1  2 


(A— 19) 


Proceeding  the  same  for  K  times,  one  can  obtain 


[1-A 


K  t- 1 

I.  (JB) 


1=1 


CD 

•/I]  I 

v  i  m= 


E 

2  m 


»  (JB) 


“  m-1  i  , _  K  o_i 

I  Z  r  •  rr.  [  l  «B)  /l] 

m=2  n=l  m.n  ^/n+K  l*=l  1  \J  J 


_i 

F2 


(A-20) 


However,  | B  j  <1 .  In  the  limiting  case  where  K  the  first  term  on 
the  right  hand  side  of  Eq.  (A-20)  can  be  neglected.  Thus,  Eq.  (A-20)  is 
reduced  to 


00 


l 

m=2 


_i 

E 

m 


f 


(A-21) 


where 


f  = 


C 

TTC 


(A-22) 


and 


C  = 


K 


l 

*=1 


K  ♦  »  . 


From  Eqs.  (107)-(108),  one  can  write 


(A-23) 
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(A-24) 


I  rt>S"1  (6=0)  .  t”52  (6=0)  •  ?  . 

m=2 

Combining  Eqs.  (87)  and  (A-24),  one  has 

l  rbsm  (0=0)  =  ^Sl  (9=0)  •  { 1  +u 2 * } 
m=l 


(A-25) 
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APPENDIX  C 


TO  DETERMINE  [S21]  FROM  [Si?]  USING 
THE  RECIPROCITY  THEOREM 

As  shown  in  Figure  26,  one  wants  to  find  the  coupled  fields  inside 
the  waveguide  due  to  an  electric  source  7e  in  the  far  field.  Let  the 
electric  and  magnetic  test  sources  "J$,  M$  be  located  at  z=-L,  then, 

_  *f 

Js  =  n  x  Hn 

_  _ +  * 

MS  =  En  x  n 


__+  _+ 

En  and  Hn  represent  the  incident  modal  fields  in  the  waveguide,  and 


E  =  c  fe  +  e  1  e 
n  n  k  n  zn' 


H  =  c  (h  +  h  1  e 
n  n  v  n  znJ 


-j6nz 


(A-27) 


where  en  (or  hn)  is  the  transverse  electric  (or  magnetic)  field,  and 

+ 

ezn  (or  hzn)  is  the  longitudinal  electric  (or  magnetic)  field.  Cn  can 
+  + 

be  either  Anm  or  Bnm  to  denote  the  TE  or  TM  case.  The  test  sources 
_  _  _r  __r 

Js  and  Ms  will  generate  E  ,  and  H  in  the  far  field  which  are 

_  _fi 

assumed  to  be  known  (see  Eqs.  (131 ) - ( 1 33 ) ).  Now,  let  ue  generate  E 
_e 

and  H  internal  and  external  to  the  guide.  From  reciprocity,  one 
obtains 
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//  \fx  TTe  -  Fe  X  / ]  •  n  ds  =  ///  7e  *  ^  dV  • 
E+S+S0  V 


(A-28) 


However, 


-  r  a  e 

n  x  F  =  0  ;  n  x  F  =0 

on  S  on  S 


(A-29 ) 


_r  r  e  e 

Also,  (E  ,  TT  )  and  (F  ,  TT  )  satisfy  the  radiation  conditions  on  E. 


Therefore,  Eq.  (A-28)  can  be  reduced  to 


JJ  H6  •  TTS  dS '+  JJ  E*  •  7S  dS  =  Fe  *  ^  (Fp)  . 

so  S0 


(A-30) 


if  Je  =  pe  -«(| R-Rp | ) ,  a  point  source  in  the  far  field.  Rp  is  the 
position  vector  from  the  reference  point  to  point  P.  Also,  we  let 


—  —  —  ji'.*- 

E  =  i  A“  (e  -  e  )  e  p 
p  P  P  zp 


-e  -  J  Bz 

H  =  i  A-  (-h  -  hi  e  p 

p  P  P  zp; 


(A-31) 


In  the  waveguide  region  at  z=-L.  Again,  A"  can  be  either  A"  or  B- 
3  3  p  nm  nm 

to  denote  the  TE  or  TM  case.  Substituting  Eqs.  (A-26)-(A-27),  and 
(A-31)  into  (A-30),  one  obtains 


APPENDIX  D 


RADAR  RECOGNITION  OF  CAVITY  STRUCTURES 

It  was  demonstrated  in  Part  II  of  this  report  that  a  set  of  complex 
natural  resonances  could  be  associated  with  the  rim  or  aperture  of  the 
cavity.  These  resonances  are  Independent  of  any  loading  (engine)  or 
external  termination.  The  rim  or  aperture  scattering  is  most  dominant 
at  frequencies  just  below  the  cutoff  frequency  of  the  first  propagating 
mode.  The  on-axis  radar  cross  section  of  the  rim  (for  a  circular 
waveguide)  just  below  cutoff  is  some  9db  greater  than  its  nominal  value 
at  higher  frequencies  and  is  roughly  equal  to  the  average  cross  sections 
obtained  when  a  short  circuit  terminates  the  guide.  A  target 
recognition  algorithm  based  on  the  most  dominant  aperture  natural 
resonance  (oscillatory  part  -  0.35  guide  diameters  in  wavelengths)  would 
not  require  penetration  of  the  cavity.  If  the  spectrum  of  the  interro¬ 
gating  radar  signal  lies  below  the  cutoff  frequency  for  the  first 
propagating  mode  for  the  cavity  then  direct  exploitation  of  the  two 
(possibly)  complex  natural  resonances  excited  does  not  appear  too 
promising.  The  amplitude  only  approach  using  synthetically  generated 
matched-filter  response  waveforms,  Moffatt,  Rhoads  [48]  would  be  one 
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possibility.  Another  possibility  which  should  be  tested  Is  the  low 
frequency  classification  scheme  discussed  by  Kslenski  and  Lin  in 
Moffatt,  Young,  Kslenski,  Lin  and  Rhoads  [49]. 

If  somewhat  more  sophisticated  Interrogating  signals  are  considered 
then  two  procedures  can  be  suggested;  prediction-correlation  Chuang, 
Moffatt  [38]  and  the  K-pulse,  Kennaugh  [44]  both  of  which  would  be 
based  on  the  aperture  complex  natural  resonances  of  the  cavity  or  a 
combination  of  these  and  the  load  resonances  (see  below).  The  K-pulse 
seems  particularly  attractive  In  that  very  distinctive  waveforms  should 
result.  Briefly,  the  K-pulse  relates  a  single  excitation  Invariant 
Interrogating  waveform  to  a  unique  family  of  excitation-dependent 
(aspect  and  polarization)  response  waveforms.  The  integrating  K-pulse 
spectrum  has,  as  zeros,  the  complex  natural  resonances  of  the  scatterer 
and  Is  defined  to  elicit  a  response  waveform  of  minimum  duration, 
Kennaugh  [44]. 

Consider  the  solution  for  the  backscattered  field  from  a  loaded 
open  circular  waveguide  using  the  generalized  scattering  matrix 
technique  (GSMT),  Johnson,  Moffatt  [12].  The  scattered  field  is  given 
by 

’  ’ 

Ubs  =  Sn  +  S12  [(Tb2SbT2b)_1  -  S223"1  S21  u1  , 

(A-39) 

where  Sjj,  Si2>  $21  *nd  S22  are  matrices  representing  respectively  the 
rim  scattering,  radiation  characteristics  of  the  rim,  coupling 
characteristics  of  the  rim  and  reflection  of  waveguide  modes  at  the  open 


(rim)  end  of  the  guide  and  u*  is  the  incident  field.  T?b  and  Tb2  are 
diagonal  matrices  representing  respectively  transmission  down  the  guide 
(rim  to  load)  and  transmission  back  from  the  load  (load  to  rim).  Sb  Is 
a  matrix  for  the  reflection  of  waveguide  modes  from  the  load.  The 
matrices  $n»  Sl2»  $21  an<*  ^22  have  been  found  using  the  exact 
Wiener-Hopf  solution,  Johnson,  Moffatt  [12].  Based  on  these  matrices. 
Part  II  of  this  report  suggests  the  characteristic  equation 


[1  -  fn2l  »  [1  +  fn3Cl  -  fid  *  0.  (A-40) 

for  the  aperture  or  rim  complex  natural  resonances .  The  point  is  that 
each  element  in  the  matrices  $n,  Sj2»  S21  and  $22  has  (*  -  fn2)  as  a 
multiplication  factor  in  the  denominator.  If  we  write  these  matrices 
as 


$i j  * - lii-,.-  , 

[1  -  fn*] 


(A-41) 


then 


.bs  J  S11  +  S12 


S 


•  n  -1 


(T  C  T  1"^  _  22 

(b2b2b) 


c  * 
*21 


(A-42) 


or 


bs  1 


,T  s  .  ,-i  .  h? 


n-  J 


-1 


c» 


(A-43) 
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Equations  (A-42)  and  (A-43)  Illustrate  on  Important  point  which  it  is 
felt  will  be  true  for  all  cavity  structures.  The  K-pulse  spectrum 
desired  is  the  square  of  the  aperture  pole  spectrum  as  zeros.  Such  a 
spectrum  still  does  not  contain  the  poles  of  the  load,  l.e.,  the  poles 
affected  by  (Tt)2SbT2b)“**  A  K-pulse  which  included  these  would 
apparently  be  quite  long.  An  excitation  with  the  square  of  the  aperture 
pole  spectrum  as  zeros  however  would  result  in  response  waveforms  very 
highly  characteristics  of  the  load.  From  Equation  (A-42)  it  Is  clear 
that  a  K-pulse  with  the  zeros  of  l-fn^  alone  would  not  eliminate  ringing 
in  the  response  due  to  the  rim. 

The  discussion  at  this  time  is  somewhat  academic  since,  as 
discussed  in  Section  5,  of  Part  II,  we  have  not  yet  been  completely 
successful  in  extracting  the  zeros  of  l-fn2.  However,  development  of 
the  K-pulse  concept  is  In  its  infancy.  The  K-pulse  waveforms  for  a  few 
representative  scatterers  Including  the  sphere,  circular  disk  and 
circular  loop  are  presently  being  studied  on  other  programs*.  We 
envision  however  that  applications  of  the  K-pulse  concept  will  grow  in 
much  the  same  way  as  those  of  the  Impulse  response  concept,  Kennaugh, 
Cosgrlff  [42]  and  complex  natural  resonances  for  target 
Identification,  Moffatt,  Mains  [47]. 
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